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Preface 


This report presents DSN progress in flight project support, tracking and data acquisi- 
tion (TDA) research and technology, network engineering, hardware and software imple- 
mentation, and operations Each issue presents material in some, but not all, of the 
following categories in the order indicated 

Description of the DSN 

Mission Support 

Ongoing Planetary/Interplanetary Flight Projects 
Advanced Flight Projects 

Radio Astronomy 

Special Projects 

Supporting Research and Technology 
Tracking and Ground-Based Navigation 
Communicjtions-Spacecraft/Ground 
Station Control and Operations Technology 
Network Control and Data Processing 

Network and Facility Engineering and Implementation 
Network 

Network Operations Control Center 
Ground Communications 
Deep Space Stations 
Quality Assurance 

Operations 

Network Operations 
Network Operations Control Center 
Ground Communications 
Deep Space Stations 

Program Planning 
TDA Planning 

In each issue, the part entitled “Description of the DSN” describes the functions and 
facilities of the DSN and may report the current configuration of one of the seven DSN 
systems (Tracking, Telemetry, Command, Monitor & Control, Test & Training, Radio 
Science, and Very Long Baseline Interferometry) 

The work described in this report series is either performed or managed by the 
Tracking and Data Acquisition organization of JPL for NASA 
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Network Functions and Facilities 

N A Renzetti 

Office of T racking and Data Acquisition 


The objectives, functions, and organization of the Deep Space Network are summar- 
ized, deep space station, ground communication, and network operations control capabil- 
ities are described 


The Deep Space Network was established by the National 
Aeronautics and Space Administration (NASA) Office of 
Space Tracking and Data Systems and is under the system 
management and technical direction of the Jet Propulsion 
Laboratory (JPL) The network is designed for two-way com- 
munications with unmanned spacecraft traveling approxi- 
mately 16,000 km (10,000 miles) from Earth to the farthest 
planets and to the edge of our solar system It has provided 
tracking and data acquisition support for the following NASA 
deep space exploration projects Ranger, Surveyor, Manner 
Venus 1962, Mariner Mars 1964, Mariner Venus 1967, Mariner 
Mars 1969, Mariner Mars 1971, and Mariner Venus-Mercury 
1973, for which JPL has been responsible for the project 
management, the development of the spacecraft, and the con- 
duct of mission operations, Lunar Orbiter, for which the 
Langley Research Center carried out the project management, 
spacecraft development, and conduct of mission operations, 
Pioneer, for which Ames Research Center carried out the 


project management, spacecraft development, and conduct of 
mission operations, and Apollo, for which the Lyndon B 
Johnson Space Center was the project center and the Deep 
Space Network supplemented the Manned Space Flight Net- 
work, which was managed by the Goddard Space Flight Cen- 
ter The network is currently providing tracking and data 
acquisition support for Helios, a joint U S /West German pro- 
ject, Viking, for which I.angley Research Center provided the 
Lander spacecraft and project management until May, 1978, at 
which time project management and mission operations were 
transferred to JPL, and for which JPL provided the Orbiter 
spacecraft, Voyager, for which JPL provides project manage- 
ment, spacecraft development, and is conducting mission oper- 
ations, and Pioneers, for which the Ames Research Center 
provides project management, spacecraft development, and 
conduct of mission operations The network is adding new 
capability to meet the requirements of the Galileo mission to 
Jupiter, for which JPL is providing the Orbiter spacecraft, and 
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the Ames Research Center the probe In addition, JPL will 
carry out the project management and the conduct of mission 
operations 

The Deep Space Network (DSN) is one of two NASA 
networks The other, the Spaceflight Tracking and Data Net- 
work (STDN), is under the system management and technical 
direction of the Goddard Space Flight Center (GSFC) Its 
function is to support manned and unmanned Earth-orbiting 
satellites The Deep Space Network supports lunar, planetary, 
and interplanetary flight projects 

From its inception, NASA has had the objective of con- 
ducting scientific investigations throughout the solar system It 
was recognized that in order to meet this objective, significant 
supporting research and advanced technology development 
must be conducted in order to provide deep space telecom- 
munications for science data return in a cost effective manner 
Therefore, the Network is continually evolved to keep pace 
with the state of the art of telecommunications and data 
handling It was also recognized early that close coordination 
would be needed between the requirements of the flight pro- 
jects for data return and the capabilities needed in the Net- 
work This close collaboration was effected by the appoint- 
ment of a Tracking and Data Systems Manager as part of the 
flight project team from the initiation of the project to the 
end of the mission By this process, requirements were identi- 
fied early enough to provide funding and implementation in 
time for use by the flight project in its flight phase 

As of July 1972, NASA undertook a change in the interface 
between the Network and the flight projects Prior to that 
time, since 1 January 1964, in addition to consisting of the 
Deep Space Stations and the Ground Communications Facil- 
ity, the Network had also included the mission control and 
computing facilities and provided the equipment in the mis- 
sion support areas for the conduct of mission operations The 
latter facilities were housed in a building at JPL known as the 
Space Flight Operations Facility (SFOF) The interface change 
was to accommodate a hardware interface between the sup- 
port of the network operations control functions and those of 
the mission control and computing functions This resulted in 
the flight projects assuming the cognizance of the large 
general-purpose digital computers which were used for both 
network processing and mission data processing They also 
assumed cognizance of all of the equipment in the flight 
operations facility for display and communications necessary 
for the conduct of mission operations The Network then 
undertook the development of hardware and computer soft- 
ware necessary to do its network operations control and moni- 
tor functions in separate computers A characteristic of the 
new interface is that the Network provides direct data flow to 
and from the stations, namely, metric data, science and engi- 


neering telemetry, and such network monitor data as are 
useful to the flight project This is done via appropriate ground 
communication equipment to mission operations centers, 
wherever they may be 

The principal deliveries to the users of the Network are 
carried out by data system configurations as follows 

(1) The DSN Tracking System generates radio metric data, 
l e , angles, one- and two-way doppler and range, and 
transmits raw data to Mission Control 

(2) The DSN Telemetry System receives, decodes, records, 
and retransmits engineering and scientific data gener- 
ated in the spacecraft to Mission Control 

(3) The DSN Command System accepts spacecraft com- 
mands from Mission Control and transmits the com- 
mands via the Ground Communications Facility to a 
Deep Space Station The commands are then radiated 
to the spacecraft in order to initiate spacecraft func- 
tions in flight 

(4) The DSN Radio Science System generates radio science 
data, l e , the frequency and amplitude of spacecraft 
transmitted signals affected by passage through media 
such as the solar corona, planetary atmospheres, and 
planetary rings, and transmits these data to Mission 
Control 

(5) The DSN Very Long Baseline Interferometry System 
generates time and frequency data to synchronize the 
clocks among the three Deep Space Communications 
complexes It will generate universal time and polar 
motion and relative Deep Space Station locations as 
by-products of the primary data delivery function 

The data system configurations supporting testing, training, 
and network operations control functions are as follows 

(1) The DSN Monitor and Control System instruments, 
transmits, records, and displays those parameters of the 
DSN necessary to verify configuration and validate the 
Network It provides the tools necessary for Network 
Operations personnel to control and monitor the Net- 
work and interface with flight project mission control 
personnel 

(2) The DSN Test and Training System generates and con- 
trols simulated data to support development, test, 
training and fault isolation within the DSN It partici- 
pates in mission simulation with flight projects 
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The capabilities needed to carry out the above functions 
have evolved in the following technical areas 

(1) The Deep Space Stations, which are distributed around 
Earth and which, prior to 1964, formed part of the 
Deep Space Instrumentation Facility The technology 
involved in equipping these stations is strongly related 
to the state of the art of telecommunications and 
flight-ground design considerations, and is almost com- 
pletely multimission in character 

(2) The Ground Communications Facility provides the 
capability required for the transmission, reception, and 
monitoring of Earth-based, point-to-point communica- 
tions between the stations and the Network Operations 
Control Center at JPL, Pasadena, and to the JPL Mis- 
sion Operations Centers Four communications disci- 
plines are provided teletype, voice, high-speed, and 
wideband The Ground Communications Facility uses 
the capabilities provided by common carriers through- 
out the world, engineered into an integrated system by 
Goddard Space Flight Center, and controlled from the 
communications Center located in the Space Flight 
Operations Facility (Building 230) at JPL 

The Network Operations Control Center is the functional 
entity for centralized operational control of the Network and 
interfaces with the users It has two separable functional 
elements, namely, Network Operations Control and Network 
Data Processing The functions of the Network Operations 
Control are 

(1) Control and coordination of Network support to meet 
commitments to Network users 


(2) Utilization of the Network data processing computing 
capability to generate all standards and limits required 
for Network operations 

(3) Utilization of Network data processing computing 
capability to analyze and validate the performance of 
all Network systems 

The personnel who carry out the above functions are 
located in the Space Flight Operations Facility, where mission 
operations functions are carried out by certain flight projects 
Network personnel are directed by an Operations Control 
Chief The functions of the Network Data Processing are 

(1) Processing of data used by Network Operations Control 
for control and analysis of the Network 

(2) Display in the Network Operations Control Area of 
data processed in the Network Data Processing Area 

(3) Interface with communications circuits for input to 
and output from the Network Data Processing Area 

(4) Data logging and production of the intermediate data 
records 

The personnel who carry out these functions are located 
approximately 200 meters from the Space Flight Operations 
Facility The equipment consists of minicomputers for real- 
time data system monitoring, two XDS Sigma 5s, display, 
magnetic tape recorders, and appropriate interface equipment 
with the ground data communications 
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Viking Orbiter Completion Mission 

R L Gillette 

Deep Space Network Operations Section 


This report covers Viking Spacecraft activities from 1 April 1979 through 30 Novem- 
ber 1979 and continues reporting on DSN Viking Tracking support for the same period 


I. Viking Operation 

A. Orbiters 

The Viking Orbiter 1 (VO-1) spacecraft continued to 
operate normally during this reporting period, collecting and 
returning to Earth weather data and high-resolution Mars 
surface photos, as well as relaying to Earth data from the 
Lander 2 (VL-2) spacecraft Approximately 300 Mars photos a 
week have been returned to Earth during this reporting period 
The Viking Orbiter 2 (VO-2) spacecraft ceased operation on 
25 July 1978 

B. Landers 

The Viking Landers also continued to operate as expected 
dunng this reporting period All Lander 1 essential subsystems 
are healthy as the spacecraft collects imaging and meteorology 
data for weekly transmission to Earth whenever a Deep Space 
Station (DSS) is available All Lander 2 essential subsystems 
are healthy, except for the transmitter which supports 
transmission of telemetry data directly to Earth All data from 
Lander 2 are transmitted to Orbiter- 1 and then relayed to 
Earth using the Orbiter transmitter 

II. Viking Orbiter Completion Mission 

On 6 November 1979, a Viking Mars Orbit Trim Maneuver 
(MOT-20) officially terminated the Viking Survey Mission and 


moved the Viking Project into a new phase of the Viking 
Mission termed the Viking Orbiter Completion Mission The 
Viking Orbiter Completion Mission (VOCM) is the fifth phase 
of the Viking Mission following the primary mission (termi- 
nated 15 November 1976), the Extended Mission (terminated 
3 1 May 1978), the Continuation Mission (terminated 
25 March 1979) and the Viking Survey Mission (terminated 
this reporting period, 6 November 1979) The Orbiter opera- 
tions in the Completion Mission are currently scheduled to 
terminate on 1 February 1980 However, a proposal is 
currently being worked between the Viking Project and NASA 
to continue support until the Orbiter runs out of attitude 
control gas sometime in late 1980 

The objective of the Orbiter Completion Mission is to 
obtain moderate resolution photographic coverage of Martian 
surface areas not previously photographed or inadequately 
covered during the survey mission 


III. Viking Lander Monitor Mission 

On 6 November 1979, along with the start of the Orbiter 
Completion Mission, the Lander 1 Spacecraft became a sepa- 
rate mission of its own, termed the Lander Monitor Mission 
(LMM) The objectives of the Lander Monitor Mission are to 
obtain S-band Ranging data from the surface of Mars 
periodically over a long time span for the conduct of Radio 
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Science and to obtain meteorology and imaging data from the 
surface of Mars periodically over a long time span to monitor 
and dissemmate information relative to any significant changes 
with time 

Throughout the Lander Monitor Mission, the Lander 1 
Spacecraft will be in an automatic-mission mode, operating 
autonomously on the programs that have previously been 
stored in the on-board computers Lander-1 will be repointing 
its high-gain antenna and acquiring and storing imaging 
meteorology and engineering data frequently and will be ready 
to transmit these data to Earth every 7 or 8 days in response 
to a series of commands transmitted from a Deep Space 
Station The Viking Lander 1 Spacecraft should be capable of 
returning Mars data up through 1990 

IV. Radio Science 

The opportunity has arisen to make new radio occultation 
measurements of Mars with Viking Orbiter 1 during the Viking 
Completion Mission Geometrically, the spacecraft as viewed 
from Earth is occulted by Mars every 24 hours for about 
23 minutes Acquisition of data from occultation events in 
November and December will provide new information on 


seasonal atmospheric fluctuations at the 5 km level, the 
correlation of ionospheric plasma temperature with solar 
activity, and improved characterization of several Martian 
topographic features 

During the period April 1979 through October 1979, the 
only Radio Science activity has been the near-simultaneous 
Lander-Orbiter Ranging Experiment 

V. Network Support 

Table 1 shows the DSN tracking support for the Viking 
Mission from April 1979 through November 1979 Trackmg 
support started increasing in July, as expected, with the 
resumption of science acquisition activity on the Orbiter 1 
Spacecraft Prior to July, because of pressure from Voyager 
and Pioneer for Network support, Orbiter 1 was placed in a 
housekeepmg mode with all science instruments powered off 
at the end of the continuation Mission (25 March 1979) Since 
August, there has been a continual decrease in the require- 
ments for tracking support This is due to the reduction m 
overall tracking time required to return science data to Earth 
at the higher telemetry data rates made possible as Mars 
approaches opposition on 26 February 1980 


References 
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Pioneer Mission Support 

G M Rockwell 

Deep Space Network Operations Section 


This article reports on activities of DSN Operations in support of Pioneers 6 through 
11 and Pioneer-Venus (Pioneer 12) for the period October 15 through December 15, 
1979 


I. Introduction 

This article reports on the continuing Deep Space Network 
support of Pioneer Mission Operations Included in this 
article is information on the cruise phases of Pioneers 6-9, 
Pioneer 10, post-encounter coverage of Pioneer 11, Pioneer 12 
orbital operations, and other mission-related activities 

II. Mission Operations and Status 

A. Pioneers 6 through 9 

During this period, all Pioneer spacecraft 6 through 9 
appear to be healthy and functioning nominally Only Pio- 
neer 9 was tracked during this period over DSS 14 (Goldstone, 
California) See Table 1 for tracking time A Pioneer 8 track 
was scheduled, but was cancelled because the required 100-kW 
transmitter (64-meter station) was not available during the 
scheduled track time 

B. Pioneer 10 

The Pioneer 10 spacecraft continues on its course out of 
the solar system at an asymptote of 2 9 degrees above the 
ecliptic All onboard systems are functioning nominally The 
round-trip light time for Pioneer 10 at this time is approxi- 
mately 5 hours, 24 minutes That would put Pioneer 10’s dis- 
tance from Earth at 2,908,900,000 kilometers (1,807,900,000 


miles) Tracking time for September and October is shown in 
Table 1 

C. Pioneer 11 

Following its Saturn encounter, the Pioneer 1 1 spacecraft 
began its trek out of the ecliptic at an asymptote four times 
greater than that of Pioneer 10 (12 6 degrees) and in the 
opposite direction of its sister ship As is the case with Pio- 
neer 10, all systems aboard Pioneer 11 are performing nomi- 
nally Table 1 shows total tracking time for September and 
October 

D. Pioneer 12 

In order to not interfere with Pioneer 1 1 Saturn Encounter 
Operations, the Pioneer-Venus Orbiter (Pioneer 12) was con- 
figured for a minimum support activity period from 19 August 
to 8 September 1979 Following this period, normal orbital 
operations did commence and Radio Science data gathering 
was remstituted All occultation and penapsis activities have 
occurred over the Australian facilities during this reporting 
period These activities are discussed later in this article 

On 1 1 October, a successful Static Phase Error (SPE) test 
was performed on the spacecraft’s receiver 2, clearing up some 
doubt as to its condition 
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Also, during this period, a high-gain antenna (HGA) calibra- 
tion was successfully performed on 17 October, over DSS 63 
(Madrid, Spain) The overall tracking time for Pioneer 12 is 
shown in Table 1 

III. Special Activities 

Following Pioneer 1 1 ’s Saturn Encounter, the DSN began 
to support the occultation activity of Pioneer 12 again on 
9 September This began with 26-meter coverage at DSS 44 
(Australia), and 64-meter coverage was reestablished on 
5 October over DSS 43 (Australia) 

There are three ongoing Radio Science Occultation 
Experiments 

(1) Radio Science Occultation I 

(2) Radio Science Corona Turbulence 


(3) Radio Science Occultation II 
All three were discussed in detail m Reference 1 

The last hard occultation occurred on 9 December 1979, 
but data collection in support of Radio Science requirements 
will continue until 27 December 1979 The next period of 
hard occultation will begin in 1980 


IV. Summary 

At this writing, all Pioneer Spacecraft appear to be in a 
healthy state with no apparent problems that would hinder 
nominal operations The period of this report found all Pio- 
neer Spacecraft in a cruise mode, with the exception of 
Pioneer 12’s orbital activity Therefore, no real significant 
events have occurred outside of normal operational support 


Reference 


1 Howe, T W, “Pioneer-Venus 1978 Mission Support,” in The Deep Space Network 
Progress Report 42-51, pp 19-30, Jet Propulsion Laboratory, Pasadena, California, 
June 15, 1979 
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Table 1 Pioneer tracking coverage 


Month 

Pioneer 

Station 

type 

Number 

of 

tracks 

Tracking time 
(hour mm) 

September 

6-9 

26/34 m 
64 m 

0 

(PN 9) 1 

00 00 
03 00 


10 

26/34 m 

12 

48 30 


64 m 

51 

235 41 


11 

26/34 m 

21 

182 17 


64 m 

90 

749 35 


12 

26/34 m 

47 

345 49 


64 m 

13 

48 03 

October 

6-9 

26/34 m 
64 m 

0 

00 00 
00 00 


10 

26/34 m 

13 

60 35 


64 m 

48 

290 02 


11 

26/34 m 

36 

242 14 


64 m 

7 

44 32 


12 

26/34 m 

35 

225 07 


64 m 

65 

523 48 
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Radio Astronomy 

R M Taylor 

Deep Space Operations Section 


This article reports on the activities of the Deep Space Network in support of Radio 
Astronomy operations during the period September 1979 to December 1979 


I. Introduction 

The Deep Space Network 26 meter diameter and 64 meter 
diameter antennas are used to support experiments in Radio 
Astronomy 

II. Radio Astronomy Operations 

During the period of this report, support for radio astron- 
omy operations increased somewhat even though two 26 
meter stations, one in Spam and one in Australia, were taken 
out of service for upgradmg to 34 meters and the addition of 
X-band downlink capability 

A. Support for NASA-Office of Space Sciences 
Activities 

Support for activities in this category continued at previous 
levels except for the Southern Hemisphere Interferometer 
experiment which requires two stations in Canberra, Australia 
One of the two stations is temporarily not available due to 
the upgrade work mentioned above 

Pulsar Rotation Constancy observations contmue and the 
latest version of the timing and recording software from Spam 
is enroute to the Jet Propulsion Laboratory for pre-operational 
testing at Goldstone, California 


Planetary Radio Astronomy activities during this period 
mcluded observations of Jupiter, Uranus, and the Saturn rings 

Microwave Spectroscopy activity recently picked up with 
five observations, producmg a total of 37 hours support from 
the 64 meter antenna at Goldstone durmg December 

B. Support for Radio Astronomy Experiment 
Selection Activities 

The Radio Astronomy Experiment Selection Panel (RAES) 
is now chaired by Dr David S DeYoung, replacmg Dr Bur- 
bridge who resigned recently Dr DeYoung, a scientist with 
the National Radio Astronomy Observatory, is involved in 
research into the ongin and evolution of extended extragalac- 
tic radio sources, physics of galactic nuclei and quasi-stellar 
objects, evolution of dense stellar systems, solar flares, and 
interplanetary disturbances The other panel members are 
Dr Bernard F Burke of the Massachusetts Institute of Tech- 
nology, Dr W C Erickson of the University of Maryland, 
Dr Yervant Terzian of Cornell University, Dr Barry Turner 
of the National Radio Astronomy Observatory, Dr Marshall H 
Cohen of the California Institute of Technology, and Dr 
Richard M Goldstein of the Jet Propulsion Laboratory 

Currently, twelve experiments have been approved by the 
RAES Panel All of these experiments require the use of the 
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64 meter antennas plus the hydrogen maser frequency and 
timing subsystem The one exception to this is the Quasar 
Patrol experiment which can on occasion use a 26 meter 
antenna when the 64 meter antenna is not available 

Two experiments were supported during this period 

RA 175 This experiment uses the Very Long Baseline 
Interferometry capability of the Deep Space 
Network m conjunction with other non-NASA 
telescopes to observe the galactic object 
1909+04, popularly known as SS433 

RA 176 The 64 meter diameter antennas m Goldstone 
and Spam provided nineteen hours each of sup- 
port to this experiment observing the “twin” 
quasi-stellar objects, 0957+561 A, B The 
object of the observations, conducted simulta- 
neously with non-NASA telescopes in Europe 
and the U S A , is to test the current “gravita- 
tional lens” theory about the source 


These observations were not successful due 
to an error in precessed predicts generation and 
further support has been requested 

The early support of RA 175 by the 26 meter antenna at 
Goldstone and the 64 meter antennas at Goldstone, Spam 
and Australia resulted m the reduction of support require- 
ments for the momtor phase — only the 26 meter at Goldstone 
is required for that phase 

C. Support of Host Country Activities 

1 Australia During this period, host country activity m 
Australia has been almost exclusively Pulsar observations by 
the 26 meter station in Canberra These observations are 
currently bemg supported at the rate of approximately 
10 hours per week 

2 Spam In Spam the very long baseline interferometry 
development activity was supported by both the 26 meter and 
the 64 meter antennas 
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Practical Limitations on Noiseless Optical 
Channel Capacity 

S A Butman, J Katz, and J R Lesh 
Communications Systems Research Section 


Even in noiseless optical channels one must take into account the fact that the time 
resolution available is finite An optimization scheme under the constraint of a given 
information rate (in nats /second) and minimum time slot resolution is presented It is 
shown that system efficiencies in excess of 10’s of nats /photon will be extremely difficult 
to achieve due to fundamental time resolution limitations 


In a recent paper (Ref 1) some aspects of the optical 
channel using an ideal binary quantum detector are consid- 
ered In particular, it was found that using a quantized PPM 
modulation scheme, for example, the channel capacity in 
nats/photon approaches infinity in the noiseless case (except 
for the unavoidable self-noise due to the Poisson statistics of 
the photons themselves) However, the final units that are 
relevant for a communication link design are nats/second The 
purpose of this article is to show that due to physical limita- 
tions on the minimum time-resolution one can achieve, for any 
given capacity in nats/seconds there is an optimum point 
which maximizes the capacity in terms of nats/photon, thus 
minimizing the power requirement of the system 

Basically there are three types of time resolution bounds 
The first and ultimate lower bound comes from the quantum 
mechanical uncertainty principle 

AE-AT>h (1) 

If the photon source is a laser, A E corresponds to the laser 
transition linewidth, l e , AT =* l/Ai> where An is the transition 
linewidth (largest An observed is about 10 13 , which corre- 


sponds to AT ~ 0 1 picoseconds) In any case, it doesn’t make 
much sense to consider for a photon of a certain wavelength (X 
= 2 7 r c/cj) time scales that are shorter than the reciprocal of its 
frequency) Thus, we have the bound 

A T^j- (2) 

OJ 

It is true that we actually count photoelectrons, not pho- 
tons, and the arrival times of the former can be determined 
with practically arbitrary accuracy However, time accuracies 
higher than the characteristic time scales of the physical pro- 
cesses associated with the photons (which are the actual infor- 
mation carriers) are meaningless 

A second, more practical bound comes from technological 
limitations, either due to the time response of the photon 
emitter, photodetector, and electronic circuitry or due to 
computational complexity The third type of limitation comes 
from the finite bandwidth of the optical filter 

(3) 

opt 
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This limit usually lies between the first two mentioned above 

Now let us apply these limits to the optical communica- 
tions problem A model for the channel is shown in Fig 1 We 
have a Q-ary alphabet where each letter corresponds to a pulse 
in a particular time slot out of the possible set of Q slots 
Although not necessarily the obvious choice, this quantized 
PPM scheme has the advantage of being an orthogonal signal- 
ing in a channel without power and bandwidth constraints 
This is also the scheme suggested in Ref 2 

In any case it is clear that the best coding schemes call for 
concentrating the power in time slots which are as short as 
possible (1 e , performance is improved by increasing the peak 
to average power ratio) The transition probability is given by 


Here N s is the source intensity (in photons/second) and AT is 
the time slot The capacity of the channel is 

/ -N A T\ 

C s = \ \ - e s I In Q nats/channel use (5) 
In terms of nats/second we obtain 


C T =y = nats/second (6) 

or equivalently 

~N AT 

_ 1 - e s In Q . , 

C T = — -Q- nats/second (7) 


In terms of nats/photon we obtain 

C C 

q — £ — £ rg\ 

P h number of photons per channel use N A T 
or equivalently 

-N AT 

c ph = ln Q nats /photon (9) 

We will now calculate the channel performance in terms of 
nats/photon under constraints of given minimum time slot 
(AT) and fixed channel throughput rate (in nats/second) In 
this case N s = N S (Q) 

In Fig 2, results of C ph vs Q are shown for throughput 
rates of 10 4 and 10 6 nats/second Values of AT correspond to 
the uncertainty principle AT = 10’ 15 sec (= 1/co for X — 1 8 
pm) and for current technological limitations of AT ~ 10“ 9 
sec It is clearly seen that each curve has an optimum operating 
point which minimizes the energy consumption of the system 
(i e , a pomt where C ph is maximum) For example, the C T = 
10 6 nats/ second pomt is Q — 2 10 3 At this pomt C ph — 6 5 
nats/photon 

Note also that with a huge improvement in AT six orders of 
magnitude (from 10' 9 seconds to 10' 15 seconds) we get a 
much smaller increase (a factor of two to three) in C ph This 
shows that pushing the technological limits (l e , decreasing 
AT) would not gam very much in terms of C ph 

It is important to note that for fixed C T and AT, as we 
change Q, N s (transmitter intensity) also changes, as can be 
seen from Eq (7) For Q > Q OJ)t , N s approaches infinity 
Furthermore, the optimum operating pomt occurs for QAT — 
10/C r This implies a PPM block interval (needed for compu- 
tation of peak to average power ratio) of about 10 bits at all 
C T ’ s 
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The problem of estimating the angular position of a deep space vehicle moving at a 
constant velocity using two rotating tracking stations is considered This article reports on 
an initial phase of analytical studies on the optimal attainable estimation performance 
and associated receiver design Parametric dependence of the optimum attainable 
estimation performance is also studied 


I. Introduction 

Consider the situation in Fig 1 where at a reference time 
t = 0 there is a vehicle at pomt V and two stations at pomts S 1 
and S 2 relative to Earth geocenter 0 The vehicle is assumed to 
be at a distance d 0 from the geocenter 0 at time t = 0 It is also 
assumed to be moving at a constant velocity v relative to 
geocenter 0 for all time m a direction which is at an angle a 
with respect to OK The stations are both assumed to be 
rotating at a constant angular velocity to about 0 At time 
t = 0, station Sj is at an angle e Q and the vehicle at an angle y Q 
with respect to a star reference The angle between the stations 
is denoted by 17 and the distances of the stations from 
geocenter by R t and , as mdicated in Fig 1 

We assume that the vehicle is continuously transmitting a 
signal s(t) Each station receives an additive noise-corrupted 
version of this transmitted signal The received waveforms 
from both stations over a given time interval are then used to 
estimate the unknown vehicle angle 7 0 This article considers 
this estimation problem based on the assumptions that d 0 , v, 


a, T\, co, and e Q are known These assumptions are made so 
that the effect of the rotating stations and the moving target 
on angular position estimation can be studied Future studies 
will take into account imprecise knowledge of these 
parameters 

The estimation problem is defined in Section II below along 
with a discussion of the minimum attamable mean square 
estimation error performance Smce we could not determine 
the optimal estimator, a suboptimal estimator is denved m 
Section III The performance of this suboptimal estimator is 
examined relative to the optimal attainable performance 
denved in Section II In Section IV we examine the depend- 
ence of the optimal attainable performance on angular 
position, station rotation, station distance from geocenter and 
observation time duration The specific case of a sinusoidal 
rangmg signal is considered and numencal computations of the 
optimum attainable angular estimation accuracy are performed 
for several parameter values Conclusions are given in 
Section V 
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II. The Estimation Problem 

Let us first derive the equations for the received waveforms 
at each station The vehicle is assumed to continuously 
transmit a signal s(t) The received waveform at station 5, is 
assumed to be 

y 1 (t) = s[r- 0, (r,7 0 )] +«i(r) (1) 

where n,(?) is additive white Gaussian noise with power 
spectral density TV, and 0, (/, 7 0 ) is effectively the signal 
delay time Similarly, the received waveform at station S 2 is 
assumed to be 

y 2 {t) = s[t- <p 2 (t, 7 0 )] +n 2 (t) (2) 

where n 2 (t) is additive white Gaussian noise with power 
spectral density N 2 and <p 2 (t, y 0 ) is the signal delay time to 
station S 2 The signal delay times 0,(7, 7 0 ) and 0 2 (f, 7 0 ) are 
both determmistic functions of the unknown vehicle angle 7 0 
We assume that 7 0 , «,(f) and n 2 (t) are mutually statistically 
mdependent In Appendix A, the geometry of Fig 1 is used to 
derive the following expressions for 0, and 0 2 (c = velocity of 
light) 

0,0, 7 0 ) = [1 - (v/c) 2 }-' [- iy/cft- 0 d Q cos a)/c 2 
+ vR , cos (7 0 + a - e Q - c ot)/c 2 ] 


The receiver’s function is to estimate 7 0 based on observa- 
tions (y , 0), y 2O)), _ T < t < T 2 , with a goal of minimizing 
the mean square estimation error The minimum mean square 
error (MMSE) estimator is the conditional mean estimator 
when the prior distribution of 7 0 is known In this case, the 
conditional mean estimator is nonlinear Moreover, it appears 
that the problem of determining explicit estimator equations is 
not tractable An approach to overcome this problem is to 
derive suboptimum receivers that can be implemented instead 
In the next section of this paper we consider one method of 
obtammg a suboptimal receiver by using an extended Kalman 
filter estimation approach This approach results in a relatively 
simple receiver structure 

It is also of interest to determine the optimum mean square 
estimation error so that the performance of suboptimal 
receivers can be evaluated Of course, the minimum attainable 
mean square estimation error is of interest by itself Unfortun- 
ately, it appears m this case that the problem of determining 
this optimum performance value is also not tractable How- 
ever, it is possible to obtain lower bounds (Refs 1-3) on the 
minimum mean square estimation error The Cramer-Rao 
lower bound (Ref 1 , p 275) appears to be the most tractable 
to use In the case when 7 0 is an unknown but nonrandom 
parameter, the Cramer-Rao lower bound on the mean square 
estimation error of any unbiased estimator 7 Q is given by 
(Ref 1 , p 275) 


+ \ [1 - ( v/c ) 2 ] 2 [t + (vd Q cos a)jc 2 
- vR , cos (7 q + a - e Q - u>t)/c 2 ] 2 
+ [1 - (v/c) 2 ]' 1 [(d 2 0 +R\ - 2R l d 0 cos( 7 0 


£-[( 7 0 - 7 0 ) 2 ] 





(t - 0, 0, 7 0 )) 



(5) 


-e 0 -c of))/c 2 -f 2 ]| 1/2 (3) 

0 a (f.7 o ) = 0 - OV ^) 2 )" 1 [- (v/c) 2 t- vd Q cos ci)/c 2 
+ VR 2 cos (7 q + a e Q - r? - c ot)/c 2 ] 

+ j (1 - (v/c) 2 ) 2 [t + (vd Q cos a)/c 2 

- vR 2 cos (7 q + a - e Q - r? - wt)/c 2 j 2 

+ (1 - (v/c) 2 ) -1 l(d 2 0 +R 2 2 -2R 2 d 0 cos 

" (T 0 ‘ e 0 ' 77 ' wf ))/c 2 - f 2 ] \ 1 (4) 


where s(t) = ds(t)/dt In the case when y Q is a random 
parameter with known density 7>(7 0 ), the Cramer-Rao lower 
bound on the mean square estimation error of any estimator 
7 0 is given by (Ref 1 , p 275) 


E [(V 7 0 ) 2 ] 
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where the expectation in the right hand side of Eq (6) is with 
respect to the prior distribution of y 0 For a normal 7 0 with 
variance , Eq (6) reduces to 

£[(%- 7 0 ) 2 ] 



These lower bounds will be used m the remamder of this 
article to estimate the performance of the suboptimal esti- 
mator as well as the optimum theoretically attamable 
performance 

•II. The Estimator 

Consider the problem of estimatmg 7 0 in the following 
equivalent state variable formulation Let 7(f) be a variable 
state satisfying 


consideration of suboptimal estimators Since the ultimate 
performance measure of mterest is still mean square error, it 
appears more appropriate to seek approximations of the 
MMSE estimator 

There are numerous ways (Refs 4, 5) of determining such 
suboptimal estimators Our approach will be to adopt one 
version (Ref 4, p 267) of the extended Kalman filter 
algorithm This version is the Kalman filter operatmg on a 
linearization of the observation equations (9) about the state 
estimate The reason for adopting this approach over others is 
its relative simplicity In the nonlinear estimation folklore, the 
extended Kalman filter is regarded as being capable of 
performmg as well as other suboptimal schemes m most 
problems So there is a prion no reason to believe that 
constraining our approach to the extended Kalman filter is 
overly restnctive 

Let y(t) denote the extended Kalman filter estimate of y(t) 
Then a straight forward application of the equations of Ref 4, 
(p 267) shows that 7 (t) satisfies 

^ = -P(t) E (UN) Lv,«- s(t- ^,7(0))] 

1=1 


7(0 = 0 


T(-7\) = 7 0 


( 8 ) 


Then 7 (t)= y 0 , the parameter to be estimated, for all t 
Rewrite Eqs (1) and (2) as 




y(t) = 


s(t- 0J (?,7(0)) + « ! (0 


(9) 


y 2 (t ) 


s(t - <t> 2 (t,y(t))) + n 2 (t) 


■ s(f- 0 ; (f,7(O)) 


9 </>,('. 7 (0) 

by 


( 10 ) 


dP(t) 

dt 



s(t~ <t>,(t,y(t))) 


90 , 0 , 7 ( 0 ) 

by 


( 11 ) 


with initial conditions 


9 K-t x ) = 7 0 


( 12 ) 


So the equivalent problem is to estimate y(T 2 ) based on 
observations of y(t) m the interval [- T v T 2 \ 


P(-T x ) = 0 } (13) 

T o 


As we noted previously, the problem of determining the 
MMSE estimator is not tractable An alternative is to derive a 
suboptimal estimator that approximates the MMSE estimator 
Another alternative is to abandon the MMSE criterion and to 
seek estimators based on the maximum likelihood (ML) or 
maximum a posteriori (MAP) criterion However, it can be 
shown (Refs 1 , 4) that the optimum ML or MAP estimators 
are also not practically implemen table Hence, developmg 
estimators using the ML or MAP cnterion will also require 


where 7 0 and are the pnor mean and variance respectively 
of 7 q (We shall denote 


90 ( (L7 O ) 

a 7 n 


by 


7 0 =7(f) 


90, [A7(0] 
by 
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for simplicity) Also, m Eqs (10), and (1 l),.P(r) represents an 
approximation of the conditional variance of y{t) The 
solution of Eq (11) can easily be shown to be 




s (t - 0,(t,7(t))) 


30 (r, 7 (r))l 2 } 

37 J j 


Rewntmg Eq (10) as an integral equation gives 



y ; (r)- s(t- 0,(7, 7 ( 7 ))) 


30 ( ( 7 , 7 ( 7 )) 

s(t- <p ( 7 , 7 ( 7 ))) — 

ay 


dr 


(15) 


P(t ) is only an approximation to the conditional variance of 
y(t) (Refs 4, 5) Moreover, P{t) depends on the observations 
and so cannot be determined other than from simulation runs 
of the filter In spite of these pitfalls, let us examine Eq (14) 
to obtam a heuristic estimate of the best possible performance 
of the estimator 

Assume that the estimator is performing well Thus, y (t) 
will be close to 7(f) = 7 0 Assume also that P(T 2 ) is a good 
approximation of the mean square estimation error Then, 
from Eq (14) we have 

P(T 2 ) = 


s 2 ( t ~ 0 ; (t,7(O)) 


30 ( ( f . 9 ( f ))\ 2 
~ 99 / ‘ 


Thus, Eqs (14) and (15) give the estimator structure with 
7 (T 2 ), the desired estimate of y Q The only prior statistical 
knowledge of 7 0 required is its mean and variance A block 
diagram of the implementation of the estimator is given m 
Fig 3 The expression for 30, [f, 9(?)] jdy is given in Eqs (B-2) 
and (B-3) of Appendix B These waveforms are implemented 
in the receiver of Fig 3 by adjustmg the y(t) phase contribu- 
tions in the sinusoidal terms given in Eqs (B-2) and (B-3) The 
structure of the estimator is somewhat similar to the MAP 
estimator with normal prior distribution for 7 0 (Ref 1, p 
453) One substantial difference of the estimator here with the 
MAP estimator is that the gam term P(t) in the integral in 
Eq (15) is replaced by 2° 2 y 0 in the MAP estimator This is 
because the gain is updated to account for the change in the a 
posteriori variance of y Q based on the observations This is not 
performed in the MAP estimator 

In Eqs (14) and (15), 0,(r, 7) is given by Eq (3) and 
30, (t, 7)/ 3 7 by Eqs (B-2) and (B-3) m Appendix B Further 
simplification of Eqs (14) and (15) can result from using the 
simpler approximations in Eqs (24) to (27) given in section IV 
for 0, and d<pjdy Simplification of the basic estimator 
structure depicted in Fig 3 apparently cannot be done 
without specific assumptions on the signal structure 

The performance of this algorithm unfortunately cannot be 
determined analytically In evaluating extended Kalman filters, 
P(t) is often regarded as a measure of the mean square error 
However, care must be taken to adopt this conclusion since 


2 r 7 ' 2 , „ / 30 (r,7(0) \ 2 

< £ 0 !N)J S 2 {t- 0,(f,7(O)) ay J dt 

(16) 

Let 7 0 be the true value of the unknown angle So, if 7 (t) = 
7 0 , replacing 7 (f) by 7 0 in Eq (16) shows that the upper 
bound on P(T 2 ) is roughly twice the Cramer-Rao lower bound 
(Eq (5)) on the optimum mean square error Thus from the 
above heuristic point of view, the best possible performance of 
the estimator is roughly within a factor of 2 from the 
Cramer-Rao lower bound of Eq (5) 


IV. Optimum Theoretically Attainable 
Estimation Performance 

As we noted previously the Cramer-Rao lower bound gives 
a lower bound on the optimum attainable angle mean square 
estimation error In this section we shall examine the 
Cramer-Rao lower bound in a special case In particular, we 
shall assume the following set of parameters 


d Q = 8 X 10 8 km 
R l =R 2 = 6 5 X 10 3 km 

(17) 

v = 10 km/sec 
T x =T 2 = 30 mm 
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This set of parameters is consistent with the distances 
encountered in a Jupiter mission We assume m addition that 
N. = N 2 for simplicity We shall first analyze the effects of the 
relative angular positions and the rotation of the Earth on the 
Cramer-Rao lower bound (Eq (5)) This, then, gives the 
dependence of the optimum attainable performance on these 
effects 

We first consider the effect of the angular positions y 0 , e 0 
and 77 given m Fig 1 Since the problem of estimating 7 0 is 
nonlinear, the minimum attainable estimation error would 
generally depend on y g Consider first the case when to = v = 0 
for insight into this dependence Usmg the parameters in 
Eq (17) we have from Eqs (3) and (4) that 

</>,0, 7 0 ) 

= \d 2 0 + R\ - 2 d 0 R 1 cos ( 7o - e„)J "j c 

38 [l-(Vo/(< + ^ c «(V E o)] 

(18) 


<P 2 (t,y 0 ) 


dl +Rl 


2d 0 R 2 cos ( 7 0 



= [(^^‘V) ['-(Vo^ + «2» 

cos ( 7 0 - e 0 - r?)J (19) 


^[(7 0 -7 0 ) 2 ] 


C 2 

N t c 2 (d 2 0 +R 2 t )J- Ti 


t - 


(^o + ^i) 


sm 2 ( 7 0 - e 0 )dt + 


2R 2 2 d 2 0 


N 2 c 2 {d 2 0 +R\) 



sin 2 ( 7 0 - e Q - i i)dt 

( 22 ) 


Hence under the assumption that R^= R 2 and N l = N 2 , Eq 
(23) depends inversely on 


/( 5) = sin 2 6 + sin 2 ( 5 - 7 ?) = 1 - cos 7 ? cos (25 - 7 ?) 

(23) 


where 5 = 7 0 - e 0 Note that /( 5) is symmetric about 5 = 17/2, 
which corresponds to when the vehicle is halfway between the 
two stations (see Fig 1) When 0 < 77 < 90°, /(5) increases as 6 
deviates from 77/2, or when the vehicle moves toward either 
station from the midpoint So, when 0 < 77 <90°, the worst 
performance is when the vehicle is exactly halfway between 
the two stations This is shown in Fig 4 When 77 > 90° , the 
converse is true and the best performance is when the vehicle 
is exactly halfway between the two stations Since /(S) is 
mdependent of 5 when 77 = 90°, this is the best value of 17 from 
the viewpomt of uniformity of performance over a range of 
7 0 An exammation of Eq (23) shows that for 80° < 7 ? < 
100 °, the variation of performance is less than 20 % for 5 from 
0 to 77 


so 


9 7 0 ) 
97 0 


R i d o/ c + R 2 J ' 2 sm (y o " e o) ( 2 °) 


The above considerations are when u> = v = 0 Let us now 
consider when to 0 and v ¥= 0 In Appendix B, it is shown 
that approximate expressions for <j> , <t > 2 , d<p l /dy Q and 
90 2 /97 o are 


0,O,7 o ) 


90 2 (L 7 0 ) 

97 0 


R 2 d o/ c (d l + R \^' 2 sm " e o " 7 ?) ( 21 ) 



vd Q cos a 

[(d 2 0 +R])(c 2 -v 2 )} 112 


t 


Since [Rdjc{d 2 0 + R 2 ) l < 2 ]= 2 17 X 10" 2 , 0,(f. 7 0 ) is 
relatively independent of 7 0 Thus an approximation of the 
Cramer-Rao lower bound (Eq (5)) in the case when to = v = 0 
is 


(d 2 Q +R]\ 112 vd Q cos a 
V C 2 - V 2 / " c 2 - V 2 
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R.d„ 


[(d 2 +/? 2 )(c 2 -v 2 )] 


1/2 


c°s(7 0 ~ e 0 ~ 


(24) 


0 2 ('.7 o ) 


-('-m 


(7) 


vd Q cos a 


[(rf 2 +* 2 )(c 2 -v 2 )] 


1/2 


/< + « 2 \' /2 v<2 0 cos a 

\ c 2 - v 2 / _ c 2 -v 7 


Vo 


{{dl+R^ic 2 


v 2 )] 


1/2 


cos 


(7 0 * e 0 “ 77 " wr ) ( 25 ) 


90,(67 o ) 

97 0 


is of the order 10“ 2 Hence, the second term in Eq (26) can 
be neglected except when sin (y 0 - e 0 + a - c ot) is sufficiently 
larger than sin (7 0 - e 0 - tot) In an extreme case when 7 0 - 
e 0 = 0° and a = 90°, the first term in Eq (26) is zero at t = 0 
However, as t deviates sufficiently from 0, the first term will 
again dominate the second term For example, if \t I = 10 sec, 
the first term is 10 times the second in Eq (26) So, m 
instances when the observation time interval T { + T 2 is much 
larger than 10 sec, the contribution of the second term in Eq 
(26) to the Cramer-Rao lower bound will be negligibly small 
The same conclusion can be drawn for the second term in Eq 
(29) Hence, neglecting these terms results in the following 
approximation to the Cramer-Rao lower bound (Eq (5)) 


£-[(7 0 


To) 2 ! * 


2R]dl 


N x {c 2 - v 2 ) (d 2 + R 2 ) 



(t~ 0j(fi 7 0 )) 


R i d o 


l(d 2 0 +R])(c 2 -v 2 )] 112 


Sln (T 0 " e o “ wf ) 


vR. 


2 2 
C ~ V 




2 2 
C - V 


( JL ) 2 ]' 1 
\ cl ) 


• sin(7 0 -e 0 +a-wf) 


(26) 


(To " e o ” w/ )^ 


2R\d\ 


N 2 (c 2 - v 2 )(d 2 0 +R 2 2 ) 


c 


s 2 (t- <t> 2 {t, 7 „)) 


90 2 (67 o ) 

a T 0 


-1 


sin 2 (7 0 - e 0 - ??- oj t)dt 


(28) 



• sin(7 0 - e Q +a- r?- u>t) (27) 

In Eq (26), the factor in front of sin (7 0 - e 0 + a - cot) is of 
the order 10“ 7 while the factor in front of sin (7 0 - e 0 - cot) 


where <p x and <f> 2 are given by Eqs (24) and (25) respectively 

Let us now compare Eq (28) with Eq (22) when to = v = 0 
From Eqs (24) and (25) it can be seen that the dependence of 
<p and 0 2 on 7 0 is small We may assume that <p l and 0 2 are 
both essentially independent of 7 0 in Eq (28) So from the 
viewpoint of dependence on 7 0 ,the essential difference in the 
structure of Eq (28) to the structure of Eq (22) is the sin 2 
(To ~ e 0 “ w 0 an d sln2 (7o “ e o ~ r >~ factors in the inte- 
grands in Eq (28) versus the corresponding sin 2 (7 0 - e 0 ) and 
sin 2 (7 0 — e 0 - 77) factors in Eq (22) Although the earth 
rotational angular velocity go = 7 27 X 10“ 5 rad/sec, for t - 30 
minutes tot = 7 5° Hence this difference is certainly not 
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negligible This points out a significant contribution to the 
estimation performance due to the rotation of the stations 

To assess the dependence of Eq (28) on the angular 
position 7 0 we assume that 4> l and 0 2 are essentially 
independent of y 0 in Eq (28) Under the assumption that 
TVj = N 2 and = R the integrand m Eq (28) is directly 
proportional to 

sin 2 (7 0 - e Q - tot) + sin 2 (7 0 - e Q - r? - wf) 

= 1 - cos t? cos [2(7 - e Q - wf) - p) (29) 

Comparing Eq (29) to Eq (23), we see that to a first order 
approximation, the conclusions regarding the dependence of 
performance on angular position 7 0 - e Q in the case w = v = 0 
still hold here In particular, it is clear from Eq (29) that from 
a viewpoint of uniformity of performance over a range of 7 0 , 
angular positions near r\ = 90° are desirable 

Let us consider next the effect of varying the parameters 
R j and R 2 on the optimum attainable estimation perfor- 
mance Recall that R l and R 2 are the distances from the 
stations to geocenter In Eq (17), the values of R 1 andf? 2 are 
for ground-based stations The other case of interest is when 
the two stations are in geostationary orbit with R x = R 2 We 
shall consider an orbital radius of up to 10 s km Hence we 
need to examme the dependence of Eq (5) on R x = R 2 for a 
range of these parameters from 6 5 X 10 3 km to 10 s km We 
still take d Q , the distance of the vehicle to geocenter, to be as 
in Eq (17) Hence d Q is still much larger than R 1 and R 2 In 
this case, an examination of the derivation in Appendix B 
shows that we can still use the approximation in Eq (28) to 
Eq (5) with the expressions m Eqs (24) and (25) for <p l and 
02 , respectively Since d 0 is much larger than R j = R 2 , an 
exammation of Eqs (24), (25) and (28) shows that the 
Cramer-Rao lower bound is proportional to l//? 2 = 1/7? 2 In 
other words, the root mean square estimation error is directly 
proportional to 1/7? . = 1 /R Increasing, R^= R 2 from 6 5 X 
10 3 km to 10 s km will decrease the optimum root mean 
square estimation error by two orders of magnitude 

Finally, let us consider the effect of varying the observation 
duration T y + T 2 on the optimum attainable estimation 
performance We assume that T { + T 2 is large compared to 10 
seconds and that the other parameters are given as in Eq (17) 
Then Eq (28) is again a valid approximation of Eq (5) with 
0 and 0 2 approximated by Eqs (24) and (25), respectively 
We also assume that the frequency of the ranging signal s(f) is 
much higher than l/(T + T ) and also much higher than 
w/27t (w = rotational angular velocity of the stations) It is still 
difficult to assess the dependence of Eq (28) on T and T 2 in 


general because of the sin 2 (7 - e Q - cot) and sin 2 (7 - e Q - 
r\ - cot ) terms in the integrals in Eq (28) These terms change 
the value of the integrands as T and T 2 are varied To a first 
order approximation it appears that the right-hand side of Eq 
(28) is inversely proportional to 

(T l + T 2 ) [sin 2 (7 0 -£ 0 -u 7’ 2 )- sin 2 (y Q - e Q + coTJ 

+ sin 2 (7 0 - e Q - i) - coT 2 ) - sin 2 ("y 0 - e Q - 1? + w^)] 

= (T l +T 2 ) \2- cos t? cos [2(7 0 - e Q - w T 2 )~ 77)] 

- cos r? cos [2(7 0 - e 0 + coT ,) - tj] } 

In the case when n = 90°, Eq (30) reduces to 2(7’ 1 + T 2 ) 
Thus, when 17 = 90°, the optimum attainable root mean square 
error performance is approximately inversely proportional to 

y/(r‘ + t 2 ) 

Finally, we consider a specific ranging signal s(t) and 
perform numerical computations of the Cramer-Rao lower 
bound 

Example 

Consider a sinusoidal ranging signal of frequency f Hz 
That is, 

s(t ) = x /2Scos(27r/ c f) 

We assume that for 1 = 1,2, the demodulated ranging signal 
power to noise spectral density ratio is 

— = 10 dB 

N 

l 

This signal-to-noise ratio is consistent with X-band carrier, 
20-dB vehicle antennae gain, 53-dB station antenna gains, 
receiver noise temperatures of 50°K, 20 W vehicle transmitted 
power and a 3-dB modulation loss We also assume that 

d 0 = 8X 10 8 km 
v = 10 km/sec 
/?, =R 2 = 6 5 X 10 3 km 
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These parameters are consistent with that encountered in a 
Jupiter mission with ground-based stations We also assume 
that T x = T 2 Numerical Monte Carlo integration was used to 
compute the value of the Cramer-Rao lower bound for vanous 
values of signal frequency / and observation time duration T x 
+ T 2 The numerical computations are within a 1% accuracy 
These numerical results are summarized in Tables 1 and 2 
below The listed angle estimation accuracies m these tables 
are the square root of the Cramer-Rao lower bound 

Table 1 shows that the optimum angle estimation accuracy 
is inversely proportional to the frequency of the sinusoidal 
ranging signal Although this particular relation between 
estimation accuracy and signal frequency does not hold in 
general, it can be easily seen from Eq (5) that signals of higher 
frequency give a smaller Cramer-Rao lower bound Also note 
that Table 2 shows that the estimation accuracy is approxi- 
mately inversely proportional to T x + T 2 , as we would expect, 
since rj = 90° 

We note that the above angle estimation accuracy was 
obtained using the Cramer-Rao lower bound, Eq (5), which is 
valid when 7 0 is an unknown but nonrandom parameter 
Suppose instead that 7 0 is a random parameter and can a 
priori be assumed to be normally distributed Then the 
relevant lower bound on mean square estimation error is 
Eq (7) We claim that if the a priori variance of y Q is much 
larger than the lower bound Eq (5), then the above estimation 
accuracy calculation is still valid This follows because Eq (5) 


is essentially independent of 7 in this case since 17 = 90° 
Hence, the expectation term in Eq (7) is 1 /(lower bound 
Eq (5)) 

V. Conclusion 

This work has considered the problem of estimating the 
angular position of a moving vehicle using two rotating 
stations The optimum attainable angle mean square estima- 
tion error was derived along with an implementable sub- 
optimal estimation algorithm A situation comparable to that 
encountered in a Jupiter mission was further analyzed In this 
situation it was shown that the optimum angle between the 
two stations from a viewpoint of uniformity of estimation 
performance is 90° It was also shown that the optimum 
attainable estimation accuracy varies inversely with the dis- 
tance of the stations from geocenter and approximately 
inversely with the square root of the observation time 
duration The optimum attainable angular estimation accuracy 
was numerically computed for a sinusoidal ranging signal 
These computations show that the optimum attainable estima- 
tion accuracy is 0 02 qrad for a 2 MHz signal and an 
observation time of one hour 

This work has only considered the problem of estimating 
one angle with the range, velocity and other angles of the 
vehicle known Further work should be done to include some 
or all of these parameters as parameters to be estimated along 
with the angle considered in this work 
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Table 1 Estimation accuracy vs signal frequency 


Signal frequency 

Optimum angle 

fc 

estimation accuracy 

2 MHz 

0 018 inad 

5 MHz 

0 0072 jurad 

10 MHz 

0 0036 qrad 

20 MHz 

0 0018 qrad 


Table 2 Estimation accuracy vs observation time 


Observation time 

Optimum angle 

duration (Tj + T 2 ) 

estimation accuracy 

10 min 

0 046 qrad 

30 min 

0 026 qrad 

60 mm 

0 018 qrad 

90 min 

0 015 /urad 


23 




Fig 1 Vehicle and station locations 



Fig 2 Relative locations for t = 0, t,, t. 


24 





Fig. 3 Estimator block diagram 
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Appendix A 

Derivation of Expressions for fa and fa 


It follows from Eqs (1) and (2) that the signal value 
received at station S at time t is the signal value transmitted 
by the vehicle at time t - <t>j(t, y Q ) for r = 1,2 Define 

r0) = t- <t>'(t,y 0 ), i = 1,2 (A-l) 


For time instants t and t let 

i / 

r,(f,r) = distance between the vehicle at time t and 

l v f y i 

station S, at time t 
i / 

r 2 {t j ,t / ) = distance between the vehicle at time t { and 
station S 2 at time t 


It is then clear from Eqs (1), (2) and (A-l) that 


4>,( f .7 0 ) 


r, 

c 


(A-2) 


= d\(t)+R\-2R x d x (t) 
cos [fo " e o " ~ #(',)] 

r lW =d]{t) + R\-2R 2 d x {t) 

cos [h + e 0 - >0 + Q(t) + cjt) 

Next, using Eqs (A4) and (A-5) we get 

cos 8 ( t ) = [1 - sin 2 0 (^) 1 1/2 
d n + vt cos a 

_ o I 

d(t) 


(A- 6 ) 


(A-7) 


(A- 8 ) 


or equivalently, 


r 


(0 


c 


(A-3) 


Now, substitution of Eqs (A-5) and A- 8 ) into Eqs (A- 6 ) and 
(A-7) yield the following expressions for r 2 l (t p t j ) and r 2 (t F t) 

r\ (t r t') = d 2 (t t ) + R 2 - 2 R l d Q cos ( 7 0 - e 0 - wf ; ) 


for i = 1,2 Hence if Eq (A-3) can be solved for T t (t), then 
0|(t,7 o ) may be obtained using Eq (A-l) In order to 
accomplish this we must first obtain expressions for r ] (f,,^) 
and r 2 (f t ,t^ Smce r ( (r) < t in (A-3), we are interested in the 
cases 0 <t l <t J , t l <0<t ) and t t <t < 0 only in determining 
these expressions Consider first the case when < 0 as 

given in Fig 2 Also we shall assume that 0 <a < 180 in Fig 
2 Let 

d{t t ) = distance between the vehicle and Earth geocenter 
at time t 

l 

Oit') = relative angle at geocenter between the vehicle 
positions at t = 0 and t = t 

From the geometry of Fig 2 we can write the following 
relations 

d 2 {t) = v 2 t 2 +d 2 0 +2d Q vt t cos a (A4) 

d{t) ~(vt) 

— - = ^ (A-5) 

sin a sin d(t) 


- 2R X vt t cos (a + 7 0 - - w t) (A-9) 

r 2 2 (t,d,) = d 2 (t) + R\ - 2 R 2 d Q cos (rj + e Q - y Q + ojt) 

- 2 R 2 vt t cos (a - 77 - e Q + 7 0 - u>t) (A-9) 

So, Eqs (A4), (A- 8 ) and (A-9) define r 2 (t f t f ) and r \(t p t ) for 
the case t i <f / <0and0<a< 180° Arguments similar to 
the above can be used to estabhsh Eqs (A4), (A- 8 ) and (A-9) 
when 0 < t < t or t < 0 < t and also when 180° < a < 

o 'll i 

360 Hence, we may conclude that these equations hold for 
all cases of interest in solvmg Eq (A-3) for T ( (f) For simplicity 
let us denote t ) t) by r f for the remainder of the appendix 
Squaring Eq (A-3) yields 

c 2 (r ( - r) 2 = r 2 (r p t), i = 1,2 (A-10) 

Next, substitution of Eqs (A4), (A- 8 ) and (A-9) mto Eq 
(A-10) yields the following quadratic equations satisfied by r 
and r 2 
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(c 2 - v 2 )r 2 - 2 [c 2 / + vc/ 0 cos a - cos (a + y Q - e - co/)] r 

- [c/ 2 +2? 2 - 27? 1 rf 0 cos(r Q - e Q - w/)- c 2 / 2 ] =0 

(All) 

(c 2 - v 2 )t 2 - 2 [c 2 / + i>c/ Q cos a - vR 2 cos (a - t? - e Q + y Q 

- wt)]r 2 

- [< + ^2 - 2 R 2 d 0 C0S ^ + e 0 " ^0 + 

- c 2 / 2 ] = 0 (A- 12) 

The constraint that r < / identifies the roots of Eqs (A-ll) 
and (A-12) that give t 1 and t 2 respectively These are 

T 1 = (1 - (v/c ) 2 ) 1 [/ + (vd Q cos a) fc 2 - vR l cos (y Q +a - e Q 
- c ot)jc 2 ] - | (1 - (v/c) 2 ) [/ + (vd Q cos a)/c 2 


2 2 

- vR j cos ( 7 0 + a - e Q - co/)/c ] 

+ d - (v/c) 2 )^ 1 [(d 2 +/? 2 - 2 R t d 0 cos( 7 0 - e 0 
-t ot))/c 2 - t 2 }\' 12 (A-13) 

r 2 = (1 - (v/c) 2 )" 1 [/ + ( vd Q cos a)/c 2 - v/? 2 cos ( 7 0 + a - e Q 

- tj - co/)/c 2 ] +| (1 - (v/c) 2 ) 2 [/ + (vrf 0 cos a)/c 2 

- v^ 2 cos ( 7 q + a - e Q - tj - cof)/c 2 | 2 

+ (1 - (v/c) 2 )"' [(^o +/? 2 ' 2R 2 d o cos ^o ' e o - * 

- w/))/c 2 - / 2 ] J 1/2 (A-14) 

Finally, Eqs (A-l), (A-13) and (A-14) give the expressions in 
Eqs (3) and (4) for 0j(/.7 o ) and $ 2 (/.7 0 ), respectively 


28 



Appendix B 


Approximations to <f> u (f> 2 , 


3</> 


9yo 


1 and 


d<t> i 
9yo 


Consider first 0 given by Eq (3) Ignoring terms smaller 
than 10 -9 , we have 


jO - ( y/c ) 2 ) 2 [? + (vd Q cos a)/c - vR t cos ( 7 0 + a - e Q 
- ccf)/c 2 ] 2 +(1 - (v/c) 2 )" 1 [(d 2 + * 2 - 2*^ 

c °s (7 0 ‘ e 0 " wr ))/c 2 - C 2 ] } 1/2 

/c? 2 +/? 2 \'^ 2 2vc? 0 cosa 

_ \c 2 -v 2 / L 1+ (d 2 +J R 2 )(l-(v/c) 2 ) r 


-wr)/c 2 -? 2 ]} ; 

• |(1 - (v/c) 2 )" 2 [ f + K> cosa)/c2 
- v/? J cos (7 0 - e 0 + a - w/)/c 2 ] 

• (vR j /c 2 ) sin (7 q - e 0 + a - c ot) 

R i d o ) 

+ sin (7 - e - ut) l (B-2) 

c - v 2 \ 


2 Vo 

d 0 +R t 


cos (y 0 - e 0 - wf) 


1/2 


&l> 2 0,7o) 

97 0 


vR 2 

— r sm (7„ " e Q + a- r? - cot) 

c -v 


2 2 
C ~ V 


1 + 


2vd Q cos a 


(c/ 2 +J R 2 )(l-(v/c) 2 ) 


Vo , 

Z^ C0S(T “' e » 


CJ/) 


(B-l) 


_ 2 

+ |(l-(v/c) 2 ) [/ + vc? 0 cos a/c 2 

2 

- v/? 2 cos (7 0 - e 0 + a - T? - w/)/c 2 ] 

+ (1 - (v/c) 2 ) 1 [(</ 2 +/? 2 - 2R 2 d Q cos (7 0 - e Q 


where the second line in Eq (B-l) is obtained by using the 
approximation (1 + x) 1 ^ 2 = 1 + (l/2)x with accuracy of the 
order 10“ 10 Substitution of (B-l) into Eq (3) and ignoring 
terms of order 10~ 7 and smaller gets Eq (26) A similar 
argument using Eq (4) gets Eq (27) Consider next 30,/37 o 
and 30 2 /07 o Differentiating Eq (3) and Eq (4) yields 


d0j(C 7 0 ) 



sm (7 0 “ e 0 + 0! " w 0 


-r?-w t)/c 2 -t 2 ]\ ' 

• j(l - (v/c) 2 ) [t + ( vd Q cos a)/c 2 


- vR 2 cos (7 q - e 0 + a - 77 - ojt)/c 2 ] 


• ( yRJc 2 ) sin (7 0 " e 0 + « ~ P " wf) 


+ {(1 - (y/c) 2 ) 2 [t + (vd 0 cos a)/c 2 R d \ 

+ — V sm (7 0 " e 0 - t? - ut) > (B-3) 

2 C - V 2 ) 

- vR j cos (7 - e Q + a - u>t)jc ] ' 


+ (1 - (v/c) 2 ) 1 [(d 2 0 + R 2 - 2 R l d 0 cos ( 7() - e Q 


An approximation similar to that used m (B-l) can then be 
used with (B-2) and (B-3) to establish Eqs (28) and (29) 
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Validation of Roundtrip Charged-Particle Calibrations Derived 
From S- and X-Band Doppler Via DRVID Measurements 

D W Green 

Navigation Systems Section 


In this article it is shown that roundtnp signal path calibrations for charged-particle 
effects can be computed from S- and X-band doppler data The method by which this 
is done is described, and possible error sources are examined It is shown that roundtrip 
calibrations computed from S- and X-band doppler compare to direct measurements 
using DR VID to within the DR VID data noise 


I. Introduction 

The Tracking Systems Analytic Calibration (TSAC) func- 
tion provides calibrations for spacecraft radio metric data 
based on dual frequency (S- and X-band) doppler data These 
calibrations are based on modeling the charged-particle effects 
on the uplink portion of the roundtrip path using the down- 
link only information in the S/X data and certain assumptions 
about the spatial and temporal distribution of the contributing 
charged-particle environment 

The presence of a low-density plasma along the propaga- 
tion path of a radio signal causes the group propagation 
velocity to decrease and the phase velocity to increase The 
DSN ranging system uses measurement of the roundtrip signal 
propagation time and the assumption of a constant value for 
the group velocity to imply the roundtrip distance to the 
spacecraft The presence of charged particles along the path 
causes an apparent increase in this range value The doppler 
measurement system measures the increment in the round- 
trip path length by counting cycles of the returning signal 
and assuming a constant wavelength A positive change in the 
charged-particle content along the signal path will result m a 
decrease in the measured range change from counted (inte- 
grated) doppler The apparent increase in range inferred from 


group measurements is equal in magmtude to the apparent 
decrease inferred from phase measurements taken over the 
same time period The magnitude of this effect is proportional 
to the plasma density and inversely proportional to the square 
of the frequency of the signal involved A comparison of the 
apparent increase in range from differencmg two consecutive 
range measurements with the apparent decrease in range as 
measured by doppler data counted over the same interval 
gives a measurement of the change in the charged particles 
along the signal path An implementation of this techmque 
has been mechanized in the PLanetary Operational (PLOP) 
ranging system and the resultmg data type is referred to as 
Differenced Range Versus Integrated Doppler (DRVID) data 
(Ref 1) The equivalent range change measured by DRVID is 
output in PLOP Range Units (RU) An RU is approximately 
one light nanosecond 

On December 1 , 1978, a test was run at Deep Space Station 
1 4 to determine the effects of movmg the antenna subreflector 
For this pass of data, range, doppler and DRVID data were 
taken from the Voyager 1 spacecraft at both S- and X-band 
DRVID data provide a measurement of the charged-particle 
environment along the entire signal path Thus, simultaneous 
acquisition of dual frequency doppler data and DRVID 
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allows a comparison to be made between a complete round- 
trip calibration reconstructed from the doppler data using a 
model and an actual measurement of the effect by DRVID 
This in turn allows the modeling technique used to be 
validated 

A similar pass of data durmg which both dual frequency 
doppler and DRVID data were taken occurred on April 13, 
1977 for the Viking 2 orbiter 

In the following sections, the technique used to provide 
roundtrip signal calibrations from S/X doppler data are 
described, and several possible error sources are identified 
For the passes of data processed we find that the roundtrip 
calibrations derived from S/X doppler data are equivalent to 
the same data derived from DRVID measurements to within 
the DRVID measurement noise 


P(t ) = plasma contribution computed by re- 
moving ionosphere effect from total 
downlink effect as measured by S/X 
doppler data 

P(t - TPLAS) = plasma effect evaluated at time corre- 
sponding to uplink passage of signal 
through space where plasma is assumed 
to be located 


I(t- RTLT) = ionosphere contribution at transmis- 
sion time 

X = (240/22 1) 2 = correction factor to account for in- 
creased effect on uplink due to lower 
uplink transmitted frequency 


The increment m the S-band cycle count over a particular 
interval is given by 


II. The TSAC S/X Doppler Model 

The model which TSAC uses to compute roundtrip calibra- 
tions from downlink S/X data is pictured m Fig 1 Working 
backward from the reception time, the signal is assumed to 
pass through Earth’s ionosphere at a time negligibly different 
from the reception time The solar plasma is assumed to affect 
the signal as if it were concentrated at the point on the ray 
path closest to the Sun Depending on the Sun-Earth-probe 
geometry, this can occur anywhere from zero to approximately 
500 seconds (light time to Sun) prior to the reception time 
Prior to this, the signal encounters roughly the same plasma 
configuration at a time earlier by a roundtrip light time 
(RTLT) less twice the Earth-to-plasma time Initially, the 
signal is affected by the Earth’s ionosphere at a time negligibly 
different from the transmission time, l e , the reception time 
mmus one RTLT The reliability of the assumption that the 
plasma acts as though it were concentrated at the point on 
the signal path closest to the Sun is examined in Ref 2 

The equation for the total effect due to charged particles 
on a range point received at time t is 


R (r) = I(t) + P(t) + XP(t - TPLAS) + XI (t - RTLT) (1) 


where 

R(t ) = total charged particle effect on range at 
time t 

1(f) = ionosphere contribution computed 
from Faraday rotation data 



( 2 ) 

and the corresponding X-band cycle count is 



( 3 ) 


where 

AS = mcrement in S-band cycle count 
AX = mcrement in X-band cycle count 
t = end of count time 
T = length of count time interval 
fo = S-band downlink carrier frequency 
c = speed of light 

r = time derivative of roundtrip path length to space- 
S/C craft 

r = time derivative of apparent increase in range at 
received S-band frequency due to charged parti- 
cles on uplink leg 

r dp - time derivative of apparent increase in range at 
received S-band frequency due to charged parti- 
cles on downlink leg 
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Note that 


The equation for the computed range correction becomes 


_3_ 

11 


AX- AS 


t 

t-T 



R ( t ) = SX(t) + o(t ) + X [SX(r - TPLAS - tt) - I(t - TPLAS - n) 


+ I(t - RTLT) + a(t - TPLAS TPLAS - n) 


Hence, the accumulative downlink effect during a time mterval + r?(t - RTLT)] 

t = nT is 


sm = E fa 

i~ 1 





AX- AS | 


i T 

(<-i )T 


= A r+K 

dp 


(5) 


Tliat is, the summation of the appropriate linear combination 
of S- and X-band doppler counts yields a measurement of 
the downlink contribution to S-band range due to all charge 
particles to within a constant of integration 


The true range correction being removed, the error in the range 
correction, to first order, is given by 


ER(t) = a(t) + X 


(at - TPLAS) + A(SX- I ) 


f-TPLAS 

r-TPLAS-ir 


+ A77 


f-TPLAS 

f-RTLT 


+ K 


(9) 


Expanding in terms of derivatives 


The contribution to the total downlmk range change due 
to the ionosphere is assumed to be known from Faraday rota- 
tion measurements Thus, the contribution from space plasma 
can be computed as 


ER (t) = a(t) + Xa(t - TPLAS) + n(SX -I) - T}T + K 

( 10 ) 


P(t) = SX(t)-Kt) (6) 

and substituting in Eq (1) above 

R(t) = SX(t) + X[SX(t - TPLAS) - I(t- TPLAS) 

+ I(t - RTLT)] (7) 

III. Error Equation for the Selected Model 

We next consider the following possible sources of error in 
this model 

(1) The S/X doppler measurements may be in error, 
describe the cumulative effect of this error as a(t ) 

(2) The assumed location of the plasma may be in error, 
call the corresponding error in TPLAS, ti 

(3) There may be a measurement or modehng error in the 
Faraday rotation data The cumulative effect of this 
error is 17 (t) 

(4) There is a constant of integration K m the S/X data 


where r = RTLT - TPLAS 

There are two principle uses for the roundtrip calibration 
for charged particles derived from S/X doppler data These 
are 

(1) The calibration of doppler observables 

(2) The calibration of differenced range data 

For both of these phenomena, the quantity of interest is the 
change in the roundtrip measurement over some period of 
time For doppler this period is one count time For differ- 
enced range it is the time between the last range pomt taken 
at one station and the first range point at the next station 
We are interested in the error in the range change over some 
period of time AT Differencing Eq (10) above yields, to 
first order 

A ER = AT \a(t) + \a(t - TPLAS) + n(SX- /) - tt]] 

( 11 ) 

noting that the constant of integration K in the S/X data has 
dropped out 
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IV. Results of Analysis 

A. Results of the December 1, 1978 Test 

Figure 2 shows the total downlink range change for the 
December 1, 1978 pass divided into components due to 
ionosphere and space plasma Figure 3 shows the raw DRVID 
data points at S-band and a reconstructed S-band DRVID 
curve based on S/X doppler data and the model described in 
Section II There is no absolute level information in the S/X 
doppler data, so the absolute level of the derived DRVID 
curve has been adjusted to minimize the RMS difference 
Figure 4 shows the residuals from this adjustment The RMS 
fit error of the overall fit is comparable to the raw DRVID 
noise level as determined by linear detrending of short seg- 
ments of the raw data 


Figure 10 shows the RMS residuals from these fits In 
none of these runs was there any substantial reduction m the 
apparent linear trend in the residuals 

The measurements of ionospheric electron density from 
which the assumed ionosphere contribution is computed are 
derived from data taken by a Faraday rotation polanmeter 
These rotation data are supplied modulo 180 degrees of rota- 
tion As a result, the absolute electron density level may be 
in error by an amount correspondmg to one or more 180° 
rotation mods 

Re-examination of Section III, Eq (11) shows that the 
contribution to range change from an error in the ionosphere 
model is of the order of 


Figures 5 and 6 display the correspondmg results for 
X-band DRVID Again, the RMS error over the entire fit is 
comparable to the data noise level The apparent discontinuity 
in the residuals at about 12 20 can definitely be ascribed to 
subreflector motion 


A R = r) (TPLAS - RTLT) 

The second derivative of the ionosphere error due to setting 
the night time level incorrectly amounts to 0 049 m/hr 2 per 
mod For this case, TPLAS - RTLT is -825 seconds, resulting 
in a contribution of 


B. Initial Results of the April 13, 1977 Test 

Figure 7 shows the total downlink range change for the 
April 13, 1977 pass divided into components due to iono- 
sphere and space plasma Figure 8 shows the raw S-band 
DRVID and the S/X denved DRVID for the April 13, 1977 
Viking 2 pass Figure 9 is a plot of the residuals from this fit 
There is an obvious lmear trend in the entire pass of residuals 
which a lmear least squares fit gives as -3 98 RU/hr In terms 
of the allowable time gap of six hours between differenced 
range points specified in the Voyager SIRD, this amounts to a 
contribution of about 1 75 meters out of a total error budget 
of 6 4 meters The SIRD allocates only 0 5 meters for this 
error 


C. A Detective Story 

What are the possible sources of such an error 2 Some 
likely candidates are 

(1) An error in setting the night time level for the iono- 
sphere 

(2) Errors in the S/X data Possibly frequency offset 
between S-band and X-band synthesizer frequencies 


A R 


0 04 9 m 
hr 2 


• -825 sec 


1 hr 

3600 sec 


per mod • 


7 05 RU 
m 


or 


A R 


RU 

-0 079 -j^- per mod 


Thus an error of -3 98 RU/hr would require about a +50 mod 
ionosphere error, which is at least an order of magnitude too 
large 

Finally, the minimum settable discrepancy between the 
S-band and X-band synthesizer frequencies is 10” 3 * 5 Hz The 
fractional error would then contribute 


10~ 5 Hz 48VCORU 3600 sec , + /240\ 2 
VCOHz ' sec ’ hr 1 \ 221 / 

which is exactly the error we are looking for 


= 3 76 RU/hr 


D. Final Result for the April 13, 1977 Test 


(3) An error in modelmg location of plasma 

In order to test item 3, a series of fits was made, assuming 
locations for the plasma ranging from near Earth to near 

spacecraft 


Figures 11 through 14 show the raw DRVID data and the 
derived DRVID assuming a 1 X 10“ 5 Hz discrepancy between 
the S- and X-band synthesizer frequencies and the residuals 
of the fit between the two at both S- and X-band The resid- 
uals are flat to the noise level of the data Figure 15 shows the 
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RMS fit error as a function of the assumed position of the 
space plasma 


V. Conclusions 

For the December 1, 1978 test, the roundtnp calibrations 
derived from downlink S/X data are equivalent to the mea- 
sured roundtnp effect at both S- and X-band to the level of 
observability m the DRVID data A glance at the computed 
plasma portion of the downlink range change (Fig 2) shows 
that it is within one or two range units of being a straight line 
over the entire pass In terms of the error equations derived 


in Section III, this means that the resulting roundtnp calibra- 
tion on this particular occasion is quite insensitive to the 
assumed position along the ray path of such a plasma 
contribution 

For the April 13, 1977 test, the derived data from S/X 
doppler match the measured DRVID at both S- and X-band 
to what appears to be the noise level of the DRVID data 
There is a well-defined minimum in the RMS fit as a function 
of the assumed plasma location at or near the nominal value 
(see Fig 15) These results together provide a validation of 
the assumptions upon which the S/X calibrations are currently 
based 
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Fig 1 Assumed charged-particle distribution for S/X calibration 

model 



Fig. 2 Ionosphere and space plasma contributions to downlink range change on December 1, 1978 
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Fig 5. Raw and derived X-band DRVID for December 1, 1978 
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Fig 6 X-band DRVID residuals (raw-derived) for December 1, 1978 
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Ionosphere and space plasma contributions to downlink 
range change on April 13, 1977 
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Fig 9 Initial S-band DRVIO residuals (raw-derived) for 
April 13, 1977 
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Fig 10 Initial RMS fit error vs assumed plasma location on 
April 13, 1977 


Fig 8 Initial raw and derived S-band DRVID for April 13, 1977 
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Fig 11 Final raw and derived S-band DRVID for April 13, 1977 
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Recent emphasis on energy conservation has resulted in new efforts to develop more 
efficient, long-life refrigeration techniques for the DSN maser As the new techniques are 
developed, a comparison of their reliability and efficiency to those of the closed-cycle 
helium refrigerator presently used will have to be made Measurements of the refrigera- 
tion capacity have been made on an existing refrigeration system as a function of the 
supply pressure and the return pressure in the Joule-Thomson circuit The performance 
of the refrigerator has been calculated and compared to the ideal Camot-cycle refrigerator 


I. Introduction 

A new effort has been undertaken to study efficient, 
long-life refrigeration techniques for the DSN masers Key 
areas requiring study and development are (1) highly efficient 
helium compression techniques capable of long unattended 
operation, (2) long-life seals and bearings, (3) improved 
system integration to increase the efficiency in the Joule- 
Thomson cooling process or (4) finding alternatives to the 
Joule-Thomson process, such as adiabatic demagnetization, 
to provide liquid helium temperatures at increased efficiencies 

The closed-cycle helium refrigerators (CCRs) currently in 
use in the DSN were developed long before energy conserva- 
tion became a major concern such that reliability was achieved 
at the expense of efficiency As a result over one million hours 
of operation has been logged on some thirty CCRs that have 
been installed in the DSN since 1967 (Ref 1) However, the 
CCR operates between room temperature and liquid helium 
temperatures with only about one percent of Carnot efficiency 
in producing one watt of refrigeration at 4 5 K 


As possible steps to improve upon the efficiency of the 
CCR, studies and development are presently being undertaken 
in each of the three basic subsystems of the CCR the com- 
pressor unit, the expansion engine, which provides refrigera- 
tion down to 15 K, and the Joule-Thomson (J-T) circuit, 
which provides the final refrigeration between 15 and 4 5 K 
A schematic diagram of the CCR is shown in Fig 1 


In an effort to reduce the electrical input power required 
to operate the compressor system, a new concept in the 
mechanical compression of a gas has been under development 
In the new approach, a centrifugal-reciprocating compressor 
concept has been developed which utilizes the centrifugal 
forces produced by the rotating pistons and the housing 
assembly to assist in the compression of the helium gas The 
performance of a single stage feasibility model has been 
demonstrated, and a mathematical analysis of the machine 
has been published (Ref 2) A patent application has been 
filed by NASA for the centrifugal-reciprocating compressor 
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In another area of development the regenerator matrix 
materials used in the 15 K expansion engine of the CCR are 
being studied The 15 K expansion engine utilizes the Gifford- 
McMahon thermodynamic cycle and has proven very reliable 
It is expected that the efficiency and reliability of the expan- 
sion engine may be improved by either improving upon the 
matrix of the presently used regenerator material or by using a 
substitute material Lake Shore Cryotromcs, Inc , under con- 
tract to JPL (Ref 3), has developed a ceramic material which 
is thermodynamically superior (for T< 25 K) to the lead shot 
presently used as a regenerator material Testing of this ceramic 
material in an expansion engine will soon be conducted to 
measure its performance 

Adiabatic demagnetization refrigeration (magnetic cooling) 
is yet another area of study and development The magnetic 
cooling process is being studied as a possible replacement of 
the J-T circuit now used in the CCR to provide the final 
stage of refrigeration Two major advantages of the magnetic 
cooling process over the J-T process would be its high relia- 
bility (since it is a solid state technique) and its high efficiency 
(60-80 percent of Carnot efficiency between 4 5 K and 15 K, 
compared with approximately 20 percent of Carnot efficiency 
for the J-T expansion circuit over the same temperature range) 
After initial testing of the paramagnetic salts, an experimental 
magnetic cooling system will be built and tested 

The three areas of study and development mentioned above 
will be the subject of individual reports to be presented at a 
later time and thus will not be discussed further in this report 
Rather, the scope of this report has been focused on the opera- 
tion of the present CCR to determine its refrigeration capacity 
under various operating conditions It is felt that the study of 
a more efficient refrigeration system should begin with a 
complete understanding of the refrigeration system as it 
presently exists This then provides the baseline for a compara- 
tive study as new materials or methods of operation are inte- 
grated into the system 

The refrigeration capacity for the present CCR was mea- 
sured in two separate tests First, the capacity was determined 
for various J-T supply pressures and various refrigerator 
supply pressures while holding a constant J-T return pressure 
Second, the capacity was again measured, this tune as a 
function of the J-T return pressure while maintaining constant 
supply pressures to both the J-T circuit and the expansion 
engine In the second test particular attention was paid to the 
capacity measurements for return pressures above the critical 
pressure of helium (2 27 X 10 s N/m 2 (33 psia)) to determine 
how well the Joule-Thomson circuit would regulate when no 
liquid was being produced The performance of the CCR under 
these conditions was then calculated and compared to the 
ideal Carnot-cycle refrigerator Before discussing these tests 


and the test results, the method of thermometry used will first 
be described 

II. Temperature Monitors 

Temperature determination at the 4 5 K station was handled 
by two methods resistance thermometry and vapor pressure/ 
gas thermometry A Lake Shore Cryotromcs digital thermom- 
eter (Model DRC-70C) provided direct temperature readings 
using a silicon diode (Model DT-500-Cu) placed at the 4 5 K 
station (see Fig 1) A vapor pressure transducer, located in 
close proximity to the diode, connects to a pressure gauge 
outside the refrigerator and provides the other means of 
temperature sensing The temperature of this station dictates 
the pressure of the gas within the gauge At very low tempera- 
tures, a portion of the gas within the transducer will liquify 
In the small range of temperatures for which this liquid and 
its vapor will exist in equilibrium, the pressure gauge readings 
may be compared to a known temperature-vapor pressure 
curve to determine the temperature at the cold station Above 
this temperature range the device behaves as a gas thermom- 
eter, with the pressure readings being proportional to the 
temperature 

The relationship between the two temperature sensing 
devices at various temperatures was measured by adjusting the 
J-T return pressure and by applying resistive heating to the 
4 5 K station The results are shown in Fig 2 with the resist- 
ance thermometer readings plotted on the y-axis and vapor 
pressure readings plotted along the x-axis Also shown is the 
curve for the accepted values for the temperature-vapor 
pressure relationship as tabulated in the American Institute of 
Physics Handbook (Ref 4) The 0 1 K difference between the 
curves at the low pressures may be attributable to the calibra- 
tion error of the digital thermometer The measured data were 
not corrected for this error The diode has been temperature 
cycled several times and has shown good repeatability in the 
data The 0 1 K sensitivity limit of the digital readout adds 
some additional uncertainty to the data at the low pressures 
The vapor pressure device was unable to maintain a liquid/ 
vapor equilibrium to the 2 27 X 10 s N/m 2 (5 2 K) limit as 
shown by the deviation of the shape of the measured curve 
from the AIP curve near 1 86 X 10 s N/m 2 (27 psia) Thus, the 
vapor pressure thermometer should not be relied upon for 
accurate temperature measurements above 1 86 X 10 s N/m 2 
(4 9 K) where the device changes to a gas thermometer 


III. Refrigeration Capacity Measurements 

For a fixed pressure operation, the refrigeration capacity 
at the 4 5 K station of the CCR is determined primarily by 
the flow rate in the J-T circuit, by the temperature of the 
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15 K station, and by the efficiency of the heat exchanger 
between the 15 K station and the J-T valve This capacity is 
equal to the heat of vaporization of the liquid helium pro- 
duced by the J-T expansion plus the heat capacity of the 
cooled gas that had not liquified by the J-T expansion Elec- 
trical input power levels below the cooling capacity will 
vaporize a fraction of the liquid reserve before reaching a 
new hquid/vapor equilibrium, power levels in excess of this 
cooling capacity will deplete the liquid reserve and warm the 
refrigerator It has also been observed that above 5 2 K, the 
critical temperature at which liquid helium can be maintained, 
the cooled gas from the expansion process also has a signifi- 
cant cooling capacity This will be discussed shortly 

Determination of the refrigeration capacity of the CCR was 
made by applying an electrical load to the resistive heating 
element located on the 4 5 K station An electrical power load 
above the assumed capacity was applied to quickly deplete the 
reserve liquid in the system, causing the temperature and 
vapor pressure at this station to rise By incrementally decreas- 
ing the input power, the temperature and the vapor pressure 
were kept from rising above a nominal 4 8 K and 1 72 X 10 s 
N/m 2 (25 psia), respectively (well below the critical tempera- 
ture and pressure of 5 2 K and 2 27 X 10 s N/m 2 (33 psia), 
respectively), so that the J-T expansion process would con- 
tinue to produce liquid When the power was decreased to 
below the cooling capacity, the temperature and the vapor 
pressure dropped, indicating liquid was being stored This 
provided a quick method of approximating the capacity A 
more accurate capacity determination could then be made by 
applying power levels about this point for periods as long as 
one hour to determine the highest input power for which the 
temperature and the vapor pressure would stabilize In this 
manner the refrigeration capacity has been measured for 
three J-T supply pressures and two different refrigerator 
supply pressures while holding the J-T return pressure con- 
stant The results are tabulated in Table 1, where the refrigera- 
tion capacity has been denoted as <2 meas and the temperature 
of 4 5 K station as T c (The remainder of the data m the table 
will be discussed later ) The capacity measurements were made 
throughout the day, during which the ambient temperature 
surrounding the compressors would change The effect of the 
daily temperature cycling has been observed to produce up to 
a 50-mW capacity change on previous CCRs However, tem- 
perature changes of 7°C during the measurements on this 
CCR showed no noticeable change in the refrigeration capacity 

The temperature variation as a function of the J-T return 
pressure and the various heat loads was also monitored The 
results have been tabulated in Table 2 and are shown in 
Figs 3 and 4 In Fig 3 the temperature at the 4 5 K station 
has been plotted as a function of the applied heat loads The 
values for the constant pressure curves are given for the no 


load condition since at the higher heater power levels the 
return pressure will decrease slightly Each pressure curve ends 
at the highest heater power for which the temperatures and 
pressures would stabilize, increasing the heater power beyond 
this pomt would cause the refrigerator to warm excessively 
(The power levels could be incremented in approximately 
20-mW intervals ) The maximum power level at which the 
stabilization occurs may be considered the cooling capacity 
for that particular J-T return pressure (or for that particular 
temperature) 

Figure 3 shows that there is a nearly linear degradation in 
cooling capacity as the temperature is increased from 4 5 to 
8 K, and that there is still a significant cooling capacity for 
temperatures above 5 2 K, the critical temperature above 
which liquid helium cannot exist Above this critical tempera- 
ture the refrigeration capacity must be obtained non- 
lsothermally from the sensible heat of the gas rather than from 
the isothermal heat of vaporization Although temperature 
and pressure regulation at the 4 5 K station should be diffi- 
cult to maintain because of the low thermal capacity of the 
helium gas, the data show very little temperature change with 
change in heater power for temperatures below 6 K 

While the present CCR relies on the liquid helium reserve 
to insure additional temperature stability, the operation of the 
CCR at increased J-T return pressures (so that the cold station 
temperature is above 5 2 K) also provides reasonable tempera- 
ture stability This depends on the provision that the heat 
load of the device being cooled is sufficiently less than the 
refrigeration capacity of the CCR under these operating 
conditions 

For the CCR tested, operating the J-T return pressure as 
high as 3 10 X 10 s N/m 2 (3 atm) increased the temperature 
by only 1 1 K over the nominal 4 5 K achieved for the normal 
operating return pressure of just over an atmosphere This 
result is significant as there are numerous applications (such 
as cooling superconducting magnets or Josephson devices) 
that could be satisfied with an efficient and reliable 5 5 K 
refrigerator Another application could be the cooling of a 
maser since the refrigeration capacity of the CCR at 5 5 K is 
much greater than the heatload produced by an operating 
maser Although it may be desirable to optimize maser perfor- 
mance by operating the maser as cold as possible (Ref 5), 
there are other factors, such as energy requirements and 
reliability, which determine the necessary trade-offs to estab- 
lish the temperature that provides the best overall performance 
The higher return pressure would permit a reduction in the 
electrical power requirements by operatmg with smaller com- 
pressors, which may also have the added advantages of pro- 
longed maintenance-free operating times and an increased 
overall lifetime of the system 
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Figure 4 shows the temperature at the 4 5 K station as a 
function of the J-T return pressure Portions of several con- 
stant power curves have also been shown The vapor pressure 
curve has also been placed on the graph to show the similar 
shape between the vapor pressure curve and the J-T return 
pressure when no heater power is applied The return pressure 
curve appears to be a direct extension of the vapor pressure 
curve showing no discontinuity or change in slope as it extends 
beyond 2 27 X 10 s N/m 2 (33 psia) The pressure differences 
between the two curves occur because the J-T return pressure 
gauge monitors the pressure at the room temperature gas line 
rather than at the 4 5 K station Restriction in the heat 
exchanger causes the pressure at the external pressure gauge 
to be lower than indicated by the vapor pressure gauge (It 
should be noted that the pressure drop due to this restriction 
has a small effect on the available refrigeration capacity and 
directly affects the temperature at which the refrigeration 
takes place (Ref 6) ) The appropriate test points for pressure 
and temperature measurements are shown in Fig 1 


compressor during one cycle (1-2-3-4-1) is 


W 

-f - T H (S, - ^ - V + - *3> - <V V 

where 

x = (h 1 - h 2 )/(h i - h f ) 

is the fraction of the gas liquified and 

h , = xh,+ (\- x)h 
3 / g 

If no liquid is produced, the cycle follows path (l-2-3'-4'-l), 
and the net work required is 


W 


m 


= W 5 2 )-(V a 2> 


IV. Refrigerator Performance Analysis 

When a heat load corresponding to the cooling capacity of 
the machine is applied at the 4 5 K station, the temperatures 
of the incoming and outgoing streams of gas just above the 
third exchanger are equal Were they not, the system would 
either cool down or warm up To determine the performance 
of the J-T loop below the 15 K station, the remainder of the 
refrigerator may be thought of as a single compressor operat- 
ing at 15 K This refrigerator may then be analyzed as if it 
were a Linde-Hampson refrigerator as shown in Fig 5 The 
Linde-Hampson refrigerator has been discussed thoroughly 
elsewhere (e g , see Refs 7, 8, and 9) so only a brief review of 
the analysis will be presented here 

It will be assumed that the final heat exchanger is 100 per- 
cent effective, that there are no heat leaks into or out of the 
system, and that the J-T expansion represents the only irre- 
versible pressure drop throughout the cycle shown in Fig 5 
Under these conditions the heat absorbed from the source 
during one complete cycle may be expressed as 



that is, the net refrigeration per unit flow rate depends solely 
on the enthalpies of the incoming and outgoing gases, and is 
independent of what happens below the top of the heat 
exchanger If liquid is produced, the net work required by the 


The performance of the refrigerator is expressed as the coeffi- 
cient of performance by 

Q nf 

COP = — 

w , 

net 

The efficiency of this refrigerator is determined by using the 
figure of merit of the refrigerator, defined by 

COP QJW net 

FOM = = — - — 

COP c t c kt h -t c ) 

where COP c is the coefficient of performance for a reversible 
Carnot cycle, and the FOM varies from zero to unity The 
coefficient of performance and the figure of merit for the 
measured refrigeration capacity may be evaluated in the same 
manner 

The results for the analysis have been included in Tables 1 
and 2 The values for the entropies and the enthalpies used 
to determine W net and Q ref were estimated from the 
temperature-entropy diagram given in Ref 10 As the only 
other non-measured parameter, T H was assumed to remain 
constant at 15 K throughout the measurements since there 
was no provision made to measure the temperature of the 
15 K station The 15 K temperature is representative of the 
temperature of the 1 5 K station of the CCR when cooling a 
maser amplifier Heat loads in excess of this amount would 
raise the temperature of this station by an unknown additional 
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amount, resulting in a smaller difference between Q mea s an d 
Q ref (which depends on this temperature) 

The results shown in Table 1 indicate that operating at the 
higher refrigerator supply pressure increased the refrigeration 
capacity at the 4 5 K station The increased supply pressure 
results in a lower temperature at the 15 K station In turn, 
the temperature of the J-T supply gas, which is in a recupera- 
tive heat exchange with the 15 K station, is lowered Expand- 
ing the gas at the lower temperature results in the increased 
cooling capacity The data presented in Table 2 show that 
the capacity increases with decreasing return pressure This is 
also consistent with the temperature-entropy diagram How- 
ever, to insure contaminants will not be drawn into the gas 
lines, a positive pressure is maintained for the J-T return gas 

V. Summary 

Several important facts have been established with the 
testing of this CCR, which may be considered representative 


of all DSN CCRs First, the vapor pressure gauge worked very 
well at the extreme low temperatures, but should not be 
relied upon beyond 1 86 X 10 5 N/m 2 (27 psia) (4 9 K), 
where it prematurely begins to function as a gas thermometer 
In contrast, the resistance thermometer showed a very small 
but constant error across the 4 through 8 K range, which could 
be corrected to give accurate temperature readings Second, 
the maximum cooling capacity of this particular CCR was 
determined to be 1 00 W when operating with a refrigerator 
supply pressure of 2 20 X 10 6 7 N/m 2 (305 psig), a J-T supply 
pressure and return pressure of 1 82 X 10 6 N/m 2 (250 psig) 
and 111 X 10 s N/m 2 (16 1 psia), respectively, and a com- 
pressor temperature of 25°C Finally, the refrigeration cap- 
acity and the temperature stability for J-T regulation above 
the critical temperature and pressure of 5 2 K and 2 27 X 10 s 
N/m 2 (33 psia) remained quite high for temperatures and 
return pressures as high as 6 K and 3 44 X 10 s N/m 2 (50 psia), 
respectively This cooling capacity at the increased return 
pressures is significant and may permit the utilization of 
smaller compressors with reduced power requirements 
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Table 1 Refrigeration capacity as a function of refngerator 
supply pressure and J-T supply pressure 


Refrigerator Supply, N/m 2 X 10 5 (psig) 

22 04 
(305) 

22 04 
(305) 

22 04 
(305) 

20 31 
(280) 

20 31 
(280) 

Refrigerator Return, N/m 2 X 10 S (psig) 

8 60 
(110) 

8 60 
(110) 

8 60 
(110) 

7 91 
(100) 

7 91 
(100) 

J-T Supply, N/m 2 X 10 s (psig) 

2031 

(280) 

18 25 
(250) 

16 18 
(220) 

18 25 
(250) 

16 18 
(220) 

J-T Return, N/m 2 X 10 s (psia) 

1 11 
(16 1) 

1 11 
(16 1) 

1 11 
(16 1) 

1 11 
(16 1) 

1 11 
(16 1) 

J-T Flow, SCFM 

1 48 

1 39 

1 25 

1 38 

1 28 

r c ,K 

4 7 

4 6 

4 5 

4 5 

4 5 

COP c 

0 456 

0 442 

0 429 

0 429 

0 429 

"net.* 

10 82 

9 87 

8 61 

9 79 

8 85 

<W w ( i0 02 W > 

1 037 

1 000 

0 834 

0 896 

0 795 

^^meas 

0 0958 

0 101 

0 0969 

00915 

0 0898 

«™meas 

0 210 

0 229 

0 226 

0 213 

0 209 

2 r ef 

1 81 

1 60 

1 31 

1 59 

1 35 

COP le{ 

0 168 

0 162 

0 152 

0 162 

0 153 

™ V ref 

0 368 

0 367 

0 354 

0 378 

0 357 


T h = 15 K 

CCR tested SN 78005 
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Table 2 Refrigeration capacity as a function of 
J-T return pressure 


P , N/m 2 x 10 s (atm) 

1 11 
(1 1) 

1 82 
(1 8) 

2 33 
(2 3) 

2 74 
(2 7) 

3 24 
(3 2) 

3 95 
(3 9) 

4 96 
(4 9) 

5 98 
(5 9) 

J-T Flow, SCFM 

1 39 

1 22 

1 20 

1 17 

1 13 

1 05 

0 99 

0 90 

T C K 

4 6 

5 2 

5 7 

6 0 

65 

7 3 

8 1 

8 6 

COP 

c 

0 442 

0 531 

0613 

0 667 

0 765 

0 948 

1 17 

1 34 

w ,,w 

net 

9 87 

7 11 

6 32 

5 45 

4 81 

4 15 

3 34 

2 58 

Q , W (±0 02 W) 

^meas 

1 000 

0 896 

0 803 

0 703 

0616 

0 489 

0 396 

0 294 

COP 

meas 

0 101 

0 126 

0 127 

0 129 

0 128 

0 118 

0 119 

0 114 

FOM 

meas 

0 229 

0 237 

0 207 

0 193 

0 167 

0 124 

0 101 

0 0851 

e re r w 

1 60 

1 31 

1 21 

1 13 

1 05 

0 884 

0 726 

0 606 

CW ref 

0 162 

0 184 

0 191 

0 207 

0 218 

0 213 

0 217 

0 235 

FOM , 
ref 

0 367 

0 347 

0 312 

0 310 

0 285 

0 225 

0 185 

0 175 


Refrigerator supply 22 04 X 10 5 N/m 2 (305 psig) J-T Supply 18 25 X 10 s N/m 2 (250 psig) 

Refrigerator return 8 60 X 10 S N/m 2 (110 psig) T ^ = 15 K 

P 2 = 18 25 X 10 5 N/m 2 (18 atm) CCR tested SN 78005 
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Fig 1 Schematic diagram for CCR 
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Fig 3 Temperature as a function of the applied heat load for 
several J-T return pressures 
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Fig 4 Measured cooling capacity as a function of 
J-T return pressure 




Fig 5 Thermodynamic cycle and flow diagram 
for Linde-Hampson refrigerator 
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The International Comet Exploration spacecraft will spend 2 4 years in a solar electric 
propulsion stage (SEPS) cruise mode, during which the health of the spacecraft will be 
monitored at least once a day Although a daily two-hour DSN pass is now planned, the 
same surveillance function can be carried out by a much smaller, cheaper system The 
system studied here uses a commercial fast Fourier transform spectrum analyzer for 
noncoherent, synchronous detection of an alternating-tone signal received by a small 
parabolic antenna (0 8 to 3 mm diameter) Integration times vary from 3 7 to 9 5 min 


I. Proposed Surveillance Method 

The proposed International Comet Exploration (ICE) 
spacecraft will be capable of monitoring its own health during 
the solar electric propulsion stage (SEPS) cruise mode It 
would transmit one of several signals (five, perhaps) whose 
meanings range from “I’m OK” to “Help 1 ” The signal would 
be detected by a small, dedicated ground receiving system If 
the spacecraft says that it has a problem or if no signal is 
detected, then a regular DSN tracking session would be 
scheduled 

To reduce size, weight, and pointing accuracy requirements 
we shall try to make the ground antenna as small as possible 
The data rate can be as low as one word (two or three bits) per 
hour or per day Under these circumstances we choose non- 
coherent detection of a multiple frequency shift keyed 
(MFSK) signal by a spectrum analyzer The main problem is 
acquisition of the signal, which is corrupted by oscillator phase 
jitter, displaced by Doppler error, and smeared by Doppler 
rate error Consequently there is an upper bound on integra- 
tion time, which in turn puts a lower bound on signal power, 
no matter how much time is nominally available for detection 


Section II and Fig 1 describe the overall surveillance sys- 
tem and set down the assumptions used in the analysis Sec- 
tion III and Table 1 give four numerical examples The rest of 
the article (up to Section VIII) consists of back-up material 
Sections IV and V derive certain mission-dependent param- 
eters (signal power and Doppler rate error) needed by the 
analysis, which is carried out in Sections VI and VII Conclu- 
sions are offered in Section VIII 


II. Design Assumptions 

A conceptual sketch of a surveillance system is given in 
Fig 1 We list and discuss our assumptions 

(1) The RF signal consists of alternating sine-wave tones 
with frequencies f 0 and f l , each held for the word time 
T d (but see item (8) below) Such a signal is proposed 
in Refs 1 and 2 for MFSK acquisition The difference 
/i - f 0 , having one of several possible values, consti- 
tutes the message word These values should be much 
less than a priori Doppler uncertainty, discussed below 
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(2) Perfect word synchronization is assumed The word 
time T d is long enough (several minutes), the range is 
known well enough, and spacecraft time is accurate 
enough for the ground system to know when the transi- 
tions between / 0 and f x occur Word sync can be 
refreshed whenever the regular DSN tracks the space- 
craft 

(3) The short-term frequency stability of the spacecraft 
oscillator that generates the signal is specified by a plot 
similar to Fig 2, which plots the square root of Allan 
variance, A/// 0 , vs integration time r The behavior 
shown is typical of quartz crystal oscillators At r x 
there is a transition between the behaviors A f/f 0 = 
const/r and A f/f 0 = o y = const To avoid signal power 
degradations above 1 dB we require 

/oVi< 013 (0 

which comes from Eq (13) below Typically, Tj = 0 5 
sec, in which case 


o < 10 4 5 * * * * 10 S-band 

.v 

a <3X10“ 11 X-band 

y 


(4) The transmitter power is 5 W for S-band, 10 W for 
X-band All the radiated power goes into the surveil- 
lance signal At present, X-band is actually not suitable 
for this system (which is supposed to reduce the need 
for DSN tracking) because X-band downlink requires 
an uplink for pointing the spacecraft antenna Never- 
theless, we have provided X-band designs in case a 
stand-alone X-band downlink is developed 

(5) The ground antenna is a circular paraboloidal dish This 

assumption is intended only as a starting point for 
other designs, a circular beam shape may not be best 
for this application We shall design both tracking and 
nontracking systems In a tracking system the antenna 

stays pointed at the spacecraft for the length of time 
required by the detection process In a nontracking 

system the antenna is pointed at the correct decimation 
on the meridian and the spacecraft sweeps through the 

beam Thus detection time is limited by the beam- 

width The antenna has to be tilted slightly from day to 
day or week to week 

A GaAsFET receiver is used See Table 2 for noise 
temperatures 


(6) Doppler and Doppler rate uncertainty are critical to a 
MFSK system The signal must hold still long enough 
for it to be located Let v be the component of vehicle 
acceleration provided by the SEPS We assume a Dop- 
pler rate error f given by 

f'T (2) 

where X is the wavelength This is an upper bound for 
the error that would occur if the SEPS were to quit 
entirely or even start blasting in the wrong direction 
Our Doppler uncertainty bandwidth W is computed by 
letting /accumulate for one day 

W = 86400/ (3) 

We compute a maximum SEPS acceleration of 1 8 X 
10 -4 m/sec 2 during the SEPS cruise II, 1986 0 to 
1988 4 (Section V below) This gives 

/ = 0 0028 Hz/sec, W = 240 Hz at S-band (4) 

/ = 0 010 Hz/sec, W = 870 Hz at X-band (5) 

It is assumed that the prediction errors of all other 
contributions to Doppler and Doppler rate are less than 
these numbers To be conservative we shall make the 
bandwidth of the detection system at least W Actually, 
the system might be able to get by with less bandwidth 
because a large frequency deviation is itself cause for 
alarm Therefore, a DSN tracking session would be 
scheduled whether or not the surveillance system finds 
the signal 

(7) A commercial FFT spectrum analyzer with 256 chan- 
nels and a band analysis mode is assumed No time 
window shaping is used, we found that the uniform 
window performs better than a hanmng window In 
the system of Fig 1 we use a heterodyne frequency f h 
outside the nominal uncertainty band to avoid foldover 
of the spectrum 

(8) The detection strategy uses “L-look” spectral analysis 
and threshold detection (Ref 3) as follows The band- 
width b of an analyzer channel must satisfy 256 b > W 
Each data segment of duration l/b results in a spec- 
trum, the system sums L of these spectra and declares a 
signal present if the summed output of some channel 
exceeds a certain threshold above the average channel 
noise power Thus T d = L/b Once / 0 or /, is found, 
the analyzer bandwidth can be reduced for faster detec- 
tion of the other frequency In fact, the spacecraft 
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might simply spend less time on /, than on / 0 Acquisi- 
tion would be made on / 0 (perfect word sync as- 
sumed), then f l would be detected with a narrower 
bandwidth Other schemes suggest themselves, here we 
shall examine only the acquisition problem 

Other assumptions about the mission appear in Sections IV 
and V 

III. Sample Design Parameters 

Table 1 gives the numerical parameters that define four 
possible surveillance systems For both S-band and X-band we 
have designed a tracking system and a nontracking system The 
tracking systems, in which the antenna follows the spacecraft 
for as long as the detection system needs to see the signal, 
represent the lowest antenna gams consistent with a detection 
error probability of 10 -4 This minimum gain is approxi- 
mately proportional to the square root of the Doppler rate 
error / (Eq (29) below) Antenna diameter is 2 3 m for 
S-band, 0 8 m for X-band The nontracking systems, in which 
the spacecraft sweeps across the stationary antenna beam, use 
slightly larger antennas (3 m for S-band, 12m for X-band) 
because the reduced detection time requires more signal 
power While the nontracking systems are mechanically less 
complex than the tracking systems, they are also less flexible 
and reliable because they have only one brief shot per day at 
the spacecraft 

No margins have been built into these designs except for 
those created by our approximations (positive, we hope) 
Actual systems might have to use antennas with higher gains 
and narrower beams 

IV. Power-to-Noise Ratios 

Using a design control table format (Table 2), we derive the 
ratio P/N 0 for an ideal isotropic ground antenna, where 

P = carrier power 

N 0 = one-sided noise spectral density 

Maximum range, 5 8 X 10 8 km (3 9 AU), is assumed 

The results are 

(-£-] = 5 4X 1(T 4 sec -1 S-band 


= 6 0 X 10 4 sec 1 X-band 


V. SEPS Acceleration 

Standard rocket equations (Ref 4) give 

2P_ 

V Mgl 

where 

v = vehicle acceleration due to rocket propulsion 
P = power expended on propellant 
M = vehicle mass 

/ = specific impulse = V/g, V = exhaust speed 
g = Earth gravity, 9 8 m/sec 2 

Near the end of the long SEPS cruise, just before the Tempel 2 
rendezvous, 3 thrusters will be used with 10 kW available to 
them We estimate P = 5630 W, M = 2200 kg, /= 2940 sec 
Therefore, 

v = 1.8 X 10 -4 m/sec 2 

At the time of greatest range, 1987 7, we estimate 
v = 6 3 X 10 -5 m/sec 2 

To be conservative we shall design for a combination of 
maximum acceleration and maximum range 


VI. L- Look Spectral Analysis 

This section is about threshold detection by a spectrum 
analyzer of a corrupted sine wave in white Gaussian noise For 
a pure sine wave, see Ref 3 Here, we attempt to account for 
the degradations due to (1) position and drift of the signal 
relative to nearby analyzer channels, and (2) short-term oscilla- 
tor phase jitter Our analysis is approximate and incomplete, 
representing the best we have been able to achieve in the time 
available A rigorous treatment would be a research project 
Throughout, we neglect the effects of sampling 

A. Analyzer Channel Outputs 

Let X be the (suitably normalized) sum of L successive 
outputs of a particular analyzer channel of width b Then X is 
approximately a scaled version of a noncentral chi-squared 
random variable with 2 L degrees of freedom (Ref 3) Its mean 
and standard deviation are approximately 
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E(X) = N 0 b+P' = N Q b(l +/■') (6) 

N n b , 

aW =^f- (7) 

where r' = P'/(N 0 b ) and P' is the average signal power in the 
channel during the L looks If the channel contains only 
noise then r - 0, else 


i <8) 

' " »J> ) 

where G pd is the loss (<1) from position and drift of the 
signal during the L looks, and G j is the loss from short-term 
jitter of the spacecraft oscillator Because these two losses are 
actually intertwined, they deserve a unified treatment For the 
present, we simply compute them separately and multiply 
them 

B. Position and Drift Loss G ^ 

This loss depends on the Doppler rate error / In time Lib 
the signal frequency drifts by Lf/b, in units of the channel 
width b the drift width is 


mean loss” between 0 3 dB and 1 5 dB According to Ref 5, it 
is difficult to include this loss in performance evaluations 
because little is known about the exact characteristics of phase 
jitter We shall nevertheless use the information in the Allan 
vanance plot (Fig 2) to obtain a theoretical estimate of the 
jitter loss G. The result will put an upper bound on the 
long-term Allan variance a 2 and a lower bound on the channel 
bandwidth b 

If the oscillator were modulated by either a stationary 
phase process 0(f) or a stationary frequency process 0(f) the 
situation would be in hand because the oscillator output 

x(f) = cos [2 rr/ 0 f + 0(f) + 0 o ] 

would be stationary (We assume 0 O is uniformly distributed 
and independent of 0(f)) Assuming that 0(f) is Gaussian, 
Middleton (Ref 6) and others compute the spectrum of x(t) 
given the spectrum of 0(f) or 0(f) 

Unfortunately, the typical Allan vanance behavior A f/f 0 = 
a y = const for large r indicates the presence of nonstationary 
“flicker of frequency” modulation Previous theory does not 
apply Since the behavior A f/f Q = const/r for small t indicates 
that a stationary phase process is also present, we adopt the 
following model for 0(f) 


5 = ^7 (9) 

b 2 

Let the frequency ongin be centered in the channel Let mb be 
the signal frequency at the midpoint of the L-look time 
interval If m = 0 the midpoint of the frequency drift is at the 
center of the channel, if m= 1/2 the dnft midpoint is at an 
edge of the channel We call these cases the center case and the 
edge case If the drift during a single look is much smaller than 
b, then G pd is approximately the average of sine 2 (itf/b) over 
the drift interval, thus 


0(f) = 0(f) + 0(f) 

where 0(f) and 0(f) are independent Gaussian processes, 0(f) is 
stationary, and 0(f) has stationary second differences The 
relevant second moment properties of 0(f) and 0(f) are charac- 
terized from Fig 2 by 

° 2 (<P) = ^/ov.) 2 00 


pd 


G pd 5 ) = J 


fm+6 /2 / \ 

/ “ 

Jm-b/2 ' ' 


du (10) 


W 27r/ ) = 


(In 4 )/ 3 


(12) 


C. Jitter Loss G ; 

Using a spectrum analyzer, first observe a perfect sine wave 
whose frequency is centered in one of the channels Then 
replace the sine wave by a real oscillator having the same 
average power Because of short-term phase jitter, the height 
of the spectral peak will decrease According to Ref 2, mea- 
surements of S-band transmitters have placed this “spectral 


where S ^ is the one-sided spectral density of 0(f) The l// 3 
spectrum has a rigorous meaning in the theory of processes 
with stationary nth differences Eq (11) does not depend on 
the shape of the spectrum of 0(f) as long as 1/r, is much 
smaller than the bandwidth of 0(f) One obtains Eqs (1 1) arid 
(12) from conventional stationary -process formalism by inte- 
grating the Allan-vanance filter response times the spectral 
density of 0(f) 
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In Appendix A we denve the following approximation for 
the jitter loss 


Since N Q b can be estimated as the average channel power, Eq 
(15) determines the threshold 


In 



= (2 nf Q a y ) 


1 


(36 In 4 )b 2 


(13) 


Suppose that a signal is present The probability that we 
miss it or that a spurious channel output exceeds the threshold 
is approximately 


Although the analysis is not rigorous, the author believes that 
Eq (13) shows the main dependencies of the loss on the 
Allan-vanance turnover point (Tj, o y ) and the channel width 
b The first term on the right comes from the stationary noise 
<p(t) and is already well understood The second term comes 
from the flicker-noise jitter i p(t), as b gets smaller, the integra- 
tion time gets longer, and the effect of the low-frequency 
phase variations predominates 

Within the confines of the analysis, Eq (13) actually gives 
an upper bound for the loss We suggest using it as long as it 
indicates a loss below about 2 dB 

0. Detection Criteria 

We set a detection threshold N 0 br, u where 1 < t? < r' A 
signal is detected if the outputs of a set of adjacent channels 
are more than N Q bri, and all other channel outputs are less 
than N Q br\ 

Fix an error probability e Anticipating the discussion 
below, we define z Q ,z l , and z 2 by 

0( Zo )=^G(z 1 ) = e,Q(z 2 )=VT (14) 

where n is the number of channels and Q(z) is the probability 
that a standard Gaussian is more than z 

Suppose that no channel has a signal The probability of 
false detection is approximately 


II P{X<N 0 bn\r' =r' k }+P FD (17) 

k 

where r' k is the effective SNR in the Arth channel We require 
the first term of Eq (17) to be at most e This term depends 
on the drift center m and the drift width 6 = Z,//h 2 For a 
given /, 6 is at our disposal but m is not Thus our detection 
strategy has to work for all m Without loss of generality we 
can assume 0 < m < 1/2 Lack of time forces us to consider 
only the center case m = 0 and the edge case m = 1/2 Also, we 
assume 5 < 2 so that at most two channels have significant 
amounts of signal 

For the center case we consider only one channel in Eq 
(17) Setting 

P{X<N Q b V }=e 
We obtain from Appendix B that 

x V 1 + 2 r 

r? = l+r' — (18) 

s/r 

where 


x, = 2] v i - M4Z.) -jc^ 2 - i)vwm 09) 


P FD = nP{X>N 0 b n \r' = 0} 
Setting P FD = e we obtain from Appendix B that 


0 = 


1 +2 r 
(l+r ') 2 


( 20 ) 



and r' is calculated using G pd { 0, 6) From Eqs (15) and (18) 
(15) we get the detection criterion 


where 


* 0 +*j y/T+2r' 

vT 


( 21 ) 


x 0 = z 0 V 1 - l/(4L) + ^(z 0 2 - l)V r T7(4l) (16) 


for the center case 
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In the edge case the two channels on either side of the drift 
midpoint m = 1/2 see an equal amount of signal Therefore we 
set 

[P{X<N 0 b V }] 2 = e 
which leads as before to 

x V 1 + 2 r' 

t? — (22) 

sfL 

where 


In the two examples below it turns out that the point (S opt , 
r m in) 1S ^e intersection of the r = /^(S) and r = r 2 ( 5) curves 
(Fig 3) For S < 5 opt> r 2 ( 8) is larger, for 5 > S opt , rj(S) is 
larger For both examples, 5 opt = 16 approximately 

We can explain this from another point of view For a given 
r there is a range of L that will support detection with a given 
error probability If L is too small there is not enough time to 
build up the signal over the noise If L is too large the drift 
spreads the signal over several channels The minimum r is the 
one that shrinks this /--range to a single point 

Suppose now that r > r min Assume that for 5 < 5 opt the 
edge case applies r 2 ( 6) > t-j( 8) The required L can be found 
by iterating 


x 2 = z 2 V1-/3/(4Z,)-|(z 2 2 - Y)y/W L) (23) 

and r' is calculated using G pd (1/2, 6) The detection criterion 
for the edge case is 


L = 


x Q +x 2 V 1 +2 r' 2 

/ 

r 


whose right side depends slowly on L 


(28) 


, x Q +x 2 y/l+2 r' 
r = 

VL 


E. Approximations 

Put e = 10- 4 , 8 opt = 1 6, G pd = G pd (l/2, 0) = 0 41 , and 
compute G } from Eq (13) Then 


Let e, /, and b be fixed Given 5 (or L ) the right sides of 
Eqs (21) and (24) depend slowly on r', so they can be iterated 
to get the required r for the center and edge cases Computing 
G pd (m, 5) (m = 0 and 1/2) and G ] we get the channel SNR 


P 

N 0 b 


G H G 

pd / 


as two functions of 5, say r,(5) (center case) and r 2 ( 5) (edge 
case) Since we don’t know which case holds we must choose 
the worst one by setting 


r(8) = max [/■ 1 (6),r 2 (5)] (25) 

Finally, we get the optimum values of 8 and L, and the 
minimum required r, by minimizing Eq (25) 


r = r(8 ) = minr(S) 

min v opt 7 v 7 

6 <2 


(26) 


8 ,b L 

r _ Opt 

"° pt f 


(27) 


W,) n 


'mm* 


G .G sfd 

pd i v 


-y/f (29) 


opt 



for P/N q > (P/N q ) i 


8 


2 b 


. G „G , 
L Pd i J 



(30) 


VII. Tracking and Nontracking Systems 
A. Design Equations 

From Ref 7 we extract some nominal relationships for 
parabolic antennas Let d = diameter, X = wavelength, 6 = 
angle off-axis Then 

On-axis gain G a = (777-/4) (</ 2 /X 2 ) 

Half-power beamwidth Q a = 63 X/d deg 

Off-axis gain G(d ) = G fl exp (-In 16 d 2 /df) 

We get the weakest possible tracking system by specifying 
G a such that 
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where r IS0 = (P/N 0 ) lso /b The only a pnon restriction on the 
word time T d = L opt /b = 5 opt b 2 lf is the tracking time 
available 

If the antenna is fixed, the spacecraft simply passing 
through the center of the beam, then word time is limited We 
arbitrarily designate the available antenna time by 


B. Derivation of the Examples 

1 S-band Recall that / = 0 0028 Hz/sec, (/ > /A r 0 ), so = 
5 37 X 10~ 4 sec -1 From Eq (4) we choose an analyzer 
bandwidth of 250 Hz (the next bandwidth above 240 Hz for 
the analyzer we have in mind) Since we assume n = 256 
channels, we have b = 1 Hz Equation (13) with a y = 10~ 10 , 
r j =05 sec gives G / = 0 8 1 


T = 240 6 sec 

a a 


the time it takes a point on the celestial equator to traverse an 
angle dj 2 If the r mm system happens to have T d > T a then 
we must increase G a until T d < T a This can always be done 
because T d = Lib is roughly proportional to 1/r 2 , whereas T a 
is proportional to 1/x/T We also have to account for the 
varying gain as the spacecraft sweeps across the beam Since 
the average value of G(0)/G a for \0\ < 0 a /2 is 0 81, we shall 
simply insert this extra loss into Eq (8) This is a makeshift 
adjustment, of course 

The system that makes T d = T a can be obtained from the 
following iteration on 6 , starting with 6 = 1 


L 


5 b 2 

f 




17700 z r b 2 

ISO 

r = 

L 2 

r' = 0 81 GG pd (1/2, 6) r 


y (3i) 


6 =-£- 


r'L 2 

L * 0 +x 2 V 1 +2/ 


, 2/3 


J 


A T d = T a system leaves little room for error, there are at 
most two chances to acquire the signal We have computed this 
option simply to find out what can be done with a nontrack- 
ing parabolic antenna One might design an r mln system using 
a nontracking fan-beam antenna with the wider dimension of 
the beam oriented east to west The beamwidth would allow 
the spacecraft to stay in the beam for several word times The 
antenna would probably have to be adjusted daily in declina- 
tion 


To derive the r min (tracking) system we plot r vs 5 for the 
center and edge cases (Fig 3) The intersection gives 5 opt = 
1 59, r min = 0 925 Then L = 5 opt b 2 /f= 568 looks, T d = 
Lib = 568 sec The antenna parameters are G a = r min /r, so = 
1723, d = 2 30 m, 6 a = 3 56 deg, and T a = 427 sec 

The iteration Eq (31) gives the T d = T a (nontracking) 
system Parameters are 5 = 0 909, L = 325 looks, G a = 2967, 
d = 3 02 m, 6 a = 2 71 deg, and T a =T d = 325 sec 

2 X-band Start from /= 10“ 2 Hz/sec, (F > /A^ 0 ) 1SO = 6 03 X 
10~ 4 sec -1 The analyzer bandwidth is 1000 Hz, so b = 4 Hz 
Equation (13) with a y = 3 X 10 -1 1 , Tj = 0 5 sec gives = 
0 81 From Fig 2 we get the r min system 5 opt = 1 63, r m ln = 
0 409, L = 2608 looks, T d = 652 sec, G a = 2715, d = 0 791 m, 
6 a = 2 83 deg, and T a = 340 sec 

The parameters of the T d = T a system are 6 = 0 554, L = 
886 looks, G a = 6387, d = 1 21 m, 0 a = 1 85 deg, and T a = 
T d = 222 sec 


VIII. Conclusions 

Provided that certain assumptions about Doppler rate error 
and spacecraft oscillator phase stability are met, an alternating- 
tone signal from a spacecraft 4 AU away can be acquired by a 
ground antenna a few meters in diameter, a receiver with noise 
temperature of order 200 K, and a commercial 256-channel 
spectrum analyzer Although our design examples assume para- 
bolic antennas, other types such as fan-beam antennas or 
electronically steerable arrays should be studied, especially if 
one wants to avoid mechanical tracking 

Our analysis has included an estimate of the effect of 
oscillator flicker noise on MFSK detection A preliminary 
theory with a certain amount of heuristic argument has 
yielded the approximation Eq (13) for the spectral mean loss 
More work on this theory is needed to improve its accuracy 
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Table 1 Design parameters for ICE surveillance systems 



S-band, 

tracking 

S-band, 
non tracking 

X-band, 

tracking 

X-band, 

nontracking 

Doppler rate error 
/(Hz/sec) 

0 0028 

0 0028 

0 010 

0 010 

Spectrum analyzer 
bandwidth (Hz) 

250 

250 

1000 

1000 

Word time T d 
(sec) 

570 

330 

650 

220 

Antenna gain 
G a (dBi) 

32 4 

34 7 

34 3 

38 1 

Antenna 
diameter d (m) 

2 3 

3 0 

0 79 

1 2 

Half-power 
beamwidth & a (deg) 

3 6 

2 7 

2 8 

1 9 


Table 2 Design control table 




S-band 

X-band 

1 

RF power to antenna, dBm 




Transmitter power, dBm 

37 0 

40 0 


Circuit losses, dB 

-0 7 

-0 7 (’) 

2 

3 

High-gain antenna gain, dBi 
Space loss, dB 

28 5 

39 7 


loss = \ 2 I(4itR) 2 ,R = 5 8 X 10 8 km 
\ = 013m or 0 356 m 

-275 0 

-286 2 

4 

Power received by isotropic antenna, 
dBm (1+2 + 3) 

-210 2 

-207 2 

5 

Noise spectral density, dBm/Hz 

N 0 = kT,k = 1 38 X 10' 20 mW sec/K 

-177 5 

-175 0 


GaAs FET front end, K 
Noise from ground, cosmos, 

120 

220 


atmosphere, K 

10 

10 

6 

P/N 0 , isotropic antenna, dB sec"^ 
(4-5) 

-32 7 

-32 2 



(ALLAN VARIANCE) 
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GROUND SYSTEM 


Fig 1. Surveillance system block diagram 
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Fig 3 Finding the minimum SNR 
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Appendix A 

Loss From Oscillator Phase Jitter 


The foundations of the second-moment theory of processes 
with stationary nth differences are set out by Yaglom (Ref 8) 
These processes provide rigorous models for phenomena which 
appear to have spectral densities behaving like 1 //“ as /-*0, 
where a < 2n + 1 For our case, a = 3, so n = 2 suffices while 
n = 1 does not 

Let 8(f) (the phase of an oscillator) be a process with 
stationary second differences Its spectral representation is 

0(t)=a o +a l t+^ct 2 + [ l -F(u)icjt]^) 

1 J o° -W 2 

(A-l) 


6(t) has so little high-frequency power that 8(t) exists When 
one attempts to do calculations, one gets results that depend 
strongly on the high-frequency components of 8(f), which we 
know intuitively to be irrelevant to the situation The instanta- 
neous oscillator frequency is perhaps meaningless It is essen- 
tial to use the correct version Eq (A-l) with a cutoff function 
F(u>) 

We now compute an estimate of the mean loss incurred 
when measuring the spectrum of an oscillator with phase 8(f) 
We can assume that the output has been heterodyned to 
baseband and that we are measuring the power in a frequency 
band of width b centered at zero frequency Let bh(bt) be the 
time window of the measurement, where h satisfies 


where c is a constant, a 0 and a 1 are random variables, F(w) is 
a bounded function satisfying 

F(u>) = 1 + 0( |col), co -*■ 0 

= 0(l/|co|), co->±°° (A-2) 

and Z(co) is a complex-valued process with orthogonal incre- 
ments, related to the two-sided spectral density l h S gg (co) 
by 


E\Z(u’)- Z(co)| 2 



(A-3) 


h(x) = h(~x) > 0 

fh(x)dx = l,fx 2 h(x)dx < 00 
For later use define 

H(y ) = J e ,xy h(x) dx 

The mean spectral loss is given by 

G i = f$E e ,[0(s) - 6(t) h 2 h (bs) h (bt)dsdt (A-5) 


The spectral density must satisfy 


Assume that 8(f) is a Gaussian process Then 


L 


co 


1 + CO 




4 “00 


(A4) 


We have assumed that the measure given in Eq (A-3) (the 
spectrum of the generalized process 8 ( t )) is absolutely contin- 
uous 

Changing F( co) while retaining Eq (A-2) simply adds a new 
component to a } , which is just a constant of integration, a 
frequency offset If one formally integrates the spectral repre- 
sentation of a stationary process twice, one obtains Eq (A-l) 
with F(u>) = 1 Such a representation for 6(f) is valid only if 


where 


Ee ,{e(s)-e(t)\ = e X p ~D(s,t) 


D(s,t) = £'[0(*)- 0(t)\ 


Since b 2 h(bs)h(bt) dsdt is a positive measure with total mass 
1 , Jensen’s Inequality gives 


In pr < ff D(s,t) b 2 h(bs)h(bt) dsdt (A-6) 
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To first order in D(s,t ), the two sides of Eq (A -3) are equal 


(A-9) 


In computing D(s,t) from Eq (A-l) we assume that c = 0, 
tfj =0 The assumption c = 0 means that there is no average 
frequency drift The assumption a, = 0 means that we don’t 
want a constant frequency offset to affect the measurement 
Also, it is convenient to replace F(w) by F(oj/b) 

Let 

W(t ,<o) = e‘ u 1 - 1 - icof F 

Then 


D(s,t ) = I \W(s,io)~ W(t^)\ 2 S ge (o:)^ 

Jo 

and 

ln !f e l < , {f (A - 7 > 

where 

V(y) = j \e lxy - H(y) - ixyF(y)\ 2 h(x) dx 


We still consider the function F to be at our disposal 
Fixing y we find the number F(y) that minimizes V(y) This is 
done by projecting e lx v - H(y ) onto the vector ix in the space 
L 2 [/i(x)] The answer is 

^ * WW) <A4!) 


It happens that the function F given by Eq (A-8) satisfies Eq 
(A-2) The function V(y) given by Eq (A-9) is a high-pass 
frequency response, and the right side of Eq (A-7) is a 
“high-pass variance” of d(t) with cutoff frequency propor- 
tional to b 

Concerning the way we have played around with F, we can 
say only that it seems to be the right thing to do Since 
changing F only changes the offset frequency a, and the effect 
of frequency offset is covered by the loss G pa , we feel free to 
minimize the effect of the offset on G ; We want the average 
frequency during the measurement to be zero in some sense 
More work on this part of the theory is needed 

Now set 

h(y) = 1, \y\ < 1/2 
= 0, elsewhere 

■W") - 

co 

where 0(f) is a stationary process with bandwidth much 
greater than b, so that F(c o/b) passes most of the energy in 
0(f) Then Eq (A-7) becomes 

In ^ < a 2 (0) + — f V(y)^- (A-10) 

G , 2irb 2 J 0 y 3 

By numerical integration we find that the last integral is 1/36 
(approximately) Substituting for a 2 (0) and A from Eqs (11) 
and (12) we get Eq (13) 
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Appendix B 

Chi-Squared Tail Estimates 


Equation 26 4 3 1 of Ref 9 is used to express the tails of a 
noncentral (or central) x 2 in terms of the tails of a standard 
Gaussian For our situation there are v = 2 L degrees of free- 
dom and a noncentrality parameter X = 2 Lr (or zero) The 
approximation requires L > 50 Let X be a scaled version of 
the noncentral x 2 , define X* = [X - E{X)] /a(X), and let Z be 
a standard Gaussian Given a probability p define x and z by 



P = 


1 + 2 r' 
(1 + r ') 2 


(B-l) 


P [X* > x) = P {Z > z} = p 

Straightforward manipulation of the approximation in Ref 9 
gives the approximation 


To get an approximation for the lefthand tail P {X* < x) , 
replace x by ~x and z by -z in Eq (B-l) The righthand tail of 
X is heavier than the Gaussian tail, the lefthand tail of X is 
lighter 
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The Practical Limits of Photon Communication 
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We show that for photon communication, the rate R 0 = 1 nat per photon is the rate 
beyond which one encounters an explosive increase in both the required ratio of peak-to- 
average signal power and in the required bandwidth expansion On the basis of these 
results we conjecture that no practical photon communication system can be designed to 
operate above 10 nats per photon 


I. Introduction 

In a recent report (Ref 1), it was shown that the R 0 - 
parameter associated with noiseless optical communication 
using photon-counting techniques (hereafter we call this 
“photon communication”) is one nat per photon Now for any 
channel, R 0 is widely believed to be the rate above which the 
implementation of a reliable communication system becomes 
very difficult, but there is no really sound mathematical 
support for this belief In this paper, however, we will give 
rigorous mathematical substantiation to this “ R 0 conjecture,” 
for the special case of photon communication 

Roughly speaking, we shall prove that for photon com- 
munication, the rate R 0 - 1 is the rate beyond which one 
encounters an explosive increase in both the required ratio 
of peak-to-average signal power and in the required bandwidth 
expansion 

Precisely speaking, what we shall prove is this Let p denote 
the rate (in nats per photon) of a given reliable photon com- 
munication system, let a denote its required ratio of peak-to- 
average signal power, and let 0 denote its required bandwidth 


expansion factor Then necessarily, as we will show in Sections 
II and III, 


a > e p 1 - 1 

(1) 

e p_1 - 1 


P > 

(2) 


P 


Thus as p increases linearly, both a and /3 must increase expo- 
nentially On the basis of Hqs (1) and (2), we conjecture that 
no practical photon communication system can be designed 
with p > 10 On the other hand, in Ref 2 it was shown that 
one could design a practical system at about p = 3 using pulse 
position modulation and Reed-Solomon codes Thus, the gap 
between what is presently practical and what may someday be 
practical is reasonably small This is in spite of the fact that 
channel capacity (i e , the largest p that is theoretically pos- 
sible) is extremely large (Ref 3) 

II. The Poisson Channel Model 

We assume that any photon communication system works 
as follows The time interval during which communication 
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takes place is divided into many subintervals (“slots”), each of 
duration t 0 seconds The transmitter is a laser which is pulsed 
during each time slot, it may be pulsed with a different intens- 
ity in each slot At the receiver is a photon counter, which 
accurately counts the number of photons received during each 
time slot We denote by jc, the expected number of photons 
received during the ith time slot, jc ( will be called the intensity 
of the zth pulse 

It may be that “noise photons” are present in such a 
system, but in many cases of practical interest, noise photons 
are extremely rare (For example, in a careful analysis of a 
potentially practical system, Katz estimated the rate of arrival 
of noise photons to be around 1CT 3 per second ) In any event 
we shall make the assumption that no noise photons exist In 
this case, because of the Poisson nature of photon arrivals, 
the probability that exactly k photons will be received during 
a slot in which the laser was pulsed with intensity x is e~ x x k /k ] 

Thus described, the optical channel is a discrete memoryless 
channel with input alphabet equal to the set of nonnegative 
real numbers (the possible values for the intensities x t ), and 
output alphabet equal to the set of nonnegative integers (the 
possible outputs of the photon counter) If a real number x 
is transmitted, the probability that the integer k will be 
received is given by 

p( k\x) = e~ x ~ (3) 

We call the channel described by Eq (3) the Poisson channel 

A code for this channel is a set of vectors = (jc z , , 
x m ), i = 1, , M, of length n Each component x l; is a non- 

negative real number and represents an intensity of the trans- 
mitting laser Assuming that each component of a code word 
requires one time slot for transmission, the rate of such a code 
is 

log 4/ 

R = nats per slot (4) 

n 

On the other hand, each component x tJ represents an average 
number of (received) photons, and so the code’s rate in nats 
per photon is 

p = — nats per photon, where (5) 


p = x^jnM, photons per slot (average) (6) 

The rate R in Eq (4) is a measure of “bandwidth expan- 
sion ” If, for example, we are transmitting at a rate of A nats 
per second, using a code of rate R nats per slot, it follows that 
we require A/R slots per second Thus the slot rate is equal to 
the nat rate multiplied by the factor l/R We thus define 

B — — — — slots per nat (7) 

log M 

and call the bandwidth expansion factor 

In Eq (6), p represents the average pulse intensity, i e 
the average number of photons received per time slot At a 
given frequency this number is proportional to the average 
received power On the other hand, the quantity 

L = max x (8) 

u ‘i 

represents the maximum required pulse intensity required by 
this particular code, and it is proportional to the maximum or 
peak received power We denote the ratio of L to p by a 

max (x ) 

/ IJ ' 

a = — = (9) 

P 1 

n • M 2 -j x u 

i,i 


We have now precisely defined the quantities a, (3, and p 
mentioned in Section I Our proof of inequalities in Eqs (1) 
and (2) rests on certain information — theoretic results about 
the Poisson channel, which we now describe 

We would like to compute the information — theoretic 
capacity of the Poisson channel, whose statistics are given in 
Eq (3), that is, the maximum possible mutual information 
I(X,K), where X is a nonnegative random variable, and A” is a 
nonnegative integer-valued random variable related to X by 
the conditional probabilities However, it is very easy to see 
that this maximum is infinite (take X to be a discrete random 
variable which assumes a very large number of values which 
are very far apart) To get a meaningful problem, we must 
restrict X somehow The most natural restriction (Eq (6)) 
is to fix the mean of X, and so we define 

C{p) - sup [I{X,K) E(X)=p] (10) 
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(17) 


According to Shannon’s noisy-channel coding theorem (see 
Ref 4, Chapter 7), C(p ) represents the maximum possible rate 
(in nats per channel use) of a reliable photon communication 
system which is restricted to operate at an average of /a photons 
per slot or less 

A second possible restriction is on the maximum value that 
X can assume (Eq (8)) Thus we define 

C(ji,L) = sup [I(X,K) E(X) = n,X<L] (11) 

Again, according to Shannon’s theorem, C(p,L) represents 
the maximum possible rate for any photon communication 
system which is restricted to operate at an average of <p 
photons per slot, and a maximum of L photons per slot 

In Section III, we will prove the following results about 
C(p) and C(p,Z, ) 

C(ji) < log (1 +p) + plog (l +-M (12) 


p < (e p ~ 1 - If 1 

Now since the bound on the nght of Eq (15) is easily seen 
to be an increasing function of /a, it follows from Eqs (15) and 
(17) that 


but7?=^ _1 This proves Eq (2) 

To prove Eq (1) observe that Eq (13) tells us that a >e? 
for L < 1 This is stronger than Eq (l),thusEq (1) can only 
fail for L > 1 But if L > 1 , the ratio a = L/p is by Eq (17) 
greater than e^ -1 - 1 This proves Eq (1) 

It now remains to prove Eqs (12) and (13) This we do in 
the next section 


C(ji,L) < p log — if L < 1 (13) 

P 

In the remainder of this section, we will show how Eqs 
(12) and (13) can be used to prove our main results, Eqs (1) 
and (2) 

First, note that by the converse to the noisy-channel coding 
theorem, the rate R of a reliable communication system which 
operates at an average of p photons per slot is bounded by 

R < COi) (14) 


III. Proof of Eqs. (12) and (13) 

In this section we shall give proofs of the important in- 
equalities of Eqs (12) and (13) We begin with Eq (12) 

By definition, C(p) is the largest possible value which can 
be assumed by the mutual information I(X,K ) when the test 
source X is restricted to satisfy E(X) = p But by well-known 
results on mutual information (see Ref 5, Chapter 1), 


I(X,K) = H(K)-H(K\X) 
<H(K ) 


(19) 


Now C(p) itself is bounded by Eq (12), and since log (1 + /a) 
< p, we have 


R < p 


1 +log 



(15) 


The rate of this system measured m nats per photon is, by Eq 
(5), R/n, and so 


where in Eq (19) H(K) denotes the entropy 2 k p k log p k l of 
the random variable K Since for the Poisson channel, 
E{K\X) = X, it follows that E(K) = E[E(K\X)] = E{X ), and 
so the random variable K has the same mean as X, viz , p 

The problem of maximizing the entropy of a nonnegative 
integer-valued random variable with a given mean is easily 
handled with standard information-theoretic techniques 
Omitting the details (see Ref 5, Problem 1 8), the result is 


p < 1 + log 



(16) 


H(K) < log (1 + p) + p log 



( 20 ) 


for p > 1 a simple algebraic manipulation of Eq (16) yields 


(provided E(K ) = p, K assumes only nonnegative integer 
values) Equations (19) and (20) combine to give Eq (12) 


65 



k > I 


(25) 


We turn now to Eq (13), which lies somewhat deeper The 
problem is to maximize I{X,K ) given that the distribution of 
X is restricted to [0,L] and has mean p 

We introduce the notation (see Eq (13)) 

x* 

p k (x) = e~ x — (21) 

to denote the probability that k will be received given that x 
is transmitted We define 


g k ^ g k - 1 g k+l’ 

The function F(x) of Eq (23) can be written in the form 

F(x) = x(log x - 1) + ^2 p k (x ) log j- + X(x - p) - C 
/t=o g x 

(26) 

If we differentiate F(x) twice, using the formula p' k (x ) = 
P k - , (x) ~ p k (x), we get 


P k = E \P k (*)\ 


(22) F"(x) = ~ + lP k - 2 ^~ 2 P k -\^ +p k ^ log 

k k 


the expectation being with respect to the distribution of X 
The quantity ~p k represents the probability that k will be the 
channel’s output, given that the input is the random variable 

Using standard techniques of convex analysis, it is now 
possible to show that a particular distribution confined to 
[0,L] with expection p maximizes I(X,K) if and only if for 
some constants C and X, 


x 


+ 


£ 


P k (x) log 


2 

s k+ 1 


g k g k+ 2 


e x F"(x) 


- E 




k = 0 


k°k + 2 


b k+ 1 


(27) 


(28) 


” P k ( x ) 

F(x) = 2J P k ( x ) l og — — + X(x-p)- C<0,0<x<h 

k = 0 P k 

(23) 

where equality holds in Eq (23) at all points of support of 
the distribution (Eq (23) is essentially the same as Theorem 
4 5 1m Gallager (Ref 4) The only differences are that our 
channel has a countably infinite alphabet, rather than a finite 
one, and that we have an extra constraint E(X) = p, which 
necessitates the Lagrange multiplier term X(jc - p) However, 
the modification of Gallager’s analysis needed to arrive at 
Eq (23) is quite easy, and we omit it ) 

Equation (23) is a very strong condition that must be 
satisfied by an extremel distribution For example, we use it 
to show the following 

Lemma For any L,p, a maximizing distribution can have mass 
at at most one point on (0,1) 

Proof Define 

g k = k'P k = E(e~ x X k ) (24) 

Then by Schwarz’s inequality [E{X iX 2 )] 2 < E(X J) E(X^) 
applied to X t = e ~ x / 2 X ( F~ l '>! 2 , X 2 = e~ x ! 2 X ( - k+ 1 )/ 2 , we 
have 


Each coefficient in the series Eq (28) is by Eq (25) non- 
negative Since e x /x is a decreasing function for 0 < x < 1 , 
it follows that e x F”(x) is decreasing in this range also Since 
e x F"(x ) is positive at x = 0+, e x F"(x), and so F"(x) also, can 
be zero for at most one value of 0 < x < 1 

An extremel distribution by Eq (23) must have F(x) < 0 
for all 0 <x < L Since F(x) is differentiable, it follows that 
F’(x) = 0 whenever F(x) = 0 Now if say /^(Xj) = F(x 2 ) = 0 
with 0 < x, < x 2 < 1, then there exists x, < x 3 < x 2 with 
F'(x 3 ) = 0 By the above remarks F’(x l ) = F\x 2 ) = 0 as well 
This in turn implies the existence of x 4 , x 5 x t < x 4 <x 3 < 
x s < x 2 with F”(x 4 ) = /’"(xj) = 0 But we argued above that 
F"(x) could vanish at most once on (0,1) Thus, / r (x) can 
vanish at most once on (0,1), i e , the optimizing distribution 
can have mass, at most, at one point in (0,1) This completes 
the proof of the lemma 

We now use the lemma to prove Eq (13) Since L ^ 1, the 
lemma tells us that a distribution on [0,L] with ^(x) = p can 
have mass only at x = 0 and one other point x = C < L We 
shall complete the proof of Eq (13) by showing that for any 
distribution concentrated at x = 0 and x = L the resulting 
mutual information I(X,K ) satisfies 

I(X,K) < p log ~ (29) 
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(It is sufficient to take C = L because the right side of Eq (29) 
is an increasing function of L ) 


F{q) = (l-Q) log Y^+(G-^)log|-L(q-0 


Thus, let X have distribution 

P(X = 0) = p 

P(X = L) = q, p + q = 1 


and wish to show that f(q) < 0 for 0 < q < 1 This is easily 
seen, given the following (whose straightforward verifications 
are omitted) 

/(0)=/(l) = 0 (32) 


Then a straightforward calculation yields 

I{X,K) = (1-0 log + Q log | - L(q - Q) 

Q=q( \~e~ L ) (31) 

Our goal is to show that this quantity is < p log (L/p), but 
p = E(X) = qL So we must show that the right side of Eq (31) 
is < qL log 1 /q Subtracting this quantity from Eq (31), we 
define for a fixed L > 0 


/'(0) = -»/'(!) = « (33) 

f"(q 0 ) = 0 °n!y for q = ~ - j (34) 

1 - e~ L L 

Because /(0) = 0, /'(0) = f(q) is negative for all suffi- 
ciently small q If now f(q) = 0 for 0 < q < 1 , it would neces- 
sarily follow that /' would vanish at two interior points of 
(0,1) Since also /'(l) = 0,/” would vanish at two pomts of 
(0,1), contradicting Eq (34) Thus f(q) < 0 for all 0 < 1 , 

and this completes the proof of Eq (13) 
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The classic framework for analysis ultimately follows a senes of rules whether or not 
those rules are consciously applied The rules usually take the form of define the 
problem ( system j, define a measure of effectiveness, generate alternatives based on 
analysis, and, weigh and decide Many methods have been developed and many analytical 
and quantitative techniques and tools designed for general applications to different 
models This paper descnbes the approaches used in developing the basis for a system 
design evaluation model Two well known evaluation cases are modeled, leading to a 
general model that is then detailed for application specifically to the Navigation Network 
conceptual designs Identification of the system attributes and their associated perspec- 
tives and the method of scoring and ranking a candidate design are discussed 


I. Introduction 

In the conceptual phase of development of any new system, 
several approaches are taken to ensure the best possible 
solution of the design problem, several because there is simply 
no a prion way of knowing which approach will best meet the 
user’s need within a defined set of constraints Since it is not 
cost effective to expend resources in constructing different 
prototypes to physically (actually) test each of the 
approaches, evaluations need to be performed during the paper 
stage of design 

There are two factors that make evaluation difficult First, 
although attributes such as growth potential and operability 
tend to be common to all purpose-designed systems, all such 
systems are not alike Therefore, the capacity to view each of 
the attributes in their true relations or relative importance 
assumes a system uniqueness that makes their definition 
elusive and time-consuming, especially smce both the attri- 
butes and their perspectives must have the support of a 


consensus of those persons representing the organizations that 
ultimately have a stake in the new system These are among 
the matters that must be dealt with at the engineering level of 
the design organization 

The second factor is that decisions regarding the choice of 
approach to be used in building the system have to be hved 
with for many years after implementation Modifications are 
expensive and the downtime, for example, may be unaccept- 
able to the system user In this context, the choice, or 
decision, is made at the management level but the quality of 
the choice is heavily dependent upon the degree of profes- 
sionalism applied in the development and application of the 
system evaluation criteria at the engineering level 

II. Considerations and Requirements 

One of the tasks defined in the pre-Project phase of the 
Navigation Network Study was the determination of a set of 
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evaluation criteria for the candidate designs, which were to be 
developed under yet another task 

The principle purpose of the evaluation criteria was to 
provide Project Management with an acceptable means of 
judging the merit of each candidate design and an acceptable, 
uniform method of scoring and ranking each Subsidiary 
purposes were to 

(1) Assist the Design Study Team by suggesting topics to 
consider in their efforts to produce acceptable candi- 
date designs 

(2) Assist in identifying the level of work necessary and the 
amount to be done in completing a candidate design 
through the research/reference stage to provide the 
data needed to respond to relevant evaluation criteria 

(3) Give DSN Management, in association with the Project 
Report on how designs are ranked, another input to use 
in reaching a decision whether or not to proceed with 
the Navigation Network implementation 

The requirement that the evaluation criteria had to be 
developed separately and independently from the candidate 
designs themselves posed an initial difficulty in getting a useful 
word picture, or model, of just what would constitute a 
candidate design No written, formal definition was available, 
so a definition was developed and given for practical purposes 
as a “functional design ” To expand this a bit further, the 
planning configuration design used to develop future budget 
line items for the Navigation Network served as a useful 
example The planning configuration is described here simply 
as a set of configured generic hardware, and for the sake of 
completeness as receiving instruments consisting of a 10- or 
20-meter antenna, a cryogenic front-end amplifier, the VLBI 
receiver and sampler, a frequency and time standard 
(H-maser), a control computer, and the communications 
equipment (DOMSAT) 


III. Model Development Approaches 

Three different approaches were investigated, initially, in an 
attempt to find an acceptable model of the evaluation process 
an engineering approach, a library research approach, and a 
strawman approach 

A. Engineering Approach 

In the engineering approach, a survey was made of the work 
other persons in the DSN had done in the area of evaluation 
criteria The results of the survey were so meager they could 
not be used 


B. Library Approach 

Concurrent with the engineering approach, a library 
research effort was initiated to determine the availability of 
useful literature on evaluation criteria that could be obtained 
and used within a reasonable time period This research also 
proved fruitless in revealing how an evaluation criteria process 
might be developed 

C. Strawman Approach 

The criteria for the strawman approach were developed 
from suggestions made by various DSN personnel The criteria 
were partitioned under topic headings, for example, Design 
and Operation, using very flexible rules to afford the sem- 
blance of order that was needed because the randomness of 
criterion input was not conducive to their falling into cogni- 
tive relationships In retrospect, the strawman criteria can be 
viewed as a set of unstructured attributes mixed-in with sys- 
tem requirements 

The strawman criteria were sent to a number of managers 
who had fiscal, techmcal, or user mterests in the proposed 
Navigation Network, their comments were solicited, and 
responses received were carefully reviewed and were inter- 
preted to mean that 

(1) A series of filters or gates were needed to weed out 
candidate designs that were unresponsive to the desired 
attributes 

(2) Cost and risk were elements of the evaluation process 
that should be treated as parameters and not as 
attributes 

An analysis of the comments received pointed to the 
inadequacy of the strawman approach to develop evaluation 
criteria, however, the comments nourished the thought that a 
generalized case model would be feasible 

IV. Definition of a Case Model 

The basic ground rules for a generalized case model were 
established as simply as possible The case model 

(1) Should not be associated with the DSN 

(2) Had to be sufficiently complex to require more than a 
superficial effort to develop 

(3) Had to be one that most people could readily identify 
with and understand in order to sustain their interest 

The methodology was to distil the underlying philosophy 
of evaluation from the case model and to adapt it to the 
current need With the ground rules in mind, two cases were 
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identified where individuals make large expenditures in rela- 
tion to their personal income that require prudent considera- 
tion at each stage of their involvement the purchase of an 
automobile and the purchase of a home 

It was decided to look at both cases to determine what 
concepts were involved in the respective evaluation processes 
Since the Navigation Network could conceivably involve the 
use of new equipment at new locations, it was further decided 
that the cases to be defined should be the purchase of a new 
automobile and a new custom home In citing the need for 
new and custom items, it was felt that a greater effort would 
be required and thus result in a more definite process than 
cases involving the upgrading of an existing automobile or 
home 

Although both case study results are available to the 
interested reader, only one — the purchase of a new auto- 
mobile — is described further for the purposes of illustration 

A. Purchase of a New Custom Home 

Not described further due to time and space limitations 

B. Purchase of a New Automobile 

1 Requirement. The requirement for a new car has to be 
formulated in general terms as a first step in satisfying the 
need As an example, the requirement could be for a sedan in 
which a family may be transported on short urban journeys 
with economy of operation, after incurring a certain initial 
cost 

The first logical step is to obtain automobile marketing 
data, which may be accomplished by obtaining brochures of 
various models handled by local dealers representing several 
different manufacturers At this point the prospective pur- 
chaser goes through a preliminary elimination or weeding-out 
process of candidates by reviewing the contents of the 
brochures This will eliminate automobiles that 

(1) Are not sized for that particular family 

(2) Have a fuel consumption that delivers fewer miles per 
gallon than required 

(3) Exceed the upper cost limit 

2 Technical constraints Models that meet the initial 
requirements are then subjected to further scrutiny based 
upon technical constraints Here the criteria will reflect 
personal experience and may perhaps include 

(1) Number of years the model has been available to the 
public without significant modification or manufac- 
turer’s recall to rectify design defects 


(2) Resale value of earlier version of the model 

(3) Opinions of other users of a particular model who are 
qualified in the eyes of the prospective purchaser 

Additional models may be eliminated at this point through 
their subjective failure to meet all, or most, of the technical 
attributes considered relevant by the prospective purchaser 
Those candidates that survive are given some preliminary order 
of preference, based on how well they subjectively meet the 
technical attributes overall 

3 Optimization evaluation Each of the remaining candi- 
dates then goes through an optimization evaluation to deter- 
mine how closely certain needs are satisfied These needs 
might include 

(1) Driveability 

(2) Passenger comfort and convenience 

(3) Performance improvement 

(4) Appearance 

(5) Safety 

The candidates are then ranked accordingly to how well all the 
criteria have been satisfied by different manufacturer’s models 
and option combinations 

4 Cost In parallel with the optimization process, the cost 
of ownership is being developed Costs include 

(1) Cost of basic vehicle and selected options 

(2) License fees and taxes, annual 

(3) Insurance, annual 

(4) Estimate of annual fuel costs, based on some arbitrary 
distance 

(5) Maintenance costs, annual 

5. Additional considerations The final choice of vehicle, or 
even a decision to buy one at all, may depend upon additional 
considerations such as 

(1) Availability of the model in the desired configuration 
at the needed date 

(2) Dealer concessions on the total price 

(3) Availability and cost of financing under different 
purchasing plans 

(4) Ability of the prospective buyer to dispose of his 
existing vehicle at what he considers to be a fair price 
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(5) A willingness to relinquish or defer one or more 
options to meet cost constraints 

(6) Opportunity costs 

C. Workability of the Case Model 

The workability of the two case models must be demon- 
strated, the best way would be to use them in the respective 
process of new car and custom-home purchasing Since neither 
activity is feasible, the next best thing is to obtain the opinion 
of experts as to the usefulness of the case models 

1 Expert opinion, model 1, new custom-home This case 
model was first reviewed separately by two practicing civil 
engineers They both agreed that the process was acceptable 
but recommended changes and additions to the evaluation 
criteria, which were included in the model 

The revised model was then reviewed by a mechanical 
engineer who is in the process, personally, of constructing his 
own home He agreed that the process was essentially correct 
and that although he thought that some more evaluation 
criteria might be included, the model was adequate for the 
present purposes 

2 Expert opinion, model 2, new automobile This case 
model was reviewed separately by two managers who had 
recently purchased new vehicles for themselves Each one 
agreed that the process was essentially the one they had gone 
through, although neither one had written it down on paper 
Both remarked that they would use some additional criteria 
and these were added to the model used in this report 

Based upon the interview statements, it was concluded that 
the evaluation processes of the two case models were workable 
and useful 

D. Word Representation 

A word representation of the evaluation process was 
developed for each of the two case models A comparison 
revealed a close degree of correspondence m their respective 
structures The result of that comparison may be characterized 
as follows 

(1) Requirements have to be expressed in terms of 

(a) What has to be done 

(b) How well the purpose has to be done 

(c) A ball-park cost limit expressed m terms broad 
enough to invite effective competition 


(2) A wide range of candidate designs, products, prop- 
erties, and services already exist to satisfy the most 
discriminating of needs The market has to be re- 
searched for those that might quahfy for candidacy 

(3) An initial evaluation is necessary to eliminate those 
candidates 

(a) That do not meet the functional and performance 
requirements 

(b) That exceed the upper cost bound 

(4) Survivors of the preliminary requirements evaluation 
then have their capabilities subjected to a technical 
attributes evaluation where some feehngs as to order of 
choice are expressed and some final candidate rejec- 
tions may be made because the nsks and inconveni- 
ences outweigh the benefits and conveniences of 
ownership 

(5) The initially ranked candidates are studied in detail and 
the optimization of the attributes within the physical 
constraints imposed by a particular design are evalu- 
ated Then for each feature of design, the evaluation 
results of the candidates are compared and a final order 
of technical preference is established 

(6) In parallel with the optimization evaluation process, 
the fixed, one-time costs and the continuing costs are 
estimated to an accuracy that is consistent with the 
available data 

(7) The decision to procure any of the final candidates is 
developed from weighing different financing plans, for 
example, outright purchase, deferred payments, leasing 
or rental, and constraints (a set of evaluation criteria 
reflecting the state of family affairs) imposed by the 
business climate agamst the technical capabilities to be 
offered in each case The choice may also be influenced 
by the ability of any of the final candidates to tolerate 
changes of certain feature optimizations in order to 
meet cost limitations while still being able to perform 
the essential task 

Alternative representations are often useful in developmg 
additional perspective on a problem Figure 1 shows a top-level 
block diagram of the generalized evaluation criteria process, 
derived from the word representation of the case models 

In practice, the expansion of each of the block diagram 
boxes to give increasing detail of the internal processes will be 
dependent upon the point in the diagram that the action has 
reached Initially, it is desirable to minimize the work to be 
done because the number of candidates will be greatest 
Therefore, it would be expected that the processes at the top 


71 



of the block diagram would be simplest, le , the go/no go 
kind 

At the next level, more judgement and consideration would 
be needed in evaluating each attribute of the candidate, with 
rejection still being a more or less probable result At the 
technical optimization evaluation and life cycle costing stage, 
however, more detailed processes would be expected since 
more precise data need to be generated to satisfy the desire to 
“get the best for the least ” It is also conceivable that the 
block diagram process as a whole may be repeated more than 
once if some of the criteria applications reveal marginal 
technical capabilities that can be relieved only by changing the 
functional, performance, or initial cost requirement 

V. Definition of the Evaluation 
Criteria Process 

Figure 2 shows the top level flow chart that evolved from 
the task of developing the Navigation Network evaluation 
criteria process from the general model It has the same 
essential flow of data as shown in Fig 1 The major change has 
been to differentiate between the processes and decisions to 
emphasize the progressiveness of the acceptance or rejection of 
candidate designs 

The flow chart (Fig 2) may be regarded as being divided 
into four stages in series The first stage consists of the func- 
tional requirements, performance requirements, and upper 
cost limit gates Simple decisions are made here on whether a 
candidate design will be able to do all things it is supposed to 
do and do all of those things as well as it should at an 
affordable cost 

In the next stage, the remaining candidates are evaluated 
and scored for their compliance with Technical Attributes and 
given a preliminary ranking If there is an obvious dichotomy 
in that ranking, those candidates in the lower subset will be 
eliminated at the Technical Ranking gate Candidate design 
rejections occurring at this point will further reduce the 
amount of detail work that has to be done in the third stage 
Flere the survivmg candidate designs go through two parallel 
processes 

In the one, they are evaluated and scored on how well the 
competing needs of Implementation and Operations are 
satisfied, and the Intermediate Ranking that follows is 
critically reviewed to determine if there are any obvious 
inconsistencies in that ranking, caused by sconng subjectivity, 
that need to be corrected Any subsequent ranking changes are 
justified and documented In the second parallel process, the 
life cycle cost of each candidate is estimated to a level of 
accuracy consistent with the available data 


These two data components constitute the output require- 
ments of the Navigation Network study and are two of the 
mputs to the Management Review process Other inputs 
include the alternative funding plans that are created to meet 
the different costs of each of the candidate designs and 
business constraints, such as the choice of performing the 
detail design of the selected system in-house, or to let a system 
mode contract 

To define the system attributes needed at the Technical 
Evaluation stage, the various constraints to be placed upon 
each candidate design were considered They were loosely 
identified as Implementation, Operation, Management, and 
User, a general question was formulated in each case to relate 
them to the Navigation Network Typically, for the Operation 
constraint the question was, “How will this design affect 
existing operational activities and procedures 9 ” Answering this 
question made it easier to see the relationship between the 
constraints and the system attributes, as well as between the 
constraints and the general model, so as to be able to place 
them m the model accordingly 

Each attribute needs one or more perspectives in the shajie 
of a question that will act as an evaluation criterion or 
discriminant This part of the task is difficult to do because a 
thorough knowledge of each of the constraints is needed So 
both the attributes and the perspectives have to be developed 
and refined through iterative discussions with the numerous 
persons possessing the appropriate knowledge and an interest 
in the quality of the final design 

To define the system attributes and perspectives for the 
third or Technical Optimization stage (see Fig 2) the con- 
straints considered were Implementation, Operation, and Cost 
The question used, for example, to relate the Implementation 
constraint to the general model was “Is this candidate system 
easy to design 9 ” 

An essential part of the evaluation process is the recording 
of how effectively a candidate satisfies a particular criterion 
The purpose of setting score values for the evaluation criteria 
is to provide the evaluators with a first indication of how the 
various designs rank against each other To decide on the score 
values, the set of evaluation criteria has to be complete, that is, 
the system attributes to be evaluated and the wording of the 
associated perspective statements must have had consensus 
approval before development of the sconng can begin The set 
of evaluation criteria have to be searched for a theme that can 
be used to partition them into subsets If more than six subsets 
result from this process it is probably too many and another 
theme should be considered In the Navigation Network 
evaluation task, four different levels of risk were conceived 
and used to score the evaluation catena as follows 
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(1) Criteria that must have a positive response smce failure 
to satisfy them could mean failure of the design to 
meet the user need 

(2) Criteria related to those areas of relative inflexibility, 
and therefore high risk, that could impact costs 
severely if it turns out later that the wrong decision had 
been made For example, 

(a) Key equipment charactenstics that, once fixed, 
cannot be changed without impacting operational 
requirements, commitments, and costs 

(b) User costs 

(3) Criteria related to Operations where the risks may be 
considered as next lower, because of the potential for 
flexibility and work-around 

(4) Criteria that are considered to be related to the lowest, 
most controllable nsks 

Clearly, four scormg values are needed There is, however, a 
second condition to be satisfied Some of the perspectives are 
worded to determine whether the candidate design will be 
advantageous, neutral, or disadvantageous in terms of the 
attribute being assessed In such instances the scormg value is 
given a plus and minus range 

Those criteria concerned with the highest risk were given a 
yes/no response The next level was given a score of 10 No 
attempt has been made to differentiate between the criteria in 
this subset in order to establish the range of seventy of risk It 
is felt that if any such differentiation becomes necessary it will 
occur naturally dunng the ranking review when each case will 
be argued on its own specific merits The subset of criteria 
whose level of nsk is associated with “flexibility” and 
“workaround” is given a scormg value of 6 while those criteria 
associated with the lowest level of risk are given a scormg value 
of 3 

The ranking of the Navigation Network design evaluations 
is to be carried out in three phases Candidates that reach the 
Technical Attributes Evaluation stage are evaluated and 
scored Their scores are then individually totaled and the 
preliminary rank is decided on the relative total score so that 
the highest scoring candidate will be ranked first A similar 
activity takes place at the Technical Optimization Evaluation 
stage Here the scores for a candidate design are added to those 
it obtained at the previous level of evaluation to produce a 
total score on which to base the intermediate ranking This 
ranking is not the final arbiter A review of how each score is 
obtained and how it influences the ranking may uncover 
doubtful or unacceptable conditions existing among the rank 
leaders and underrated or highly attractive conditions lurking 


among the rank trailers This situation can only be resolved by 
a round of evaluation that concentrates on the perceived 
inequities created by the first round This review will be 
conducted by the Navigation Network Project team Their 
determination of the ranking with the reasons for any changes 
that they make will be presented to the Techmcal Steering 
Committee with an open forum for comment Further reasons 
for ranking changes made at this meeting will be considered by 
the Project Manager before he makes his ranking recommenda- 
tion to the DSN management 

Some dissimilarities have been noted between the general 
model evaluation process and the Navigation Network evalua- 
tion process, these are 

(1) Detailed designs for new homes and production models 
of new automobiles are readily available from which 
the individual may make a choice, whereas the Naviga- 
tion Network functional designs have yet to be 
developed This dissimilarity may be disregarded, how- 
ever, because if no candidate designs are develojaed the 
Evaluation Criteria will serve no useful purpose 

(2) The intent to procure a new home or automobile has to 
be of a high order or the individual will not undertake 
such a time consuming task, whereas the decision to 
build a Navigation Network is an institutional, rather 
than an individual, decision Smce the decision is 
outside the scope of the Evaluation Criteria this 
dissimilarity may be disregarded 

(3) There is an opportunity m the Navigation Network 
evaluation process to modify candidate rankings by a 
review procedure that was not observed in the general 
mode There are two possible reasons for this First is 
that an individual who is evaluating automobiles or 
home designs is most hkely doing that job with little 
formality, the “scoring” being done mentally So the 
need for the ranking review would not arise Second, 
one or two persons with fairly identical interests will be 
concerned in an automobile or new home purchase 
whereas a number of users having considerable diversity 
of interests are concerned in getting the greatest group 
use out of a unique facility such as the Navigation 
Network 

(4) The costing of a new home or automobile at the detail 
level is a very much simpler task than costing a new 
design for a Navigation Network at the conceptual 
level Also, the differences in cost are of several orders 
of magnitude greater for the Navigation Network 
compared with the car or home 

Table 1 lists the complete set of Evaluation criteria that will 
be used in the evaluation process of the Navigation Network as 
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illustrated in Fig 2 Considerable care was exercised in writing 
the various perspective statements As a result they are 
concise, and to stimulate discussion during their development 
a short note was written about each one to explain the 
background and to suggest how a candidate design under 
consideration might be scored A subset of these notes that 
refer to the Implementation constraint are appended below 

A. Maturity 

The purpose of the question is to find out (a) if the theory 
of the basic principles of the design has been developed to the 
satisfaction of the scientific community concerned with such 
matters, and (b) if these principles have been adequately 
demonstrated in practice through the use of purposefully 
designed hardware/software A “Yes” answer to the question 
should get the maximum score There is also a case for 
considering the question to rate “Yes” or “No” only, with 
acceptance of the design for a “Yes” response and rejection 
for a “No ” This approach, however, has not been selected for 
these Evaluation Criteria 

B. Buildability 

If the design under consideration is considered to be so 
specialized that a qualified electronic equipment manufacturer 
would need special facilities, tools, or skills to do the job, then 
the design should be given a zero score If, on the other hand, 
it is considered that the design can be built commercially with 
a very low probability of a cost overrun or late delivery, or 
both, then the questions should receive the highest possible 
score 

C. Practicability 

If the baseline geometry permits view periods of sufficient 
duration to enable the ground station equipment to acquire 
necessary and sufficient data to satisfy the needs of all the 
expected concurrent users, then the question should receive 
the highest score Otherwise, the specific design under con- 
sideration should receive a lower score 

D. Testability 

The purpose of the first question is to ensure that the 
design, if built, can be tested simply and easily whenever 
necessary to ensure that the capabilities are within specifica- 
tion In addition, the testing has to be accomplished with 
ordinary test equipment used by trained maintenance techni- 
cians in a short time span at low cost There is a strong desire 
to be able to test the design performance by independent 
means This could be accomplished by using other DSN 
equipment or the facilities of other organizations These 
independent verifications, therefore, will have to be done 
infrequently because of the coordination effort required and 
the extra cost involved If the answer to either question is 


“Yes,” then the design should be awarded the highest score 
For a “No” answer the score should be zero 

E. Compatibility 

Integration of the Navigation Network with the DSN could 
disturb the planned or implemented compatibility of the 
major interfaces These interfaces involve equipment, software, 
RF Imks, technical information, data and operations Some of 
these interfaces will be of interest only if a Project is in a flight 
phase, if a Project is in any pre-flight phase then all of the 
interfaces would be of interest If any interface 

(1) Will remain unchanged, then a zero score should be 
awarded 

(2) Is changed in a beneficial way, then it should be given a 
positive score 

(3) Is changed in a way that creates any new burden or 
imposition, then it should be given a negative score 

F. Growth Potential 

Past experience had shown that DSN equipment has been 
used in a capacity or has served a purpose that was not 
perceived at the time of its original design On other occasions, 
equipment has been modified to do additional tasks without 
detriment to the original task But here we are concerned with 
equipment in the conceptual stage and yet to be built If we 
know that the capability has been deliberately limited for 
reasons of time or cost, we know that there is growth 
potential If we can observe a recognizable innovation rate in 
the technical field of interest, we can expect growth potential 
Those designs that have an identifiable growth potential 
should receive a positive score, whereas those with obvious 
growth limitations should receive a negative score Designs 
with no identifiable growth potential and no obvious limita- 
tions should be given a zero score 

G. Schedule 

At the conceptual level, the only schedule data available is 
the PERT/CPM chart maintained by the Project Office If that 
chart shows completion on schedule for a particular design 
then a maximum score should be given 

H. Safety 

The question here is not whether the equipment can be 
designed to meet federal or state safety standards, since that is 
mandatory The question is how fault tolerant can the 
equipment be designed to localize failure damage, and will the 
risk to life be no worse than already experienced within the 
DSN, keeping m mmd the need to implement a highly 
automated system A design meeting such requirements should 
be given the maximum score 
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Table 1 Navigation Network evaluation criteria 


Navigation Network evaluation criteria requirements 

Attribute Perspective Response 

Functional Can the measurement types required to navigate all missions in the mission set be obtained 9 Yes/No 

Performance Can the most stringent accuracy requirements of the ground based system be attained 9 Yes/No 

Cost Is the preliminary cost estimate withm the limitation set by management 9 Yes/No 

Technical attributes evaluation process 

Implementation ^ 

Attribute Perspective Score range 

Maturity Is the technology of the design well understood 9 0-+10 

Buildability Can the design be built in a production environment to meet the requirements 9 0 - +10 

Practicability Can the design acquire data often enough and long enough to meet user requirements 9 0 - +10 

Testability Can the design performance be validated 

Simply and repeatably 9 0-+10 

By an independent technique 9 0-+10 

Compatibility Will implementation be favorable to ground/flight interface compatibility at the following 

points for any approved mission 9 

Tracking and data system (TDS)/mission design 9 -10 - +10 

TDS/spacecraft 9 -10- +10 

TDS/launch vehicle 9 -10- +10 

TDS/mission operations system 9 - 1 0 - + 1 0 

TDS/integrated systems 9 -10 - +10 

Will the implementation have a favorable interface compatibility with other TDS data types 9 -10 - +10 

Growth potential Can the design 

Serve additional purposes without modification 9 -10 - +10 

Be modified to serve additional purposes without affectmg the original purpose 9 -6 - +6 

Can the system capabilities improve with anticipated advances in technology 9 -10 - +10 

Schedule Can the design be implemented in time to meet the need date 9 0 - +3 

Safety Under fault conditions can the design be made to fail without damaging equipment or 0 - +3 


injuring personnel 9 

Technical attributes evaluation process 
Operation 


Operabihty Is the approach consistent with unattended operation 9 0 - +6 

Availability Can the system operate for an extended time period without the need for maintenance 9 0- +3 

Restorability Can the system operation be restored after failure by using DSN resources alone 9 0 - +3 

Repairabihty Can all equipment be repaired by using DSN resources alone 9 0 - +3 

Capability Will implementation increase DSN capacity to meet operational commitments 9 -10 - +10 

System validation Can the system performance be validated within the TDS 9 -10 - +10 

Can the user product be validated within the TDS 9 -10 - +10 

Technical attributes evaluation process 
Management 

Communications Will implementation be favorable to existing and planned DSN communications arrangements 9 -6 - +6 

Data Processing Will implementation be favorable to existing and planned DSN computing capabilities 9 -6 - +6 

Energy Will implementation be favorable to DSN mean energy consumption 9 -6 - +6 

Manageability Will implementation be favorable to existing management arrangements 9 -10- +10 
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Table 1 (contd) 


Attribute 

Perspective 

Score range 


Technical attributes evaluation process 



User 


Spacecraft 

Will implementation be favorable to existing spacecraft operations and the spacecraft design 
of all approved but unflown projects 9 

-10 - +10 

Mission operations 

Will implementation be favorable to MOS arrangements for all existing spacecraft and all 

-10+10 

system 

approved but unflown projects 9 


Data processing 

Will implementation be favorable to existing and planned computing arrangements 9 

-10 - +10 

Costs 

Will implementation be favorable to user costs 9 

-10 - +10 


Technical optimization evaluation process 



Implementation 


Ease of design 

Is the task within design organization experience 9 

0-+10 


Are all the functions easy to understand 9 

0 - +6 

Ease of installation 

Can installation be done without any 



Site preparation 9 

0 - +6 


Special equipment 9 

0 - +6 


Specially trained personnel 9 

0 - +6 

Ease of modification 

Are all the functions mechanized simply 9 

0 - +6 


Are the functions nearly independent 9 

0 - +6 


Are the function interfaces simple 9 

0 - +6 


Technical optimization evaluation process 



Operation 


Ease of operation 

Are the periods of operator activity short 9 

0 - +6 


Are there few periods of operator activity in a duty cycle 9 

0 - +6 

Ease of maintenance 

Is there sufficient time in the duty cycle to perform necessary maintenance 9 

0-+10 

Ease of repair 

Will failed equipment be repairable in a reasonable time period without the need for high skills, 
long experience, or special equipment 9 

0 - +6 

Life cycle cost 


Perspective 

Response 


What are the life cycle costs for each design 9 

Dollar value 


What is the DSN differential life cycle cost for each design 9 

Dollar value 


What is the annual differential cost to each user for each design 9 

Dollar value 


Are there any costs that are soft or fuzzy 9 

Yes/No 


Is the implementation favorable to DSN cost-effectiveness in terms of user hours per 
dollar of DSN funding 9 

Yes/No 
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INITIAL CANDIDATES 



SELECTION DECISION 

1 


Fig. 1 Top level block diagram of the generalized evaluation criteria process 
for selecting a design candidate 








INITIAL CANDIDATES 



SELECTION DECISION 


I 


Fig 2 Top level flow chart, NAVNET project evaluation criteria process 
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Wet tropospheric path length fluctuation will, at some level of system sensitivity, 
begin to interfere with the search for gravitational waves using the spacecraft doppler 
method This article investigates radiosonde data from Edwards Air Force Base and radio 
metric scintillation data collected over a long, nearly horizontal path in Hawaii Utilizing 
a previous hypothesis that wet tropospheric path length fluctuation is proportional to 
total wet tropospheric signal delay, the two types of data are shown to be in reasonable 
agreement for averagings times (tJ of approximately 3000 seconds The two-way 
modeled tropospheric fractional frequency fluctuation at r = 1000 seconds is 1 6 X 

IQ-14 


I. Introduction 

Proposals have recently been advanced to search for gravita- 
tional waves in ultraprecise two-way doppler data Preliminary 
estimates of gravitational wave characteristics indicate that a 
total measurement system fractional frequency fluctuation 
( ofAF/F ), where F is an S- or X-band frequency) of 1 X 10“ 1 5 
over the time scales of interest (50 to 5000 seconds) will be 
required (Ref 1) At this level, fluctuations in the wet (water 
vapor) component of the tropospheric signal delay (flj will 
surely constitute a major error source 

In a previous article (Ref 2), Berman modeled wet tropo- 
spheric fluctuation under the assumption that wet tropo- 
spheric path length fluctuation is proportional to the total wet 
tropospheric signal delay In this article, a very low frequency 
wet tropospheric path length fluctuation spectrum is con- 
structed from radiosonde data measured at Edwards Air Force 


Base, and a high frequency fluctuation spectrum is constructed 
from tropospheric radio metric scintillation data taken by 
Thompson (Ref 3) in Hawaii It is shown that the two mean 
spectra are reasonably consistent for fluctuation frequencies 
3 X 10 4 Hz under the assumption that the wet tropo- 
spheric path length fluctuation is proportional to the total wet 
tropospheric delay 

II. Radiosonde Data 

Wet zenith tropospheric signal delay values ( R wz , cm) we re 
computed from radiosonde measurements made at Edwards 
AFB during June through November 1977 and April 1978 
through March 1979 Edwards AFB is located in the Mojave 
Desert about Fifty miles north of Los Angeles and has a 
climate similar to that of the Deep Space Network (DSN) 
Goldstone tracking complex The radiosonde measurements 
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were made at irregular intervals ranging from periods of several 
hours to approximately one day Changes in zenith delay 
(A7? lvz ) as a function of separation time (r a ) were computed 
for all measurements withm 3 days of each other These were 
then sorted into “bins” of 1/ 10-day each and the mean value 
of each bin calculated Figures 1 and 2 display these mean bin 
averages (in cm) 


Whereas the radiosonde data show greater r fl = 1000 second 
(17 mmute) fluctuation that the water vapor radiometer data 
(7% versus 2%), the one-day modeled fluctuations are compar- 
able at 27% and 29%, respectively 

In the remainder of the article, the 1977 and 1978—1979 
data fits will be combined into a single model as follows 


As previously noted, Berman (Ref 2) has hypothesized that 
fluctuation in wet tropospheric delay is proportional to the 
total wet troposphenc delay To test this hypothesis, the 
changes in wet zenith delay (A R wz ) were correlated with 
mean zenith delay (7? wz ) for each of the bins in both data sets 
Computed correlation was significant in almost every case, 
with the average computed correlation coefficient ( r ) being 


IT \° 29 

’* w fr.)- 0 07«w, (lo5o) 

A fractional frequency fluctuation may be defined for this 
data 


1977 data (6 months) r = 0 59 o(AF/F) = o R ( T a )/cT a 

1978-1979 data (12 months) r = 0 26 


Based on this observed correlation, it is considered that 
the heunstic assumption of proportionality between wet tro- 
pospheric fluctuation and mean wet troposphenc delay con- 
tinues valid 


where 

c = velocity of light, cm/second 
T a = doppler averagmg time, seconds 1 * 


To utilize this relationship, each of the changes in wet 
tropospheric delay was normahzed (AR WZ /R WZ ) by the aver- 
age value of the delay Again, these were sorted into 1/10-day 
(2 4-hour) bins Figures 3 and 4 present the mean bin values 
for the normalized (fractional) fluctuations The lines in 
Figs 3 and 4 represent least squares linear curve fits to these 
data The equivalent parametric forms for these data fits are 
(j a in seconds) 

1977 data (6 months) 


The fractional frequency fluctuation for the Edwards AFB 
radiosonde data thus becomes 


o(AF/F) = 20 X 10' 14 



when a (yearly average) value of 8 cm is assumed for R 


o„ (r ) = 0 0797? 

R v a J wz 

wz 



1978-1979 data (12 months) 


III. Tropospheric Radio Metric 
Scintillation Data 

In Ref 3, Thompson, et al , describes the phase spectral 
density of wet tropospheric signal delays over a 64-km Hawai- 
ian (nearly horizontal) range The fluctuation frequency (a) 
range of vahdity for this data is 


a* (t ) = 0 069R wz 

WZ 



3 X 10' 3 Hz < p < 3 X 10' 1 Hz 


These can be compared to the estimate Ref 2 made based on a 
very prehminary examination of tropospheric fluctuation 
obtained from water vapor radiometer data 


To compare these data to the Edwards AFB data, one must 
scale the results by the appropriate total wet delays The exact 
troposphenc water vapor distribution at the time of the 
Thompson experiment over the 64-km Hawau radiometer 


(r a ) = 0 027? 

wz 


wz 



1 For the remainder of this article, “measurement separation time” will 

be equated with “doppler averaging time ” 
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range is unknown Using the U S Standard Tropical Atmos- 
phere, 15° N (Ref 4), it is estimated that there were a total of 
87 gm/cm 2 precipitable water along the 64-km path A stan- 
dard year-average zenith atmosphere at Edwards AFB has 
1 3 gm/cm 2 of precipitable water This ratio of nominal delays 
(proportional to integrated water vapor) is 66 9 and therefore 
the phase spectral density of the Thompson data must be 
scaled downward by (66 9) 2 Thus, the Edwards AFB equiva- 
lent zenith troposphere X-band phase spectral density (based 
on Hawaii measurements) is 

P 0 (u) = 1 09 X 1CT 7 jT 2 57 rad 2 Hz’ 1 


IV. Comparison of Radiosonde and Radio 
Metric Scintillation Phase Spectral 
Density 

By equating the radiosonde fractional frequency fluctua- 
tion (squared) with the previous expression for Allan variance, 
a value of A may be deduced Using this procedure, a phase 
spectral density for Edwards AFB radiosonde data is obtained 

P 0 (v) = 6 76X KT'V 58 rad 2 Hz -1 
over the range of validity 3 X 10 -6 < v < 1CT 4 


3 X 1(T 3 Hz<v<3X 1CT 1 Hz 

Reference 5 provides the following expression (with the Ref 6 
correction) for Allan variance ( a y 2 ) derived from phase spec- 
tral density (v Q is transmission frequency) 


a 2 (r ) = 4ir 2 v 2 At a 3 

y v a’ 0 a 



sin 4 (7r Z) dZ 


Figure 6 shows the phase power spectra for both the Hawaii 
microwave data and the radiosonde data taken at Edwards 
AFB The microwave data have been normalized to year- 
average zenith values for Edwards AFB, thus, both curves 
represent phase spectra for the troposphere in temperate semi- 
and regions As stated in Section III, the spectrum at very low 
frequencies flattens out due to the decrease in long-term 
troposphere variation for periods greater than a few days 
Therefore, the “combination” of the model segments as in 
Figs 5 and 6 is intuitively agreeable 


where A and a are defined from the expression for phase 
spectral density P^iy) =Av~ a (as previously given) 


Hence from the above one obtains from the Thompson data 
as corrected to Edwards AFB, for X-band frequency ( v 0 = 
8 4 GHz) and a nominal 8-cm delay 


a (t ) = 1 15 X 1(T 14 

y v a ' 



V. Tropospheric Model 

Following Ref 2, a factor of (sm 0 )' 1 is added to account 
for non-zenith delays, where 6 = elevation angle One then has 
a “combined” radiosonde/radio metric scintillation model for 
tropospheric fluctuation (cm) 


/ r \° 29 

* 0 074(m»r‘/!». (io5o) ■ 


r > 3000 s 

a 


The radiosonde fractional frequency fluctuation and radio 
metric scintillation Allan variance data 2 are plotted in Fig 5 
As can be seen, they are reasonably consistent In fact, the 
steeper slope of the very low frequency (Edwards AFB) data is 
reasonable, since one would expect the tropospheric fluctua- 
tion spectrum to begin to flatten at about 3 days (4 X 10“ 6 
Hz) due to movement of large-scale air masses 


2 For the purpose of this article, the radiosonde fractional frequency 

fluctuation is considered to be equal to the square root of the Allan 

variance 


a ( d,R ,t ) = 0 043 (sin 6) l R 

R v 9 wz 9 a 7 v / wz 

w 



r < 3000 s 

a 


Edwards AFB is a reasonable representative of a temperate, 
semiand climate, as are the Deep Space Network stations 
Table 1 gives the seasonal variation of the mean wet zenith 
tropospheric signal delay during 1977 through 1979 From 
this data, a mean yearly value of 8 cm was computed, with a 
mean summer extreme of 1 6 cm and a mean winter extreme of 
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4 cm The combined radiosonde/radio metric scintillation 
model for these conditions is presented in Fig 7 


VI. Discussion 

Previous work (Ref 2) modeled wet tropospheric path 
length fluctuation under the heuristic hypothesis of propor- 
tionality between wet tropospheric path length fluctuation 
and total wet tropospheric signal delay In this article, radio- 
sonde measurements collected over 1 5 years at Edwards AFB 
are used to demonstrate the soundness of this hypothesis 
Utilizing the hypothesis, wet zenith tropospheric fractional 
frequency fluctuation computed from the Edwards AFB radio- 
sonde data and the Hawauan radio metric scintillation Allan 


variance data are shown to be reasonably consistent in the 
region of r a ~ 3000 seconds The fractional frequency fluctua- 
tion at T a = 1000 seconds is approximately 1 1 X 10~ 14 , or 
about a factor of two higher than the value estimated in Ref 2 
from water vapor radiometer data Since this is a “one-way” 
measurement, a factor of \/2 is applied to obtain the equiva- 
lent “two-way” value of 1 6 X 10“ 14 Since a value this large 
will certainly impact any attempts to search for gravitational 
waves using the spacecraft doppler method, there exists a clear 
need to measure wet tropospheric path length fluctuation at r a 
= 1000 seconds, instead of relying on extrapolations mto this 
critical fluctuation region The water vapor radiometer is sug- 
gested as a possible instrument for measuring tropospheric 
path fluctuation at r a = 1000 seconds, and perhaps ultimately, 
providing wet tropospheric cahbration as well 


References 


1 Berman, A L , “The Gravitational Wave Detection Experiment Description and 
Anticipated Requirements,” in The Deep Space Network Progress Report 42-46, pp 
100-108, Jet Propulsion Laboratory, Pasadena, Cahf , Aug 15, 1978 

2 Berman, A L , “Parametric Modeling of Low-Frequency Water-Vapor-Induced Tropo- 
spheric Path Length Fluctuations” in The Deep Space Network Progress Report 42-47, 
pp 72-76, Jet Propulsion Laboratory, Pasadena, Cahf , Oct 15, 1978 

3 Thompson, M C , et al , “Phase and Amplitude Scintillations in the 10 to 40 GHz 
Band,” IEEE Transactions on Antennas and Propagation, Vol AP-23, No 6, 
Nov 1975 

4 Valley, S L , Editor, Handbook of Geophysics and Space Environments, 1965 Edi- 
tion, p 2-10, McGraw Hill Book Company, Inc , New York, 1965 

5 Armstrong, J W , Woo, R , and Estabrook, F B , “Interplanetary Phase Scintillation 
and the Search for Very Low Frequency Gravitational Radiation,” in The Astrophysical 
Journal, Volume 230, June 1, 1979 

6 Armstrong, J W , Estabrook, F B , and Woo, R “Corrections to Allan Variances ,” 
Astrophysical Journal (submitted) 


82 



Table 1 Mean zenith wet tropospheric signal delay (R wz ) as a 
function of season during 1977-1979 


Penod 

Delay, cm 

January-March 

6 

April- June 

8 

July-September 

11 

October-December 

7 

Yearly Average 

8 

Mean Summer Extreme 

16 

Mean Winter Extreme 

4 
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Fig 1 Zenith wet tropospheric path length change versus 
averaging time, 1977 




0 05 01 02 03040507 10 2030 

T q , days 


Fig 2 Zenith wet tropospheric path length change versus 
averaging time, 1978-1979 
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Fig 3 Zenith wet tropospheric fractional path length change 
versus averaging time, 1977 




Fig 4 Zenith wet tropospheric fractional path length change 
versus averaging time, 1978-1979 
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5 Zenith wet tropospheric fractional frequency fluctuation 
versus averaging time 




Fig. 7 Seasonal variations in combined radiosonde/radio metric 
scintillation model 


Fig 6 Phase spectral density computed from radiosonde and 
radio metric data 
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GCF-NOCC Reconfiguration 

J P McClure 

DSN Data Systems Section 


The equipment and computer programs in the Network Operations Control Center 
(NOCC) and the Ground Communications Facility’s ( GCF’s ) Central Communications 
Terminal at JPL are being rearranged and supplemented to provide an improved opera- 
tional capability 


The computer portion of the GCF’s Central Communica- 
tions Terminal (CCT) in the basement of the Space Flight 
Operations Facility was designed and implemented separate 
from the Network Data Processing (NDP) portion of the 
NOCC These independent designs, implemented at different 
times, used more computers and programs than an integrated 
design would have required The desire to reduce the computer 
count and to provide an improved operational capability led 
to the GCF-NOCC reconfiguration effort 

This effort, when completed in 1981, will 

(1) Require fewer minicomputers (16 versus 20) 

(2) Require fewer computer programs (6 versus 8) 

(3) Significantly reduce magnetic tape handling 

(4) Provide a much-improved monitor and control 
capability 

(5) Reduce operator requirements and provide the base 
for a two-operator CCT 

The present configuration of the CCT and elements of the 
NDP are shown in Fig 1 High-speed data (HSD) are accepted 
by the Error Detection and Correction (EDC) processors, 
which operate interactively with the DSSs to correct HSD 


transmission errors The corrected HSD are fed to the High- 
Speed Switch computer, which routes the data, as appropriate, 
to the Mission Operations Centers (MOCs), the Network Log 
Processor (NLP) and (not shown) the Remote Mission Opera- 
tions Centers (Ames and Germany) The NLP records both 
high-speed and wideband data (WBD) on magnetic tape to 
form a Network Data Log (NDL) The data are also forwarded 
to the Network Communications Equipment computers in 
Building 202, which routes the information to the proper NDP 
computer for analysis 

The Data Record Processors (DRPs) recall missing data 
from the DSSs (in response to data gaps sensed by the NDP 
machines) and record them on magnetic tape to generate a 
“fill” tape An off-line DRP merges the NDL and fill tapes to 
produce an intermediate data record (IDR) The IDR is the 
prime deliverable product of the DSN 

The Central Communications Monitor (CCM) senses and 
displays the performance of the entire GCF, receiving reports 
from DSS and CCT equipments 

Display information generated by the NDP’s Real-Time 
Monitor (RTM) computers is forwarded via the NLP to the 
Digital Display Processor (DDP), thence to the Video Assembly 
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Processor (VAP) for display in the Network Operations 
Control Center in Building 230 

The reconfigured GCF-NOCC, shown in Fig 2, consider- 
ably simplifies the data paths and processing Both high-speed 
and wideband data are routed through automated circuit 
switches and then accepted by the Error Correction and 
Switching (ECS) computers The ECSs error correct the 
HSD, optionally record it on a Front End Record (FER) mag- 
netic tape, and route it, as appropriate, to 

(1) The MCCC 

(2) Remote Mission Operations Centers 

(3) The Data Records Generation (DRG) computers 

The DRG records HSD on disc, senses gaps, recalls the 
missing data from the DSS, and then merges the real-time and 


recalled data as it writes the IDR on magnetic tape Wideband 
data are recorded in IDR format directly on magnetic tape If 
recalls are needed, the DRG also merges the real-time and re- 
called data to form a complete WB IDR 

The CCM, NCE, DDP, and VAP functions are largely 
unchanged, though the NCE and the DDP communicate 
directly over their own interconnecting circuits, as do the 
ECSs and NCEs 

The reconfigured GCF-NOCC capability is scheduled to be 
completed in March 1981 Subsystem-level testing begins in 
May 1980, accompanied by a phased rearrangement of the 
hardware Training and operational testing start in November 
1980 and continue until the full capability is placed on-line 
in March 1981 
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Preliminary Telemetry Operations Experience With the 
Real-Time Combiner: 1 November 1978 to 
1 November 1979 

N K Simon and C Hoynes 
Deep Space Network Support Section 


A prototype version of a two station, real-time signal combiner was installed at Gold- 
stone DSS 14/(12) and used to support the Voyager Mission at Jupiter and the Pioneer 11 
Mission at Saturn This article reports on early combined station operating experience 
with a Real-Time Combiner (RTC) in the loop 


I. Introduction 

On September 21 , 1974, the two 26-m antennas at DSSs 12 
and 13 were arrayed with the 64-m antenna at DSS 14 to 
improve the telemetry data received during the Manner 73 
Mercury encounter (Ref 1) The experimenters claimed an 
improvement in SNR of 0 7 dB over the operation of DSS 14 
alone The arraying was accomplished with a signal combiner 
that took baseband signals and phase correlated them to pro- 
duce a single additive signal The combined output signal was 
then input to a DSS 14 receiver and demodulator “for further 
processing ” The experiment was a laboratory research effort 
meant for one time use during the MVM 73 encounter It also 
demonstrated the value of arraying the antennas at a DSS 
complex A reproduced graph of BER vs time from the Sep- 
tember 21 encounter is shown in Fig 1 

The current effort began m FY78 The objective was to 
develop a permanent Real-Time Combiner (RTC) installation 
for the 64-m/34-m net A prototype unit became available in 
November 1978 and was successfully used for the Voyager 1 
and 2 Jupiter encounters and for the Pioneer 1 1 Saturn 


encounter It is the operation with the prototype unit that 
is the subject of this report 

The basic system block diagram for RTC operation is 
illustrated in Fig 2 Baseband signals from a 64-m and 34-m 
station are fed into the combmer, phase locked, and sent 
through a standard telemetry string for normal processing 
The SNR gam from combining over a 64-m station alone is 
dependent on the actual SNR spread between the two sta- 
tions For a nominal spread of 5-dB between the 64-m and 
34-m antenna the RTC theoretical gam is 1 2 dB This com- 
pares with 0 8 dB available from adding two 26-m stations to 
the 64-m signal 


II. Test Program 

A. CTA21 

The prototype RTC unit was set up in CTA 21 and initial 
testing was conducted in November/December 1978 
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Two methods of simulating a 2-station configuration at 
CTA 21 for testing RTC performance in an operational system 
environment were considered Method 1 Receiver mput ter- 
minations were isolated with a common modulated RF test 
signal fed to both receivers Levels were adjusted by means 
of S-band attenuators in each hne Method 2 Separate RF 
signals modulated by a common data-modulated subcarrier 
from the Simulation Conversion Assembly (SCA) were fed to 
each Receiver 

In both cases the baseband outputs of each receiver were 
fed into separate RTC channels Results indicated both 
methods to be acceptable Method 2 was preferred for ease 
of configuration The actual configuration used for CTA 21 
is shown in Fig 3 

With 10 dB difference between input signals the combiner 
provided a 0 4-0 5 dB SNR gam over that of the strong 
signal alone This is in the correct range for what would be 
predicted by theory 

Although the initial test configuration tried was valid, 
several tests over a period of weeks were required to obtain 
the above result This was because the following difficulties 
thwarted the experimenters during early tests 

(1) Laboratory power variations and noise interfered with 
the signal levels during the tests at CTA 21 It was 
found that the system appeared to be more stable 
after 5 pan Therefore, most testing was shifted from 
daytime to evenmgs and weekends 

(2) SNR values computed by the Telemetry Processor 
Assembly (TP A) fluctuated ±1 5 dB In trying to mea- 
sure a 0 5 dB gam, an SNR value good to ±0 1 dB is 
needed It took some time before a technique was 
evolved to do this with any confidence in the results 
The problem was exacerbated by the power varia- 
tions noted above 


provision for varying the delay in the symbol error line from 
the SCA to the Symbol Synchromzer Assembly (SSA)) Two 
alternatives remained 1) use a less desirable square wave bit 
pattern or 2) lengthen test time and use the estimated SNR 
mean for performance measurements Since combiner perfor- 
mance was measured by serially testing and comparing com- 
bined and uncombmed performance, and with system stability 
uncertain, most of the testmg at CTA 21 was performed with 
square-wave data using the symbol error count as the prime 
measurement data 

Note that the purpose of testing at CTA 21 was to mea- 
sure performance of the RTC only, therefore, all measure- 
ments were relative (combined to uncombmed), masking 
degradation of the remaming system components 

Of primary concern durmg this period was the source of 
the SNR fluctuations 

(1) Were the SNR fluctuations real 9 

(a) How accurate was the estimator 9 

(b) Was it averaged over sufficient time 9 

(2) Were there system problems causing the fluctuations 9 

(a) If so which component in the telemetry string 
varied 9 

(3) Was the signal varying 9 

Solutions included, 

(1) Running tests on swing or midshift 

(2) Switching hardware to improve the telemetry string 

(3) Running extensive BER (bit error rate) tests on the 
data to determine the appropriate data collection 
time and isolate causes of SNR variation 

It was found that approximately 10 minutes of data at or 
above 2 kb/s 1 was needed to obtain stable SNR results BER 


A Block IV receiver degradation problem also existed, and 
was eventually determined to be in the telemetry phase 
detector module Several attempts were made to improve its 
performance, with none of them successful So that testing 
could proceed, it was decided to use the Block III Receivers, 
and all successful testing at CTA 21 was performed with 
these receivers at S-band frequency 

Because of the variance in the telemetry processor esti- 
mated SNR outputs, a symbol error count was the more 
desirable criterion for measuring system performance How- 
ever, because of the delay through the combiner, a symbol 
error count at the higher data rates could not be made usmg 
the more desirable psuedo noise pattern (The system has no 


'More specifically @1 dB ± 0 1 dB = SNR 
SER = 0 056 ± 0 0025 


n 


SER(1 - SER) 

2 

a 


(0 056) (0 944) 
(0 00125) 2 


n = 350,000 bits 


3500 bits 
2000 bps 


17 5 s(@2 kb/s) 


Therefore at 2 kb/s it should only take approximately 18 s to establish 
the SNR to within 0 1 dB But with all the uncertainties of telemetry 
strmg operation, a 10 minute test is recommended to assure system 
stability and to validate proper operation of the stnng 
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tests using controlled data verified that SNR data from the 
TPA was accurate and stable over this time period It is essen- 
tial that all equipment in the telemetry string be calibrated 
and in specification, to obtain reliable data 

B. Goldstone 

The RTC was then moved to Goldstone and installed for 
a first test with live spacecraft data on December 19, 1978 
The first two tests were inconclusive A test was run on 
February 21, 1979 that showed a combined signal gam of 
from 0 4 to 1 0 dB This was still inconclusive proof of com- 
biner performance since the expected value was a steady 1 2 dB 
gain However, there were several hours of combiner operation 
with positive gain registered This was a step in the right 
direction Because station availability was extremely limited 
it was necessary to perform almost all testing at Goldstone 
durmg regularly scheduled Voyager passes on a noninterfer- 
ence basis This was highly undesirable in that there was no 
control over data rates or signal level and available configura- 
tions were limited This, coupled with several unknowns (such 
as proper telemetry receiver phasing requirements, antenna 
gain at DSS 12, CONSCAN stability, etc) and uncertainties 
(such as best procedures for configuring an array system, 
total delay between stations, Block IV receiver performance, 
etc ) was probably the greatest contributing factor to the long 
time period required to uncover many of the arraying prob- 
lems, the most significant of which are discussed below 

(1) If DSS 12 dropped lock, the only evidence at DSS 14 
was that the RTC went from 1 dB gam to 1 dB net 
loss This relates to how visible the performance of 
the RTC “system” is to station operators during 
arraying operations The subject is being studied by a 
working group 

(2) The micro-wave system could invert the signal from 
DSS 1 2 This simply required a procedural change 

(3) The baseband module in the Block IV receiver cost 
1 dB in SNR over that observed in the Block III 
receiver The baseband module problem is the subject 
of an approved design change which will be completed 
in the 64 Meter Network on 1 June 1980 As an 
interim solution the Block IV IF was cabled through 
the Block III baseband modules as shown in Fig 4 

(4) The antenna spread between DSS 12 and 14 was 5 to 
6 dB at X-band instead of the 4 2 dB expected, due to 
lower than expected antenna gain at DSS 12 DSS 12 
gain was out of specification when the antenna was 
measured and post encounter adjustments have not 
unproved it 


(5) The receiver signal level at DSS 12 would vary ±1 5 dB 
in a gusty wind (20 knots gusting to 35 knots) The 
Automatic Gain Control (AGC) variation during gusty 
winds probably was resolved by antenna servo modifi- 
cations made to reduce jitter Certainly AGC from 
both antennas needs to be monitored to track arrayed 
system performance CONSCAN at DSS 12 was first 
suspected when this condition occurred Investigation 
proved otherwise However, CONSCAN does have 
some glitches with the current software and this is 
planned to be corrected 

(6) No means existed to stop and restart the combiner 
without complete re-initialization Initialization of the 
combiner required a minimum of eight operator inputs 
and is time consuming The production version of the 
combiner will allow restart with one input and con- 
figuration change without restart Configuration and 
prepass calibration procedures required two hours 
which was seldom available with the prototype system 
The Saturn configuration should require no more 
than 30 minutes if string selection can be eliminated 
(See discussion, Section VI ) 

III. Voyager/Jupiter Encounters 

The final array configuration at Goldstone used for all 
encounters thus far is shown in Fig 4 The prototype RTC 
was used for 10 days around 5 March 1979, for the Voyager I 
Jupiter encounter, and for 14 days around 2 July 1979, for 
the Voyager II Jupiter encounter RTC gam was determined to 
be 0 8 to 1 0 dB and the combmed telemetry data stream was 
used by the project for most of each track Combiner up-time 
was > 95 percent The received carrier frequency was X-band 
and the SNR difference between DSS 14 and DSS 12 was 
approximately 5 dB 

IV. Pioneer 11 /Saturn Encounter 

A real bonus was achieved when the array system was used 
at Goldstone for 10 days around the Pioneer 11 encounter of 
Saturn on 1 September 1979 Signal levels (at S-band) were 
too low to use the RTC phase lock loop mode to combine the 
signals from DSS 12 and 14 DSS 14 is able to improve its 
gam by 3 dB at S-band by going into a “listen only” mode 
DSS 12 does not have the listen only option so that its system 
noise temperature, SNT, remains relatively high 2 For this 
reason and due to other equipment differences the station 
SNR difference goes from 5 dB at X-band to 9 dB at S-band 
With this spread the theoretical combiner gain is only 0 5 dB 


2 at DSS 14, SNT nom = 16 K 
at DSS 12, SNT nom = 26 K 
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The RTC was operated in the “dead-reckoning” mode 
whereby the signals are phased by calculating the delay time 
between the two stations using the combiner minicomputer 
Normally this calculation is made to bring the signals in range 
of the phase lock loop (PLL) In the dead-reckoning mode the 
PLL is bypassed 

RTC gain during the Pioneer-Saturn Encounter was 0 3- 
0 4 dB This was sufficient to permit the project many more 
hours of 1024b/s data at Goldstone and substantially enhanced 
the imaging mission With low RTC gain, system performance 
was difficult to measure in real time The best potential mea- 
sure of system performance with Pioneer telemetry is frame 
deletion rate of the sequential decoder Unfortunately, this 
parameter was calculated based on a running average It proved 
almost valueless in establishing short term performance of the 
telemetry string The real measure of performance was Gold- 
stone’s ability to receive 1024 b/s data while the rest of the 
network functioned at 5 12 b/s 


V. Future Use of Arrayed Systems 

Permanent RTC installations are to be made at the 64-m 
net in the first and second quarter of 1980 This will provide 
the capability to array 64-m/34-m stations at Goldstone, 
Madrid, and Canberra in time for the Saturn encounters by 
Voyager in November 1980 and August 1981 Some improve- 
ments in display capability will be implemented by that time 
to permit combined performance evaluation from a central 
location, either at the 64-m station or at JPL Ninety days of 
continued RTC operation is currently being planned for the 
first encounter and 30 days for the second encounter 

The most striking operational system problem facing the 
test investigators at the start of Goldstone array operations 
was lack of visibility into the arrayed system performance 
Voice communication provided the most useful means of 
station status monitoring and very often an investigator had to 
make trips from DSS 14 to DSS 12 and back to solve prob- 


lems No information on the microwave system or DSS 12 
status was available at DSS 14 where the combiner was in- 
stalled (see Fig 4) This lack of visibility led to the formation 
of a combiner working group at JPL This group is concen- 
trating on improving combmer operational visibility for the 
Voyager/Saturn encounters 

The RTC produces a gain of approximately 1 dB when a 
34-m signal is added to a 64-m signal in standard X-band 
configuration In order to measure performance it is necessary 
to know the gam to, say, 0 1 dB This says that the entire 
telemetry strmg must be known to ±0 1 dB Furthermore, 
each component in that string must be calibrated and within 
spec to 0 1 dB in order to make combiner gain meaningful 
During testmg, investigators found several instances of 1 dB 
variations in other components of the telemetry string (refer 
to earlier discussion of high winds, the Blk IV baseband mod- 
ule, and DSS 12 antenna gam) The message is that a care- 
fully calibrated, tuned, and stable DSS is a prerequisite to 
successful combmer operations During the encounters, par- 
allel telemetry strings were performance checked daily to 
determine optimum DSS configuration This was time- 
consuming and would be unnecessary with proper calibration 
SPTs run on telemetry strmgs bi-weekly should ensure good 
station calibration 


VI. Future System Design 

The success of increasingly complex arraying operations 
envisioned in Network Consolidation will depend on sufficient 
automation in areas of system calibration, array configuration, 
string selection, and display formatting In particular it would 
be desirable to have the combmer analyze and disconnect 
antenna receiver systems that did not meet performance specs 
As was shown earlier, an out of lock receiver can provide a 
significant negative gam to combmer system performance 
These future array systems will be sufficiently complex that 
operators will not be able to react in real time to system 
changes 


93 



Acknowledgment 


The combiner was developed by Larry D Howard under the direction of R Bruce 
Crow of Section 333 The authors were privileged, as members of DSN Operations Sec- 
tion 377, to become involved in combiner operations during development and early 
systems testing 


Reference 


1 R A Winkelstein, “Analysis of the Signal Combiner for Multiple Antenna Arraying”, 
DSN Progress Report 42-26, January and February 1975, pp 102-1 18, 15 April 1975 


94 















Fig 3 CTA 21 test configuration 



* MODIFIED FOR X-BAND 

Fig 4 Encounter configuration, real-time combining system block diagram 
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Viterbi Decoding Modified for Sources With Memory 
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This article investigates the gain in signal-to-noise ratio that may be realized by using 
redundancies inherent m TV data to modify the decoding metric of the DSN Viterbi 
decoders This modification would take into account the memory in a TV scan line to 
change the transition probabilities as originally computed using independent data The 
resulting data rate gam or error probability decrease is achieved without requiring any 
spacecraft modifications or additions A preliminary examination of this concept on a 
binary symmetric channel rather than a Gaussian channel and using a simplified Markov 
source model involving two-level or hard-clipped TV shows that very substantial decreases 
in error probability may be achieved 


I. Introduction 

Viterbi decoding for the case where binary messages from a 
memoryless source are transmitted over a memoryless channel 
implements maximum likelihood (ML) decoding In this case, 
this is equivalent to maximum a posteriori probability (MAP) 
decoding since all messages are equally likely However, when 
the source has memory, such as for video, ML decoding, and 
hence the Viterbi decoding algorithm, as it is usually imple- 
mented, is not optimal For such cases work has been done 
(see, e g , Ref 1) on encoding the source before transmission 
and then using a source decoder at the receiver This method 
utilizes the channel more efficiently However, in cases where 
the transmitter complexity needs to be minimized, such as 
when it is on a spacecraft, it is sometimes undesirable to have a 
source encoder, which may be a complex piece of equipment 
In such a case, it could be worthwhile modifying the Viterbi 
decoder at the receiving end so as to increase the transmission 
rate, while leaving the transmitter untouched 


In the present work, we consider the simple case of a 
Markov source of arbitrary order S, and a memoryless channel, 
and show how to modify the metric in the Viterbi decoder to 
implement MAP decoding Calculations of bit error probability 
(BEP) for a decoder using this new metric become message- 
dependent An upper bound on the BEP can be obtained for 
any given message by a method to be described, but elegant 
bounds similar to those in conventional Viterbi decoding (see, 
e g , Refs 2 and 3) are only possible for a few messages, and 
these are derived An upper bound on overall BEP has also 
been similarly derived for low-entropy sources For a specific 
K = 3 constraint length code, computer simulations of the 
decoder for a first-order Markov source and a binary symmet- 
ric channel (BSC) indicate that the BEP with the new metric 
are lower than the BEP with the conventional metric, as 
predicted, with the improvement being more noticeable at 
rates well below channel capacity However, many more runs 
of these simulation programs are needed to determine more 
exactly the extent of the improvement for various parameters 
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Also, we are interested in the Gaussian channel and not the 
binary symmetric channel for which the simulations in this 
work have been done 

It has been observed that the Viterbi algorithm (VA), in its 
most general form, is “a solution to the problem of MAP 
estimation of the state sequence of a finite-state discrete-time 
Markov process observed in memoryless noise” (Ref 4) But 
the Viterbi Algorithm does not seem to have been used as a 
MAP estimator in the context of decoding convolutionally 
encoded messages from a Markov source It has, however, been 
used in the problem of text recognition (Refs 5 and 6) Also, 
error bounds like the ones given here, for cases where the 
source has memory, do not seem to have been obtained, nor 
are there available comparisons between the MAP and ML 
estimator applications of the VA to decoding problems for 
sources with memory 

II. The New Metric 

A. Definition for a BSC 

Consider a BSC with channel error or transition probability 
p e and a convolutional code of constraint length K, where b 
bits are transferred at a time into the encoder and n symbols 
are output at a time (l e , a rate b/n code) Here, “at a time” 
means at each clock pulse There are bK bits in the encoder 
register 

A Markov source, which we define below, is a simple 
example of a source with memory A Markov source of order S 
is defined as one having the property that the probability 
distribution of the « th bit it outputs (for n> S) depends only 
on the values of the previous S bits output That is, if the 
sequence of bit output by the source is {u 1 ,u 2 , , w„, }, 

then for n >S, the conditional probability /?(w„ I i , , «i) 

can be expressed as 

p(u | u ,, ,u,)=p{u I u ,, ,n „) (la) 

n n- 1 ’ ’ 1 J n n - 1 n-S ' v ' 

A symmetric first-order Markov source has 

P(1|0) = p(0|l) -p s (lb) 

and 

p(0|0) = p(l|l) = < 7 s = 1 - p s (lc) 

Here p s is called the source transition probability In this 
paper, whenever we refer to a first-order Markov source, we 
mean this symmetric one 


We first consider a Markov source of order S KbK - 1 , and, 
for this case, we show that MAP decoding is achieved by using 
the Viterbi decoding algorithm with modified branch metrics 
M’ ml defined below Each branch defines a transition out of a 
state s involving b{K - 1) bits to a new state t uniquely defined 
by the b incoming bits and the b(K - 1) bits of state s Thus, if 
the Markov source memory extends over fewer than bK previ- 
ous bits, there is a unique a prion probability p ml of occur- 
rence of each branch in the state diagram given by Eq (2b) 
below 

If S> bK - 1 , we can follow exactly the same decoding 
procedure with the metrics which we define below, with the 
exception that the state diagram and the trellis diagram for the 
decoder now have 2( s- * +1 ) states, with each state represented 
by (S - b + 1) bits rather than 2 b ^ K ~ 1 ^ states, each repre- 
sented by b(K - 1) bits (Each node will, as usual, have 2 b 
branches merging at it and 2* branches emanating from it ) 
However, since S- b + 1 > b(K - 1 ) for S > bK - 1 , this 
increases the decoder complexity So, if b is small, we might as 
well increase K , which increases the encoder complexity only 
slightly, but reduces BEP considerably Thus, it may often be 
worthwhile choosing b , K , S such that S < bK - 1 holds, if 
some other constraints are not violated in doing so If b is 
large, however, increasing A even by 1 increases the number of 
states by a factor of 2 b This may increase the decoder 
complexity much more than if we were to keep K fixed and 
have 2< s ~ 6+ 1 > states in the decoder to take care of the source 
memory S being greater than bK - 1 

The branch metric of conventional Viterbi decoding for the 
r ,h branch of the m th possible path through the trellis is ~d ml 
where d is the Hamming distance between the n-dimensional 
code subvector for the branch and the corresponding n-dimen- 
sional received subvector (see Ref 2) We define the new 
metric as 


m ‘ ml (1-Pj 

In - 

Pe 


p ml can be written as the conditional probability 

p = p(u , u . „.) (2b) 

i mi m(i-l) m(i-Sy ' 

for an S* 11 order Markov source, where w , is the 1 th bit of the 

th m 1 

m possible message of block length B Note that m runs 
from 0 to 2 s - 1 We assume throughout that B » K, and 
that the last b(K - 1) bits transmitted are always 0 
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B. Definition for an Arbitrary Memoryless Channel 

The conventional branch metric , for the I th branch of 

m i 

the 777 th path is the logarithm of the conditional probability of 
the / th received subvector (Ref 2) 

M = In p(y\x ) (3) 

where x mt is the «-dimensional code subvector of the 777 th 
message sequence for the 7 th branching level, and ~y t is the 
corresponding n-dimensional noisy received vector The modi- 
fied metric in this case is 

M' = lnpCF|5c ) + lnp (4) 

mi r v ^ / mr m i 7 

where p is defined in Eq (2b) 

C. Derivation of New Metric 

We now derive the new metric defined above for the case 
where the source memory is S 


where p(m) is given by Eq (5), x ml ,~y l have the same meaning 
as in Eq (4), and p(y) is the probability of receivings This 
probability is constant for all messages m , and hence can be 
dropped from the maximization, or, equivalently, we can 
maximize the log of this probability (after dropping p(y)), l e , 

B 

,n P 'ap = ln I I + ln ( 7 ) 

1=1 

= £ {lnp(7jy mj ) + lnp mi | (8) 

i=i 

From this it follows that the branch metric for the 7 th branch 
of the 777 th message path is given by Eq (4) For a BSC with 
transition probability p e , we can further simplify Eq (4) using 
the fact that (Ref 2) 

pO 7 ,!^,) = pf m ‘ (1 -P e ) dmi (9) 


Let i , ’ “mj) be the 777 th possible message sequence 

of block length B Then the a priori probability of the se- 
quence (77 ml , ,u mB } is 

p(m)=p{u mV ,u mB ) 


Thus, Eq (4) becomes 


M = 

mi 


- d ln 
m 1 


+ 77 In (1 - P e ) + lnp , 


; (io) 


P( U mB^ U m(B-iy > U m (B-S)^P^ U m (B- 1 ) ^ U m (fl-2 ), 
U m(B-S- 1)) P( U m\\ U mO ,U m{r\y ,U m(\-S)^ 

PrnBPmB-l P m l ( 5a > 


Or, dropping the 77 ln (1 - p e ), which is constant for all 
branches, and dividing by 



where p m is defined as in Eq (2b), and we define 


we get 


{77 ,« . , ,u „,}={(), 0, ,0} (5b)* 

Then the MAP decoder should choose that f?-bit block corres- 
ponding to message m that has the maximum a posteriori 
probability given by 


P 


ap 


fl P( y . [X m) 


1=1 


{p(m)} 

p(y) 


( 6 ) 


M' = 

mi 


In p 

-d + m ‘ 

mi 


In 


i -p. 


which establishes Eq (2a) 


III. Error Probability Bounds With 
New Metric 


(ii) 


* Alternatively, we may take care of the initial bits by defining p ml , , We now derive upper bounds on the probability of error 

p mS each to be 0 5 with the new metric in a manner paralleling the derivation of 
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Eq 4 4-8 of Ref 2, which is repeated as Eq R1 below or 


E[n b m < X X *( d >oz d (ri) 

1 =\ 

where E[n b {j)\ is the expected number of bit errors caused by 
an incorrect path diverging at node /, a(d,i ) is the number of 
paths diverging from the all-zeros path (at node /) at distance d 
and with i “l’s” in its data sequence over the unmerged 
segment, d^ is the minimum distance of any path from the 
correct one, called the free distance, and Z is defined in Eq 
(R2) 


z = X VpiOOpoM ( R2a ) 

y 


Here p x (y), P Q (v) are the probabilities of receiving a given 
value of a bit y when the corresponding transmitted bit is a 1 
or a 0, respectively For the special case of a BSC, Z reduces to 


Z = sjAp e (\-y c ) (R2b) 

A Extension of the Bhattacharya Bound 

We first extend the Bhattacharya bound to the case in 
which MAP detection occurs rather than ML detection 


A . = T. f(y)p(y\x ) 


(13) 


where 


Now 


f(J) 


1, forvEA 


0, otherwise 


(14) 


m = 


o p(y\x ) 

1 < / . , forvEA 

p m p{y \x m ) ■ mm 


r p(y\x’ ) 

0 < / , otherwise 

p,„ p(y\x) 


(15) 


where the first part of the inequality follows from the defini- 
tion of A mm > and the second is trivial 

Thus, for ally, 


m < 


1 p(v \x' 

m r v - m ’ 

p m p(y\x m ) 


(16) 


If p E (m -> m') denotes the probability that message m is 
decoded when message m is sent and only two alternatives (m 
and m) exist, and if 3c m is the signal vector sent and y the 
one received, and if x' m is the other possible signal vector, then 
the Bhattacharya bound states (Eq 2 3 15 of Ref 2) that 


Using this in Eq (13) gives 


P E ( m ->• m') < XI ( 17 ) 

y V 


p E (m ^m') < X] P^y\x' m ) p(y\x m ) (R3) 
y 

We now extend this bound to the case of MAP detection We 
have 


Thus, the only difference between this inequality and the 
corresponding one for ML decoding (Eq R3), i e , the usual 
Bhattacharya bound, is the extra factor V p' m lp ln in Eq (17) 
This bound is tighter than the usual Bhattacharya bound 
whenever the decoded message m is less probable than the one 
sent, m, i e , if 


P t (m 


m) 


X p<y\x m ) 

ye A 


(12a) 




(17)' 


where 


( P m .p(y\x' m ) ) 


B. MAP Decoding Error Bound 

This bound is derived using the union bound and the 
extended Bhattacharya bound (Eq 17) The fact that the a 

(12b) 

pnon probabilities of the various possible messages are un- 
equal changes the locations of the MAP receiver decision 
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boundaries in signal space but the expression for the union 
bound is unchanged Thus, we get, for the message m actually 
being sent, that the total probability of error considering all 
other possible messages is 


P Em < E ( I8 ) 

m'=tm 

< E E /~prn m dP(j\x m ) (i9) 

y m'*m 

after using Eq (17) 

Now for memoryless channels we can show that (see 
Appendix) 

p E {m^m) < | e ] \fjjy)pjy) | mm 



In the second inequality above, the summation is over all 
tri 

We can bound the probability of bit error for the input 
message m, p(b\m), by weighting each of the terms in the sum 
in Eq (21) by the number of bit errors i mm < occurring in 
choosing that incorrect path m , 1 e , 



Accounting for all possible transmitted messages m, the overall 
BEP is given by 



(20a) p ft (MAP) = p(b\m)p m 

m 


where w m m > is the Hamming distance between messages m 
and m , and p x (y), p 0 (v) are the probabilities of receiving a 
given value of a bit y given that the corresponding transmitted 
bit is a 1 or 0, respectively Here 


< y y* & ~r 1 ,z w,nm ' (23) 

i—l V m m mm 
m m Vm 


z = E n/PiWp 0 W ( 20b ) 

y 

which, for the special case of a BSC, reduces to 


Z = Z bsc =V r 4P e ( 1 -P e ) (20c) 

We now apply these equations to the calculation of upper 
bounds on event and bit error probabilities We know that a 
necessary condition for an event error to begin occurring at 
node ] is that an incorrect path diverging from the correct one 
at node / accumulates higher total metric than the correct one 
over the unmerged segment If we denote by T the set of all 
such incorrect paths when the input message is m , using the 
union bound as in Eq (18), we get, for the probability of an 
event error occurnng at node / with m input, 


Equations (22) and (23) can, in theory, be evaluated for a 
specific code, but, except for certain cases, as in those of 
Subsection III-F, are very cumbersome in practice 

We cannot make further simplifications in Eqs (22) and 
(23), directly at least, as we can in the usual ML decoding case, 
because each sum in the bound in Eq (22) depends on the 
particular message m and the incorrect message m For certain 
specific messages m it is possible to extend the generating 
function method (as we do in Subsection III-F) to obtain 
elegant expressions for p(b\m), but that is not sufficient to 
enable calculation of the overall BEP p> 6 (MAP), since all mes- 
sages do not have the same error probabilities 

C. Comparison of MAP and ML Decoding 
Error Bounds 

If the usual Viterbi or ML decoding were used for the same 
system, we would have, for the bit error probability given that 
message m is sent, 
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p(b\m) < E zWm,n 

m'^m 


l 

mm 


(24) 


where vv mm ' and i mm ‘ are the same as in Eq (22) The overall 
BEP would be 


Pb (m = e p(^ m )p m 

m 


< 


— > w 

E p i z 1 

z -j r m mm 
m m'j=m 


(25) 


Now Eq (23) and Eq (25) both involve a double summation 
over all possible pairs of messages m, m We compare these 
term by term Consider a pair of messages m, and m 2 with 
probabilities p m { , p m2 Their contributions to the sums in 
Eqs (23) and (25) are 


W 

m .m ^ 

C = Z 1 2 1 

MAP m i m 2 


P P 
r m ,m 




D. Approximate Magnitude of Improvement 
Expected From Upper Bound Expressions 

From Eqs (23) and (25) we may obtain an order-of- 
magmtude estimate of the improvement in BEP using MAP 
decoding We assume a low enough BEP so that the actual BEP 
is replaceable by the upper bounds in Eqs (23) and (25) Since 
the BEP is low, we must have Z small, so that the terms in Eqs 
(23) and (25) with factors Z Wmm ' , where is greater 

than the free distance d f of the code, are negligible compared 
to the terms with factors Z d f Thus, we may retain only these 
terms and have 


P b ( ML) 
P b ( MAP) 


E P Z df > ' 

mm'er 


E frr,z d fi 

W m r m mm 


£ 


PrJr 


£ 


V PmPm 


mm' 


(30) 


= ^sJ~P 


p 

m m ^ 


(26) 


and 


\ w mm 

I Z 1 2 i 

2/ m i m 2 

Comparing Eqs (26) and (27), we see that 
c 

'"MAP '“ML 

for all p m , , p m 2 , with equality holding for P ml =P m2 
we get 

(MAP) error upper bound < (ML) error upper bound 



(27) 


(28) 


Thus, 


(29) 


Since the upper bounds in Eqs (23) and (25) are asymptoti- 
cally tight, we expect from Eq (29) that, for low BEP where 
the bounds are close to the actual error probabilities, the new 
metric is better than, or as good as, the conventional one 
Also, on general principles, we know that MAP decoding is 
always at least as good as ML decoding 


where I’ is the set of path pairs m, m with Hamming distance 
equal to d f between them 

Further simplification is possible only if we know the 
probabilities of the message pairs that are at the free distance 
from each other and this is specific to the code and the source 

E. Further Approximation for First-Order 
Markov Source 

In this section we make an order-of magnitude estimate of 
the improvement in BEP for the specific case of a first -order 
Markov source with low transition probability p s For such a 
source with low p s , there is a lot of redundancy m the message 
output, and we may expect the MAP method to be most 
useful in this case In making our estimate of the improvement 
in BEP in this section, we make the following assumptions 

(1) We assume that the most probable message transmitted, 
m , , has an a priori probability so close to 1 that the 
overall BEP p 6 (MAP) may be closely approximated by 
the BEPgivenWj, p(b\m l ) 

(2) Given that an error is made in decoding this message 
fflp we assume that the BEP given m t can be 
calculated by considermg only those possible erro- 
neously decoded messages that are closest m distance 
to m 1 In doing this, we are agam using the arguments 
leading to Eq (30) in Subsection I1I-D 
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(3) Among these closest-distance erroneous messages, we 
further restrict our attention to the one that has the 
highest a pnon probability There may be more than 
one such message possible but we ignore this small 
factor m gettmg our order-of-magnitude estimate 

Thus, assumption (1) enables us to consider only one 
transmitted message m l and assumptions (2) and (3) enable us 
to narrow down the set of possible erroneously decoded 
messages (given that m { is sent) to only one message, which 
we call m 2 

Smce the first-order Markov source has a low transition 
probability p s , the most probable messages it puts out will 
consist of long strings of zeros and long strings of ones with 
occasional (with probability p s ) transitions from one type of 
strmg to the other We now make the following assumption 

(4) We assume that the most probable message m l has long 
enough runs of zeros and ones that we can replace it, 
for the purpose of calculating BEPs, by one with no 
transitions at all, say, the all-zero message 

Now the messages closest to m ] will differ from it in only a 
small number of bits (like 1 or 2 depending on the specific 
code) We can see that the most probable of these, le , m 2 , 
will have at least 2 transitions m it For instance, the relevant 
portions of m, and m 2 may be -000- and -010-, or -0000- and 
-0110- Thus, for the specific code discussed m Section IV, if 
an error occurs in decodmg a portion -000- of a message, and if 
the erroneously decoded message is at the free distance from 
the correct one, it occurs because -000- is decoded as -010- 


For instance, if = 10 -3 , we have from Eq (33), 


old BEP 
new BEP 


499.5 


( 34 ) 


Thus, there is about a factor of 500 improvement m BEP in 
this case usmg the new metric 

F. Generating Function Method for Low p s Case 

In this section, we show how to calculate the improvement 
in BEP more exactly than we did in Subsection III-E The 
method to be described holds for all p s if the only message 
considered is the all-zeros or the all-ones message and also 
possibly for certam specific messages dependmg on the code 
For low p s , the most probable messages consist of strings of 
zeros and ones with occasional transitions between strings so 
that this analysis holds for most of the probable message 
sequences 

We have for 2 messages m, m , 


P E ( m 



( 20 ) 


which differs from the corresponding equation for the 
conventional decoder only in the factor P m '/P m Her e,p m 
is the probability of the transmitted message sequence m 
consideration and p' m is the probability of the erroneously 
decoded message m We can express p m and p' m using Eq (5) 
for the special case of the first-order Markov source as 


We may now calculate the improvement in BEP considermg 
only the messages m l and m 2 For symmetry, we mclude in 
p 6 (MAP), both p{b\m x ) and p(b\m 2 ), where the conditional 
bit error probability for each message is calculated assuming 
the other one as being the erroneously decoded message We 
do this usmg the method of Subsection III-D, and get 


P b (ML) 
P„ (MAP) 


P m 1 2 

^■\fP ml P m 2 


Defmmg 


% = 1 ' P s 


(31) 


(32) 


where 


PmtlPm(B-l) Pm 1 

(35) 

Pm\ 

(36) 

P { U mi^ U m(i- 1)) 

(35a) 

= '0-1)1 

(36a) 


So we have 


this becomes 


p 6 (ml) q s + p s 


p (MAP) 2 p 2 q 


( 33 ) 


P E (m -> m') < Z 


Pm ' R Pm ' 1 

m d m 1 


mB 


m 1 


(37) 


We can then rewrite Eqs (21) and (22) as 
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pfe\m) < Y, Z ' Vmm 

m i^m 



P m’ 1 
P m 1 

(21)' 


previous paragraph for the higher order Markov cases) and the 
transfer function T(D,M,f) (to take the place of T m (D,f) of 
the previous paragraph) We can then see that the event error 
probability at any node, say,/, is given by 


p(b\m) < 


E 

m 7 ^ m 


vv , 

i ,z mm 

mm 



p t (e \m) < T(D,M) 


M 



D 


Z 


and the BEP given m is given by 


(38) 


Now each of the factors \J P m ',/P m , in Eq (37) corresponds 
to a particular branch in the trellis If the input message is such 
(e g , the all-ones or the all-zeros message) that all possible 
mcorrect paths can be accounted for by tracing varying 
numbers of loops on a diagram obtained by openmg up the 
state diagram at some node, then we know that all the paths at 
various distances from the correct one can be expressed by 
terms of the series expansion of the generating function T(D ) 
(§4 3 of Ref 2) Here, D has the usual meaning that its 
exponent for any branch represents the distance of that 
branch from the corresponding code subvector of the input 
message In order to weight each of these terms by the factor 
s/Pm'JPmr a ^ we need t0 do ls to associate with each 
branch of the opened-up state diagram the appropriate factor 
V Pm’JPmi an< * t ^ en evaluate the resulting generating func- 
tion which we denote by T m (D) We then set D = Z in the 
expression T m {D ) to evaluate the bound in Eq (21)' 
Simdarly, to get the bound in Eq (22)', we obtain the 
generating function T m {D,I), where I has the usual meaning 
that 7=1 for a given branch if its bit differs from the 
corresponding bit in the input message and 7 = 0 otherwise, 
and we set D = 2, 7 = 1 in bT m /bI (D,I) to evaluate the bound 
in Eq (22)' This generating function method allows us to take 
into account all possible erroneous paths and their error 
contributions rather than just one as in Subsection III-E 

In the special case of the first-order Markov source, we get 
an even simpler expression for event and bit error probability 
bounds Each of the p ml and P m \ is either p s or q s in this 
case So each term in the sums in Eqs (21)' and (22)' has a 
factor of M k evaluated at M = \JpJq s , where k is an integer 
This integer can be determined as follows We first associate 
with each branch of the opened-up state diagram either M or 1 
or , depending on whether p m -Jp mt = pjq s or 1 or qjp s 
for the particular branch Thus, for the all-zeros or the all-ones 
message, we write M for branches corresponding to a transition 
m the message and 1 for branches with no transition, since m 
these cases, p ml is always equal to q s and p m > t can be either p s 
or q s We then evaluate the transfer function T(D,M) (this 
takes the place of the transfer function T m (D ) we had in the 


p(b\m) < 


i = i 


(39) 


M 



It is not easy to generalize these results to get the overall 
BEP considering all possible messages, because the functions 
T m (D), T m (D,f) here are message-dependent unlike the usual 
Viterbi decoding case Thus, for non-repetitive messages, a 
general expression like those in Eqs (38) and (39) that 
identifies all other possible paths and their error contributions 
are not easy to obtain A slight generalization of one method is 
possible for repetitive messages like, say, 1000, 1000, 
1000 for the code considered in the next section, by open- 
ing up, not just one node, but a path in the state diagram cor- 
responding to the message unit repeated But even in this case, 
it is necessary to consider each possible erroneously decoded 
message individually and obtain its contribution to the BEP 


As we have stated before, if p s is low, then just considering 
the all-zeros or all-ones message should be good enough to 
obtain an estimate of the overall BEP 


IV. Illustration With a Specific Code 

A. The Code 

We assume a BSC, a first -order Markov source with 
transition probability p s and a specific K = 3,b = l ,n = 2 (rate 
1/2) code The encoder is shown in Fig 1, and the state 
diagram is shown in Fig 2 Bits come m at the left in Fig 1 
and the states (written in boxes at the nodes in Fig 2) are 
represented by the rightmost 2 bits of the encoder The dotted 
and sohd lines in Fig 2 correspond to a 1 and a 0 input bit, 
respectively 
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B. Calculations 

The state diagram with the 00 node opened up and the 
branches marked with the appropriate powers of D,I,M are 
shown in Fig 3 This diagram thus corresponds to the all-zeros 
input message 


The expansion of Eq (41b) shows that the path with 
minimum distance (=5) has 2 transitions and one bit error, 
correspondmg to the term D S M 2 , of the two paths with 
distance 6, one has 2 transitions and the other 4, and each of 
these has 2 bit errors correspondmg to the 2 D s • M 2 • D(l + 
M 2 ) term, and so on 


The generating function is 


D S IM 2 

T(D,I,M ) = 

1 -DI- DIM 2 


From the expansions of Eqs (41b) and (42b), we see that if 
M 2 <l,ie , 


(40) 


Pj <0 5 


(43) 


so that the BEP given m = 00 is 


P(b\ m )\ MAP <jj 


i = i 


M 



D 


Z 


the bound of Eq (41b) is smaller than that of Eq (42b) by a 
factor of more then M 2 , since each term in the brackets in Eq 
(41b) is smaller than the correspondmg one m brackets in Eq 
(42b) 

In our example, we assume M 2 < 1 If we assume the 
bounds in Eqs (41) and (42) to be close to the actual BEP 
because of our p s « 1 assumption, there is an improvement 
m BEP by more than a factor of M 2 = pjq s Thus, for p s = 
0 1 , we have that the new BEP is lower than the old BEP by at 
least a factor of 1/9 


D S M 2 

(1 - D- DM 2 ) 2 


\M 



(41a) 


For the all-ones input message, the 1 1 node can be opened 
up If the details are worked out, we obtain, as expected, 
exactly the same expression as m Eq (41) 


C. The Simulations 


\D = Z 


= L?M 2 {1 +2D(1 +M 2 ) + 3D 2 (1 +M 2 ) 2 


+ 4£> 3 (1 +M 2 ) + } 


M 


D ~ Z 



(41b) 

For the case of the conventional ML decoder, we would 
obtam a similar expression with M = 1 , i e , 


P(.b\m) ML < 




(1 - 2D) 2 


(42a) 


D-Z 


< D* (1 + 2D • 2 + 3D 2 • 2 2 + 4D 3 • 2 3 + } 


D=Z 

(42b) 


We simulated the encoder and decoder for the code 
described above and ran several messages of 298-bit block 
lengths for various sets of parameters p s and p e , with both p s , 
p e < 0 5 Initially, we checked the program with just a block 
of zeros or just a block of ones as mput and, later, generated 
the mput messages themselves within the program by including 
a simulation of the Markov source In this case the mput 
messages are realistic ones and were found to consist of strings 
of zeros and ones as expected The simulations show that the 
BEP with the new metric is always lower (for p s , p e < 0 5) 
than with the conventional metric But a larger number of 
simulations would be needed to obtam the exact amount of 
improvement, especially for the cases with low BEP However, 
the results shown in Table 2 clearly demonstrate the improve- 
ment We also see that for low enough p e , the improvement m 
BEP is by a factor > pjq s , as predicted by Eqs (41) and (42) 

Another fact that was mdicated by the simulations (see 
Table 1), but which again needs further verification, is that the 
decoder with the new metric needs a trellis truncation depth 
of at least 20 constramt lengths for negligible truncation error, 
as opposed to the 5 constramt lengths required for the old 
decoder (Ref 3) The truncation scheme we are talking about 
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here is the simple one described in §4 7 of Ref 2, in which, as 
each set of b bits enters the registers of each state, the b bits 
which entered 20K branches earlier are removed, after the 
decoder has made a final decision on these bits by settmg them 
equal to the appropriate survivor bits of an arbitrary state In 
obtammg the results described m Table 2, we used a 40K 
truncation length rather than 20K, to be on the safe side 

D. Channel Capacity Limitations 

From the results of the simulations (Table 2) we see that 
for given p s , if the probability of channel transition p e is high, 
the BEP is high for both the new and the old decoding 
algorithm, and the ratio of the new BEP to the old BEP 
increases with p e Similarly, for a given p e , this ratio increases 
as p s increases For the given source and channel, we can 
calculate the source entropy and the channel capacity, which 
are given, respectively, by 

H(P S ) = ~P S >og 2 p s ~ (1 - p s ) log 2 (1 - pj (44) 
and 

C(p e ) = 1 - H(p e ) 


Some sets of parameters p e and p s satisfying this relationship 
with the equality sign are marked with an * in Table 2 Figure 
5 shows a plot of old BEP, new BEP, and the ratio of new BEP 
to old BEP obtained from simulations with the Markov source 
model included in the program Two sets of curves are shown 
one set has p s = 0 1 and p e varying, the other has p e = 0 1 and 
p s varying Since an insufficient number of simulations were 
made, however, the shape of the curve is not very reliable But 
it does show an increase of the ratio of new BEP to old BEP as 
we approach and exceed channel capacity 

E. Comparison of Calculations and Simulations 

None of the simulations were made for the very low values 
of p e (like 10" 5 or 10~ 6 ) that would be required for the upper 
bounds of Eqs (41) and (42) to be approached, since too 
many blocks of input would be required Even for the case of 
p s = 0 1 and p e = 0 03, which represents about the lowest data 
rate to channel capacity ratio m the simulations, there were 0 
bit errors with the new decoder in the 10 blocks of 298 bits 
used and only 1 bit error with the conventional decoder in the 
same 10 blocks, corresponding to an old BEP of 0 000335, 
and an even lower new BEP Here the bounds of Eqs (41) and 
(42) give, respectively, 


= \+p e log 2 p e + (1 - p e ) log 2 (1 - p e ) (45) 

Figure 4 shows the nature of these two curves, on the same 
plot for convenience 

Since the channel capacity represents the upper limit on the 
reliable communication rate, we should expect that no scheme 
would allow any substantial reduction in BEP above channel 
capacity Thus, if we set 

H(p s ) < C(p e ) (46) 


P 6 (ML) < 0 046 
p ft (MAP) < 0 001 


(49) 


which are both much higher than the actual BEP obtained by 
simulation Hence direct verification of the bounds of Eqs 
(41) and (42) has not been performed in the simulations, but 
the bounds are probably very useful in calculating BEPs at 
rates well below capacity where too many simulations would 
be needed if computer calculation of BEPs were to be 
performed 


we can get the outer limits on the regions of p e and p s within 
which the new scheme can be expected to be useful 

We need, from Eqs (44) through (46), 

" P e log 2 P e - (1 " P e ) log 2 (1 - p e ) - p s log 2 p s 

-(l-p j )log 2 (l-p j )<l (47) 

or 

p (i-p ) p (i-pj 

(P e ) e V~P e ) e (p s ) s (l-p s ) *> 0 5 (48) 


V. Conclusions 

We have shown that a simple modification of the metric 
used in the Viterbi decoding algorithm achieves MAP decodmg 
for sources with memory It can cause a noticeable reduction 
in BEP for sources with strong correlation between bits (l e , 
low entropy) Analysis methods to obtain upper bounds on 
the BEP obtainable with the new metric have been given The 
simulations performed for a BSC and a First-order Markov 
source verify that an improvement occurs, but to obtain the 
exact ratio of improvement, more simulations are needed We 
can conclude, however, that the ratio of new BEP to old BEP 
for a first-order Markov source and a BSC with very low p e is 
at least pjq s and that the extent of improvement reduces as 
channel capacity is approached and exceeded Simulations also 
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indicate that a trellis truncation depth of 20 constraint lengths 
is probably needed for reliable results with the new metric 

The modification of the metric required for the case of an 
arbitrary memoryless Gaussian channel is given The entire 
analysis given holds for such a Gaussian channel, but the 
simulations of the decoder have been performed only for the 
case of the BSC Except in cases where the first -order Markov 


source is specified, the analysis holds for any Markov source of 
arbitrary order More work is needed to extend the results to 
more realistic sources like video data sources Ultimately, the 
aim is to devise a reasonably simple modification of the 
Viterbi decoding algorithm that can make use of the correla- 
tion between neighboring pixels m video data of scenes to 
enable more efficient channel use than is possible with the 
normal Viterbi decoder 
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Table 1 Simulations described in Subsection IV-C to estimate effect of trellis truncation (all runs here have p, = 0 1, p e = 0 4) 


Program 

Truncation length, 
units 

Number of 298-bit 
blocks 

Conventional 

BEP 

New 

BEP 

Remarks 

01 

10K 

1 

0 4866 

0 5336 

Programs 01 through 04 show 

02 

10K 

1 

0 4832 

0 5537 

that 20K truncation length 
is needed 

03 

20K 

1 

0 4899 

0 4698 


04 

30K 

1 

0 4899 

0 4698 

Same sets of random numbers 
used for the source and 

05 

30K 

5 

0 5080 

04170 

channel in 04 as in 03, other- 

06 

4 OK 

1 

0 4765 

0 2416 

wise different sets used for 

07 

40K 

1 

0 5101 

0 4060 

different blocks of data for 
both source and channel 

08 

40K 

5 

0 5000 

0 3698 

simulations 


Table 2 Simulations described in Subsection IV-C (all with truncation length 40K) 


Program 



Number of 298-bit 
blocks 

Conventional 

BEP 

New 

BEP 

Ratio 
New BEP 

Remarks 

P s 

P e 

Conventional BEP 


1 

0 1 

04 

5 

0 5000 

0 3698 

0 7396 

All runs have different 

2 

0 1 

03 

3 

0 4508 

0 2573 

0 5708 

sets of random numbers 
for the source and 

3 

0 1 

0 2 

3 

0 3367 

0 09396 

0 2791 

different sets for the 

4* 

0 1 

0 12 

10 

0 1117 

0 02886 

0 2583 

channel 

5 

0 1 

0 08 

12 

0 02041 

0 005593 

0 2740 


6 

0 1 

0 05 

15 

0 006714 

0 0004474 

0 0667 


7 

0 1 

0 03 

10 

0 000335 

0 

0 


8* 

0 12 

0 1 

10 

0 06146 

0 013194 

0 2147 


9 

0 05 

0 1 

5 

0 07651 

0 01342 

0 1754 


10 

0 03 

0 1 

10 

0 05570 

0 002013 

0 03614 


11* 

0 2 

0 04 

10 

0 003691 

0 001678 

0 4546 


% 

Means p ^ 

p for that program define operation at channel capacity 
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Fig 4 Source entropy vs p, and channel capacity vs p e 


Fig 2 State diagram for encoder of Fig 1 
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Fig 5 BEP with new and old metrics and their ratio vs p , with 

Fig. 3 State diagram labeled with distance, bit errors, and branch P» constant and vs p, with p # constant (from Table 2) 

probability ratios relative to all-zeros input 










Appendix 

Simplifying the Extended Bhattacharya Bound for Memoryless Channel 


We show here the derivation of Eq (20) from Eq (17) for a 
memoryless channel The method used is similar to that in 
§2 9 of Ref 2 We have, from Eq (17), since p(y\x m ) and 
p(y\x m <) can be factored for a memoryless channel, 



(All) 


where u m/ , u m - l are the input bits corresponding to the code 
subvectors x m) , x m ‘ r the sum in Eq (17) runs over all possible 
y vectors consisting of B ^-dimensional subvectors, and the 
sum over y in Eq (A 1 - 1 ) runs over the two possibilities 0 and 
1 for each of these nB components 

We have 


I! Z \/pOT^ / )p(yl« mV ) 

J u i -u y 

' mi mi / 


(here, u mj = logical complement of « m; ) This bound can be 
written as 



(Al-2) 


since each sum in the first product in the previous equation is 
equal to 1 

Thus, if the messages m, m differ in w m m - bits, we get 
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